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Preface 


The  SPC  (Signal  Processing  &  Conununicadons)  software 
package  is  the  result  of  work  originated  in  Spring  1993  at  the  Naval 
Postgraduate  School  as  a  class  project  in  speech  processing.  The 
project  developed  an  interactive,  user-friendly,  tool  to  analyze 
speech  signals  using  MATLAB.  Then,  as  with  any  project  left  with¬ 
out  a  fixed  deadline,  it  grew  to  include  nunre  sophisticated  analysis 
tools  such  as  various  filtering  techniques.  Autoregressive  (AR)  and 
linear  Auto  Regressive-Moving  Average  (ARMA)  modeling  meth¬ 
ods.  Finally,  it  was  expanded  to  include  basic  communication  tools. 

The  final  product  is  a  software  package  designed  to  provide 
the  user  with  a  series  of  data  manipulation  tools  which  use  MAT- 
LAB  V.4  graphical  user  interface  controls.  SPC  can  be  used  in  the 
classroom  to  illustrate  and  to  reinftxce  basic  concepts  in  digital  sig¬ 
nal  processing  and  communications.  It  frees  the  user  from  having  to 
write  and  debug  his/her  own  code  and  gives  him/her  more  time  to 
understand  the  advantages  and  drawbacks  of  each  technique 
included  in  the  package.  It  can  also  be  used  as  a  basic  analysis  and 
modeling  tool  for  research  in  Signal  Processing.  SPC  is  well  suited 
to  reinforce  basic  concepts  presented  in  the  following  courses 
offered  at  the  Naval  Postgraduate  School; 

EC4410:  Speech  Processing, 

EC4420:  Modem  Spectral  Estimation, 

EC3420:  Statistical  Digital  Signal  Processing, 

EC34(X):  Digital  Signal  Ifrocessing, 

EC2S(X);  Communication  Theory. 

We  hope  that  users  will  find  this  package  useful,  and  we  wel¬ 
come  any  comments  and  suggestions  regarding  this  software  at 
browndw(g>ece.nps.navy.mil  (until  6/94),  or  fargues@ece.nps.na- 
vy.mU. 
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Commands 


Signal  Generation: 
antpodal 
unipolar 
sawwave 

sqwave 

triwave 

coswave 

sinwave 

bfsk 

bfskmsg 

bpsk 

bpskmsg 

ook 

ookmsg 

setsnr 

setsnrbw 


Coding; 

Its 

str2masc 

masc2str 

Modulation: 

dsbsc 

dsblc 

modindex 

envelope 

Filtering; 

avsmooth 

mdsmooth 

Data  Files: 

IdSbit 

Idltibit 

ld82bit 

loadsspi 


Generate  a  baseband,  antipodal  [-l.ll  signal. 

Generate  a  baseband,  unipolar  [0,1]  signal. 

Generate  a  baseband,  unipolar  [0,1]  or  antipodal  [-1.1]  sawtooth  wave.  This 
function  can  also  be  used  to  generate  a  ramp. 

Generate  a  baseband,  unipolar  [0.1]  or  antipo^  ['1.1]  square  wave. 
Generate  a  baseband,  unipolar  [0.1]  or  antipodal  [-1,1]  triangular  wave. 
Generate  a  sampled  cosine  wave. 

Generate  a  sampled  sine  wave. 

Generate  a  frequency-shift  keyed,  bandpass  signal  with  a  random  message. 
Generate  a  frequency-shift  keyed,  bandpass  signal  with  a  specific  message. 
Generate  a  phase-shift  keyed,  bandpass  signal  with  a  random  message. 
Generate  a  phase-shift  keyed,  bandpass  signal  with  a  specific  message. 
Generate  an  on-off  keyed,  bandpass  signal  with  a  random  message. 

Generate  an  on-off  keyed,  bandpass  signal  with  a  specific  message. 

Mix  two  signals  such  that  the  ouq)ut  signal’s  total  bandwidth  signal-to-noise 
ratio  is  set  to  the  specified  value. 

Mix  two  signals  such  that  the  output  signal’s  noise-equivalent  (in-band)  sig- 
nal-to-noise  ratio  is  set  to  the  specified  value. 


Generate  a  maximal-length,  linear  recursive  sequence. 

Convert  a  string  into  its  binary  ASCII  representation  as  generated  by  a 
modem. 

Convert  a  vector  of  1  ’s  and  O’s  representing  ASCII  into  a  string. 


Generate  a  double-sideband,  suppressed-carrier  amplitude  modulated  signal. 
Generate  a  double-sideband,  large-carrier  amplitude  modulated  signal. 
Returns  the  modulation  index  of  an  amplitude  modulated  signal. 

Returns  the  envelope  of  a  amplitude  modulated  signal. 


Smooth  a  curve  using  an  avera^  smoothing  filter. 
Smooth  a  curve  using  a  median  smoothing  filter. 


Load  a  file  of  samples  stcxed  as  8-bit  signed  integers. 

Dito  for  16-bit  signed  integers. 

Dito  for  32-bit  signed  integers. 

Load  a  datafile  store  as  ASCII  in  the  fonnat  specified  by  the  Cyclic  Spectral 
Analysis  Software  Packa^  fi'om  Statistical  Signal  Processing,  Inc. 
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Comnunris 

loadvoc  Load  vector  from  a  Creative  Labs  Soundblaster  Voice  file, 
savevoc  Save  vector  to  a  Creative  Labs  Soundblaster  Vsioe  file. 

AR/ARMA  Modeling: 

ar.corr  Compute  parameters  for  an  AR  model  using  the  autocorrelation  method. 

ar_covar  Compute  parameters  for  an  AR  model  using  the  covariance  method. 

ar_mdcov  Compute  parameters  for  an  AR  model  using  the  modified  covariance 
method. 

ar_burg  Compute  parameters  for  an  AR  model  using  the  Burg  method. 
ar_prony  Compute  parameters  for  an  ARMA  model  using  the  Prony  method, 
ar.durbn  Compute  parameters  for  an  ARMA  model  using  the  Durbin  method, 
ar.shank  Compute  parameters  for  an  ARMA  model  using  the  Shank  method, 
ar.levin  Compute  parameters  for  an  AR  model  using  the  Levinson  recursimi. 

Least  Squares  Filtering: 

Is.whopf  Compute  the  filter  coefficients  for  a  least-squares  (^timal  filter  using  the 
Wiener-Hq)f  equation. 

ls_svd  Compute  the  filter  coefficients  for  a  least-squares  optimal  filter  using  singular 
value  decomposition. 

Linear  Systems: 

minphase  Return  a  polynomial  that  is  in  minimum-phase  form, 
maxphase  Return  a  polynomial  that  is  in  maximum-phase  form, 
nonnaleq  Solves  a  system  of  Normal  equations  used  in  linear  predictive  filtering. 

Graphic  Display: 

plottime  Plot  the  time-domain  display  of  a  signal. 

Iperigrm  Display  die  periodogram  of  a  signal  on  a  logarithmic  scale, 

wperigrm  Display  the  periodogram  of  a  signal  on  a  linear  scale. 

Speech  Processing: 

sp.steng  Compute  the  short-time  energy. 

sp_stmag  Compute  the  short-time  magnitude, 

sp.stzcr  Compute  the  short-time  zero-crossing  rate. 

Dodads: 

vectedit  Interactive  cut  and  paste  vector  editing, 
vectfilt  Interactive  filter  design  and  filtering, 

vectarma  Interactive  AR  and  ARMA  modeling, 
vecttime  Interactive  speech  time-domain  analysis, 
sanalyzr  Spectrum  analyzer. 

grafiltr  Grtqihical  filter  design  through  individual  (complex-conjugate  pair)  place¬ 
ment  of  poles  and  zeros  on  the  pole-zero  plot  of  the  filter  transfer  frmction. 
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INTRODUCTION 


The  SPC  Toolbox  provides  a  number  of  cmnmands  implementing  a  variety  of  dig¬ 
ital  signal  processing  and  communication  techniques.  All  of  the  routines  are  imple¬ 
mented  using  Matlab  (version  4.0)  M-files.  Most  of  the  commands  are  stand-alcme 
or  rely  on  other  SPC  Toolbox  routines.  Some  commands  however,  rely  on  com¬ 
mands  from  the  Matlab  Signal  Processing  Toolbox  available  from  llie  Math- 
Works.  The  functions  that  rely  on  the  Signal  Processing  Toolbox  are  routines  that 
implement  FIR  windows,  analog  IIR  filter  prototypes,  or  some  spectral  estimates. 

The  SPC  Toolbox  commands  can  be  grouped  into  the  following  areas; 

•  Sampled  signal  generation. 

•  Bandpass  signal  generation  and  modeling. 

•  Information  signal  coding. 

•  Communications  systems. 

•  AR  and  ARMA  modeling. 

•  Least  squares  filtering. 

•  Linear  systems. 

•  Speech  signal  processing. 

•  Signal  display  functions. 

•  Formatted  signal  data  file  loading  and  saving. 

Additionally,  a  number  of  graphics-based  tools  (referred  to  as  “dodads”  in  this 
document  for  lack  of  a  better  name)  provide  an  interactive,  graphical  environment 
ofrering  relief  from  the  Matlab  command  line  prompt.  Some  dodads  essentially  act 
as  “front-ends”  from  which  other  commands  are  called  (i.e.  the  vector  filtering 
dodad).  Others  provide  interactive  model'uig  facilities  (i.e.  AR/ARMA  modeling 
dodad).  Still,  others,  provide  a  graphical  environment  to  perform  tasks  not  readily 
available  or  easily  performed  using  Matlab  commands  (i.e.  visual  vector  editing^ 
The  following  dodads  are  provided  in  the  SPC  Toolbox: 

•  Vector  Editor  Dodad  (visual  “cut  and  paste”  vector  editor,  veciedit). 

•  Vector  Filtering  Dodad  (interactive  filtering,  vectfit). 

•  AR/ARMA  Modeling  Dodad  (interactive  modeling,  vectanna). 

•  Speech  Tune-Domain  Analysis  Dodad. 

•  Spectnun  Analyzer  Dodad  (spectral  measurement). 

SPC  Toolbox  Conventions 

The  following  conventions  are  followed  by  all  commands  in  the  SPC  Toolbox: 

•  All  commands  have  been  tested  in  Matlab  version  4.0,  running  on  Sun 
workstations  or  on  PC’s  under  Microsoft  Windows.  Although  untested, 
SPC  Toolbox  commands  should  work  with  no  problems  on  other  platforms 
running  Matlab  version  4.0. 


SPC  Toolbox  User’s  Guide 


1 


INTRODUCTION 


•  Excq)t  fm  dodads,  no  commands  are  interactive  (i.e.  request  iiq)ut  from 
the  user  via  prompts). 

•  All  commands  taking  vectors  as  input  aiguments  accept  eitho-  Nxl  or 
IxN  vectors. 

•  Vectors  are  returned  as  Nxl  vectors  (except  functions  returning  vectors 
rqpresenting  polynomials  which  return  IxN  vecttxs  following  the  Matlab 
polynomial  convention). 

•  Aiguments  returned  after  an  eiTOT  occurs  are  returned  as  empty  vectixs. 
This  allows  errcx  checking  by  user  implemented  routines  with  the  follow¬ 
ing  code: 

y  -  bp6k(S0.10240: 

« isempty(y). 

errorCAn  error  occurred  generating  a  BPSK  signal.”); 
end; 

•  The  default  sampling  rate  used  is  that  of  the  Sun  workstation  speaker  out¬ 
put  equal  to  8192  Hz.  Hiis  allows  signals  to  be  heard  just  as  if  they  were 
being  heard  from  a  radio  speaker.  PC’s  with  soundcards  also  have  this  fea¬ 
ture  available,  except  that  the  ouq>ut  sampling  frequency  is  not  limited  to 
8192  Hz. 

•  On-line  help  is  available  by  typing  “help  command.name’’  at  the  Matlab 
prompt.  A  description  of  all  SPC  Toolbox  commands  can  be  viewed  by 
typing  “help  spdools”. 

•  In  the  documentation,  commands  are  printed  in  a  sans-serif  font  and  vari¬ 
able  names  are  printed  in  italics.  Commands  in  the  table  of  contents  are 
shown  in  capital  letters.  AH  examples  should  work  as  shown  if  the  toolbox 
is  properly  installed. 

Customizing 

Most  colors  of  objects  within  dodads  can  be  customized  to  meet  individual  prefer¬ 
ences  by  editing  the  file  spoolors.m.  On  workstations,  individual  users  can  copy  the 
spoolors.m  file  from  the  spctools  directory  and  modify  it  to  meet  their  nee^.  To 
make  individual  copies  of  spoolors.m  take  precedence  over  the  default  spoolors.m, 
the  only  requirement  is  that  the  individual’s  copy  of  spoolors.m  is  reachable  in  the 
matlabpath  before  the  spctools  directory.  This  can  be  accomplished  by  keeping  a 
copy  of  spoolors.m  in  the  current  directory  or  by  prepending  the  madabpath  with  the 
directory  containing  the  individual’s  copy  of  spoolors.m.  See  the  Matlab  User’s 
Manual  for  information  on  the  matlabpaih  variable. 

In  addition  to  coIots,  the  default  size  and  location  of  dodads  may  be  specified  in 
the  spoolors.ni  file.  Simply  “uncomment”  the  lines  according  to  the  instructions  in 
the  spoolors.m  file  and  replace  the  parameters  for  the  location  and  size  with  the 
user’s  preference.  The  default  screen  sizes  are  nine  tenth  of  the  screen  for  screens 
less  than  800  pixels  wide,  eight  tenths  of  the  screen  for  screens  from  800  to  less 
than  1024  pixels  wide,  seven  tenths  of  the  screen  for  screens  1024  pixels  wide  and 
five  tenths  of  the  screen  for  screens  greater  than  1024  pixels  wide.  The  default  size 
and  location  for  child  graphic  windows  cannot  be  specified. 
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The  file  s()lMndw.m  contains  the  default  colors  for  a  black  and  white  screen.  Using 
these  colors  for  the  default  requires  executing  the  Matlab  wtvtebg  command  before 
executing  any  dodad  comma^.  To  use  this  file,  copy  or  rename  to  file  to  spool- 
ors.m  and  follow  the  instructions  above. 

Error  Messages 

A  considerable  amount  of  error  checking  is  performed  in  the  SPC  Toolbox  routines 
in  an  attempt  to  avoid  disasters.  All  error  messages  generated  by  SPC  Toolbox  rou¬ 
tines  have  the  fcxmat: 

hmdion.name:  error  message 

If  you  receive  an  error  message  that  does  not  have  the  name  of  a  SPC  Toolbox  rou¬ 
tine,  the  error  message  was  generated  by  Matlab  or  a  command  from  some  other 
toolbox.  Currently,  the  SPC  Toolbox  relies  only  on  version  4.0  Matlab  commands 
and  a  few  commands  from  the  Matlab  Signal  Processing  Toolbox. 

If  you  receive  an  error  while  executing  a  SPC  Toolbox  command,  please  send 
email  to  “browndw@ece.nps.navy.mil”  (good  until  June  94)  or  “fargues@ece.tq>- 
s.navy.mil”.  Include  a  description  of  the  operations  in  progress  when  the  erTtM*  was 
received  and  a  copy  of  the  error  message  itself.  If  possible,  cut  the  command  and 
error  message  right  from  the  shell  used  by  Matlab  and  paste  it  into  the  email  text 
using  the  X  Windows  (or  MS  Windows)  cut-and-paste  facilities. 

Troubleshooting  Hints 

There  are  some  limitations  in  using  Matlab  as  an  interactive,  graphical  tool  due  to 
limitations  of  the  Matlab  progranuning  language.  The  foUowmg  presents  methods 
to  recover  from  the  most  common  problems. 

•  The  graphical  tools  make  extensive  use  of  global  variables.  Therefca^,  only  one 
graphical  tool  of  each  kind  can  be  open  at  the  same  time  in  a  Matlab  session  (i.e. 
two  vectflt  dodads  cannot  be  used  at  the  same  time). 

•  Care  must  be  taken  when  using  the  plot  command  from  the  Matlab  prompt  when 
a  dodad  is  open.  The  Matlab  plot  command  uses  the  current  axis  in  the  current  fig¬ 
ure  window  for  its  output.  If  a  dodad  is  the  current  figure  when  the  plot  command  is 
used  from  the  command  line,  the  axis  inside  the  dodad  is  used.  The  simplest  way 
to  recover  from  this  problem  is  to  close  the  dodad  and  to  start  over.  To  prevent  this 
problem,  open  a  new  figure  window  using  the  figure  command  b^ore  using  plot 
from  the  command  prompt  (in  fact,  before  any  other  command  that  creates  a  plot). 

•  In  the  event  an  oror  occurs  when  using  a  dodad,  try  closing  the  dodad  using  the 
close  pushbutton.  If  this  fails,  close  any  other  open  dodads  and  then  execute  the 
commands  “dose;  dear  globat”  to  recover. 

Matlab  Signal  Processing  Toolbox  Dependencies 

The  following  SPC  Toolbox  commands  depend  upon  various  commands  from  the 
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MatUb  Signal  Processing  Toolbox. 


Table  1:  Signal  Processing  Tbolbox  Dependencies 


SPC 

Toolbox 

Matlab  Signal  Processing  Toolbox 

vectfilt 

butter,  chebyl,  cheby2,  ellip,  hamming, 
hanning,  barlett,  blackman,  spectrum 

vectdme 

sanalyzr 

sp_steng 

sp.stmag 

sp.stza 

hamming,  hanning,  bartlett,  blackman, 
triang,  kaiser 

envelope 

modindex 

hilbert 

Installation  Under  Matlab  for  Microsoft  Windows 

The  following  steps  are  required  to  install  the  SPC  Toolbox  on  a  personal  com* 
puter  with  Matlab  running  under  Microsoft  Windows: 

1.  Create  a  directory  caUed  “SPCTCMDLS”  in  the  “SMATLABXTOOLBOX”  direc¬ 
tory. 

2.  Copy  the  SPC  Toolbox  ♦.m  files  to  the  ‘^MATLAB\TCX)LBOX^SPCTOOLS’* 
direcuxy. 

3.  Edit  the  "MATLABRC.M”  file  located  in  the  ‘^dATLAB”  directory.  Find  the 
section  similar  to  the  following  example  MAILABRC.M  file  and  insert  the  line: 

'drne_lener^mMekNgKt^^ 

The  location  of  this  line  can  be  critical  if  any  command  in  any  other  toolboxes  has 
the  same  name  as  a  SPC  Toolbox  command.  Care  has  been  taken  to  name  SPC 
Toolbox  commands  differently  from  those  used  in  the  Matlab  Signal  Processing 
Toolbox,  Control  Systems  Toolbox  m  System  Identification  Toolbox  as  of  9/93.  To 
ensure  SPC  Toolbox  commands  take  precedence  over  any  other  toolbox  com¬ 
mands  of  the  same  name,  enter  die  above  line  near  the  top  the  matlabpath  q)eci- 
fication.  To  ensure  SPC  Toolbox  commands  do  not  take  precedence,  enter  die 
above  line  near  or  at  the  bottom  of  the  maflabpalh  spedficatitm. 

%  Setup  the  MATLAB  search  pelh. 
mailabpathd... 

’e:\MATLAB;’.... 

’e:\MATLAB\iooftx>k\inBllab^oeneral;’,... 

‘e:\MATI.AB\toolx>idmallatMemo6;’.... 
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't 

'e 

D: 


\inflHafa\toQl)0)(\iionii;'.... 

\nwliab\lool)ok\id«nt;'.... 

vWBMDvOODOXVuOWOit  •*** 
vMtRHMDwOODOX\BpCiOOiBi 

\fradlalAmyprooB:'.... 


Installation  on  Sun  Workstations/Networks 
TDolbox  Directory 

1.  Create  a  ‘‘apdools’*  directory  in  the  “toolbox*’  directory. 

2.  Change  to  the  spdools  directory  using  the  UNIX  cd  command. 

3.  To  cq;}y  the  commands  from  an  MS-DOS  formatted  flof^y  disk,  use  the  mtools 
command  “moopy  -t  ‘a:*.m* .”  (period  inside  die  quotes  is  required).  The  “-t”  option 
is  required  to  convert  the  text  file  fomat  of  the  m-files  from  MS-DOS  text  file  for¬ 
mat  (o/H  end-of-line  markers)  to  UNIX  text  file  format  (H  end-of-line  marker). 

4.  Add  the  “spctools’’  directmy  to  the  maUabpalh. 

Disclaimer 

This  software  package  is  made  available  as  a  service  to  the  academic  commu¬ 
nity.  It  is  not  a  depository  of  Mathworks  approved  software.  The  authors  make  no 
warranty  of  any  kind,  eiqiressed  <x  implied,  with  regard  to  these  programs  <x  the 
documentation  contained  in  this  report.  Use  of  these  programs  is  at  the  users’  own 
risk  and  the  authors  (hsclaim  all  liability  for  injury,  damage,  and  losses  that  may 
result  from  their  use.  Losses  include,  but  are  not  limited  to,  loss  of  data  or  data 
being  rendered  inaccurate,  (x  losses  sustained  by  third  parties  for  a  failure  of  die 
software  to  operate.  The  U.S.  Government  assumes  no  responsibility  for  the  infor¬ 
mation  provided.  Matlab  is  a  trademark  of  The  Mathworks.  Microsoft  Windows  is 
a  trademark  of  Microsoft.  Sun  is  a  trademark  of  Sun  Microsystems.  Soundblaster 
is  a  trademark  of  Creative  Labs,  Inc. 
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Vector  Edit  Dodad 

This  tutorial  introduces  the  basic  operaticm  of  the  Vector  Edit  Dodad  (vectadM).  The  \fector 
Edit  Dodad  provides  an  interactive,  gri^hical  environment  widiin  Matlab  to  edit  vectors 
(naturally).  This  is  ideally  suited  for  editing  vectors  r^resenting  digitized  signals  such  as 
speech.  In  fact,  in  this  tutorial,  we  will  use  vecledH  to  leairange  the  words  in  a  speech  sig¬ 
nal. 

A  file  named  ‘^seatsiLvoc”  is  contained  on  the  distribution  disk  or  in  die  ‘‘spclools*’  direc- 
t(M7.  This  file  contains  the  words  “the  seat,  sit"  as  spoken  by  the  author.  It  was  recorded  on 
a  PC  using  a  microphone  equipped  Soundblaster  audio  board.  The  file  is  stored  in  the 
Soimdblaster  Creative  Voice  file  format.  To  load  the  file  into  Matlab,  the  ioadvoc  command 
from  the  SPC  Toolbox  is  used: 

seatsit  -  loadvocf'seaisir);  %  insert  pathname  as  appropriate 

Start  the  Vectm*  Edit  Dodad  with  the  command  “vectedit(seaisit)".  Once  loaded,  the  plot  in 
the  vectedit  dodad  window  will  look  like  the  following. 


'Tm*  (weend*) 

First,  note  that  a  60  Hz  interference  can  easily  be  seen  during  periods  of  silence.  Second, 
note  the  /$/  section  of  the  wwd  “seat"  contains  a  large,  low-fi*equency  transient.  Third, 
note  that  the  signal  was  sampled  at  8192  Hz  and  can  be  sent  to  the  audio  ouqiut  of  a  Sun 
workstation  (cn*  soundcard  equipped  PC)  by  using  the  “Play  All"  popupmenu  item. 

One  useful  function  of  vectedit  (even  when  the  user  doesn’t  want  to  actually  edit  a  vectm) 
is  to  use  vectedit  as  an  in^)ection  tool.  The  user  can  simply  load  a  vectm*  into  vectedit  and 
use  the  zoom  function  to  get  a  clearer  picture  of  a  particular  feature  contained  in  the  vec- 
tOT.  The  user  could  then  take  a  “sn4)shot”  of  the  zoomed  pcxtion  as  a  preliminary  means  of 
obtaining  a  hardcopy.  This  method  replaces  the  hit-and-miss  process  of  guessing  the  cor¬ 
rect  subscripts  to  use  with  the  Matlab  plot  command.  In  this  example,  we  will  use  the  vect¬ 
edit  dodad  to  take  a  better  look  at  the  wmrd  “seat." 
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To  zoom  in  on  the  word  ‘"seat,”  die  first  step  required  is  to  set  the  start  and  stop  markers  to 
surround  the  word.  Place  the  mouse  curstf  over  the  “Mark  Start”  pushbutton  and  fxess  the 
left  mouse  button.  At  this  time,  the  cursor  changes  into  a  cross-hair  ami  the  plot  title 
changes  to  “Mark  beginning  with  cursm.”  Move  the  cursm  until  it  is  smnewhere  i^ve  the 
small  “valley”  immediately  after  the  first  large  “burst”  as  depicted  by  the  dashed  line  in 
the  figure  below.  When  in  position,  press  the  left  mouse  button.  Once  the  mouse  button  is 
released,  the  cursor  changes  back  into  an  arrow  and  the  title  is  cleared.  A  votical  dashed 
line  appears  at  the  point  on  the  time  axis  where  the  mouse  button  was  pressed.  The  dashed 
line  is  the  beginning  marker.  Note  a  dotted  line  appears  at  the  end  of  die  signal  (not 
shown).  This  is  the  end  marker  and  since  it  has  not  b^n  placed  yet,  it  is  still  at  the  end  of 
the  signal. 

Using  a  similar  procedure,  place  the  end  marker  as  shown  below. 


I 


Now  that  the  word  “seat”  has  been  markmi,  all  menu  selections  that  apply  to  “..  Marked” 
will  apply  to  this  portion  of  the  signal.  Selecting  “Play  Marked”  plays  the  word  “seat” 
only.  To  get  a  better  look  at  the  word  “seat,”  select  “Zoom  Marked”  from  the  Zoom  pop- 
upmenu. 


tkn*  (aMonai) 


Ml 


■*8.3  0  38  0.4  0.45  0.8  0  88  0.8  O  M '  0  7  078 
Tklw^MSndi) 

The  low  frequency  transient  can  be  seen  more  clearly  now.  The  user  can  reposition  the 
markers  around  the  transient  and  select  “Zoom  Mark^”  once  again  to  get  an  ever  closer 
look.  Once  you  are  satisfied  with  the  capabilities  the  zoom  feature,  select  “Zoom  Full”  to 
view  the  whole  signal  and  to  reposition  die  markers  to  surround  the  word  “seat.” 


SPC  Toolbox  User’s  Guide 


7 


TUTORIAL 


In  re(»tlering  the  words,  we  “cut”  the  wcM-d  “seat”  and  then  “paste”  it  in  before  the  wwd 
“the.”  With  the  markers  surrounding  the  word  “seat,”  move  the  cursw  over  the  “Cut” 
pushbutton  and  press  the  left  mouse  button.  Selecting  the  cut  pushbutton  removes  the  por¬ 
tion  of  the  signal  between  the  two  markers  and  stores  it  in  a  cut-and-paste  buffo*.  It  then 
coll^>ses  the  signal  and  redraws  it  ^>propriately.  Both  markers  are  now  located  in  the 
same  spot  as  shown  below.  As  long  as  both  markers  remain  in  this  spot,  selecting  the  cut 
pushbutton  again  would  reinsert  the  cut  pcation  back  into  the  same  location.  This  “undo” 
feature  is  only  valid  as  long  as  the  markers  are  not  moved.  Note  that  selecting  the  cut 
pushbutton  has  no  effect  when  the  markers  are  located  at  the  beginning  and  end  of  the  vec¬ 
tor  (as  it  was  originally  the  case). . 


ao 

40 


.  1  -  _  _  ■  ■  ■  ■  ] 
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Tim*  (ttcoftdi) 

The  remaining  step  consists  in  pasting  the  word  “seat”  into  the  signal  before  the  word 
“the.”  Move  the  cursor  to  the  “Paste”  pushbutton  and  press  the  left  mouse  button.  After 
selection,  the  cursor  changes  into  a  cross-hair  and  the  plot  title  changes  to  “Mark  insertion 
point  with  cursor.”  Move  the  cursor  to  the  period  of  silence  before  the  wcNrd  “the”  and 
press  the  left  mouse  button.  The  word  “seat”  is  now  the  first  word  in  the  signal.  Selecting 
“Play  Full”  will  confirm  this  fact. 


'nmvfaiMond^ 

See  Common  Controls  and  VECTEDIT  far  further  discussion  on  the  use  of  the  vectedit  dodad. 
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Vector  Filter  Dodad 

This  tut(»ial  introduces  the  basic  operati(»  of  the  Vector  Hlter  Dodad  (vectfit).  Tte  Vector 
Filter  Dodad  provides  a  graphical,  interactive  environment  to  design  filters  and  to  apply 
the  filter  to  a  signal.  Note  that  vectfit  has  the  capability  to  display  the  filter  transfer  func¬ 
tion  and  the  spectrum  of  die  signal  the  filter  is  to  be  qiplied  to  simultaneously.  Further- 
mme,  the  filter  cutcff  frequencies  can  be  selected  with  ^  mouse. 

In  this  tutorial,  we  design  a  Qiebychev  l^pe  II  bandpass  filter  to  study  the  formant  fre¬ 
quencies  contained  in  a  speech  signal.  Load  the  “seatsit.voc**  file  as  discussed  in  the  vecle- 
M  tutorial  and  start  the  vectfit  dodad  with  the  command  vectfit(seal).  Once  loaded,  the  plot 
should  look  like  the  following. 


Tkn*  (weond,) 

Both  time  and  frequency  domains  are  available  in  the  vectfit  dodad.  As  could  be  expected, 
selecting  the  **llme’’  radiobutton  displays  the  time-domain  (as  shown  above)  and  selecting 
the  “Spectrum”  radiobutton  displays  the  frequency  domain  (as  shown  below).  The  spec¬ 
trum  shown  is  a  power  spectrum  estimate  across  the  entire  length  of  the  signal  using  the 
Welch  method  with  a  1024-point  window.  See  the  Matlab  Signal  Processing  Toolbox  spec- 
bum  command  for  further  information.  Only  the  positive  half  of  the  frequency  spectrum  is 
displayed,  dius,  the  frequency  resolution  is  half  the  sampling  frequency  divided  by  512. 
The  spectrum  is  normalized  as  discussed  in  the  spectrum  command  description.. 
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When  in  the  frequency  domain,  both  linear  and  logarithmic  scaling  are  available  by  select¬ 
ing  die  so  named  radiobutums.  Also,  it  is  not  necessary  to  be  in  the  frequency  domain  to 
design  filters.  All  controls  can  be  set  and  the  cutoff  frequencies  specified  using  the  edit 
boxes  while  in  the  time  domain.  However,  design  while  in  the  frequency  domain  does 
have  its  obvious  advantages. 

In  filter  design,  the  begin  marker  maiks  the  lower  cutoff  frequency  and  the  end  marks’ 
marks  the  upper  cutoff  frequency.  Place  the  lower  cutoff  frequency  in  the  valley  between 
the  first  and  second  formants  as  shown  below  (about  1200  Hz). 


Once  a  cutoff  frequency  has  been  set,  the  transfer  function  of  the  filter  specified  by  current 
settings  is  displayed.  In  the  case  shown,  a  lowpass  FIR  filter  was  specified  by  the  controls 
settings  at  the  time  the  cutoff  frequency  was  set.  Next,  place  the  end  marker  in  the  valley 
between  the  third  and  fourth  formants  (about  30(X)  Hz)  and  select  **Bandpass”  from  the  fil¬ 
ter  type  popupmenu.  This  produces  the  following  FIR  bandpass  filter. 
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Use  die  filter  popiqmienu  to  select  t  Chebychev  l^pe  I  filter  whidi  is  chartclerized  by  a 
flat  passband,  steqi  cutoff  and  an  equal  ripplt  in  the  stopband  (switching  to  a  linear  dis¬ 
play  better  shows  the  equal  rqiple). 


The  default  attenuation  in  the  stopband  is  10  dB.  Using  the  stopband  attenuation  popup- 
menu,  change  the  stopband  attenuation  to  20  dB. 


The  filter  design  is  now  complete.  Select  the  apply  pushbutton  to  filter  the  signal.  After  the 
signal  has  been  filtered,  its  spectrum  is  recomputed  and  redisplayed  along  with  the  filter 
transfer  function.  In  the  result  shown  below,  the  filter  appears  to  have  performed  quite 
well. 
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If  you  have  not  d(»»  so,  switch  over  to  a  linear  scale  for  a  quick  lode 


Now  switch  back  to  the  logaridimic  scale. 


Note  the  difference  in  the  display.  What  has  occurred  is  that  when  the  scaling  was  changed 
and  the  plot  redrawn,  the  limits  on  the  plot  axis  are  set  such  that  the  entire  magnitude  of 
the  spectrum  is  visible.  After  the  signal  was  filtered,  the  magnitude  of  the  frequencies  cor¬ 
responding  to  the  nulls  in  the  transfer  function  became  very  small  and,  as  such,  require  a 
larger  dB  axis  to  display.  You  might  wonder  why  this  axis  change  did  not  occur  immedi¬ 
ately  after  the  filter  was  implied.  The  answer  is  the  axis  was  held  to  its  pre-filtered  state  to 
allow  easier  comparison  ^  befexe  and  after  effects  of  the  filter. 

Sometimes  nulls  will  extend  so  far  down  as  to  make  it  difficult  to  see  details  of  die  qiec- 
trum  around  0  dB.  Fcxtunately,  the  magnitude  axis  can  easily  be  changed  to  effectively 
zoom  in  on  0  dB.  Unlike  zooming  in  on  the  time  or  frequency  (horizontal)  axis,  die  magni¬ 
tude  (vertical)  axis  can  be  zoomed  in  by  simply  clicking  the  mouse  widi  the  cursor  some¬ 
where  inside  the  axis.  On  a  logarithmic  scale,  clicking  die  mouse  widi  the  cursor  below 
zero  resets  the  lower  axis  limit  to  the  magnitude  the  cursor  was  at  when  the  mouse  was 
clicked.  The  same  operation  with  the  cursor  above  0  dB  sets  the  iqiper  axis  limit.  On  a  lin¬ 
ear  scale,  any  mouse  dick  with  the  cursor  inside  the  plot  axis  resets  the  u|^)er  axis  limit 
Note  that  the  curso*  must  be  in  a  blank  portion  of  the  screen  within  the  axis  when  using 
this  method  (not  near  any  curves  or  other  lines). 
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The  fdlowiQg  is  Ae  result  of  zoocniiig  in  oo  0  dB  using  die  shove  method. 


Switching  back  to  the  time  domain  illustrates  die  filter  effect.  Note  that  the  60  Hz  interftf- 
ence  and  the  low  firequency  transient  have  been  removed.  Haying  the  signal  using  the  play 
popupmenu  will  demonstrate  the  audible  contribution  the  f2  and  f3  formant  frequencies 
make  to  this  sample  of  speech. 


The  restore  pushbutton  will  reload  the  original  signal  allowing  further  experimentation. 

As  with  vectedit,  a  ‘'zoom”  feature  is  suppcated  by  setting  the  beginning  and  end  markers 
and  selecting  *^oom  Marked”  from  the  zoom  pqiupmenu.  In  vectfit,  this  feature  has  been 
extended  to  the  frequency  domain.  Markers  in  die  time  and  frequency  domains  are  inde> 
pendent  and  thus,  when  t^  display  is  zoomed  in  on  a  feature  in  the  time  domain,  the  same 
feature  is  still  zoomed  in  on  after  switching  to  the  frequency  domain  and  back. 

See  Common  Controls  and  VECTFILT  for  further  discussion  on  the  use  of  the  vectfit  dodad. 
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AR/ARMA  Modeling  Dodad 

This  tutorial  introduces  the  basic  operatitm  of  the  AR/ARMA  Modeling  Dodad  (vectarma). 
The  AR/ARMA  Modeling  Dodad  provides  a  graphical,  interactive  envirtnunent  for 

•  Autoregressive  (AR)  modeling  via  the  autocoirelatitm,  covariance,  modi- 
fied-covariance,  and  Burg  methods. 

•  Autoregressive-Moving  Average  (ARMA)  modeling  via  the  Rrony, 

Durbin,  or  Shank  methods. 

•  Performing  spectral  analysis  using  the  above  AR/ARMA  modeling  meth¬ 
ods. 

In  this  tutmial,  we  will: 

•  Create  an  AR  model  using  all  available  data. 

•  Create  an  ARMA  model. 

•  Create  an  AR  model  using  a  pmtion  of  the  available  data. 

•  Create  an  ARMA  model  of  a  speech  phoneme  and  then  reconstruct  the 
phoneme  from  the  model. 

•  PerfcMm  spectral  analysis  on  a  speech  signal. 

Contained  on  the  distribution  disk  (or  in  the  spctools  directory),  are  three  files  named 
seatsiLvoc,  uuu.mat,  and  tOx.mat  which  will  be  used  in  the  tutmial.  Each  will  have  to  be 
loaded  into  the  matlab  workspace  befcsre  using.  The  file  seatsit.voc  is  loaded  using  the 
SPC  Toolbox  command  **8681811  -  loadvoc('8ealsit');”  and  the  two  *.mat  files  are  loaded  using 
the  Matlab  commands  **load  uuu”  and  *‘load  tOx”.  The  seatsitvoc  file  contains  the  words 
**the  seat,  sit”  as  spoken  by  the  author,  the  uuu.mat  contains  the  voiced  vowel  Ai/,  and  the 
file  tOl  contains  data  sets  TOl,  T02,  T03  and  T04  from  [1]. 

Start  the  AR/ARMA  Modeling  dodad  using  the  command  **vectarma(t01,1)”.  The  second 
argument  tells  vectfit  to  start  with  the  sampling  frequency  set  to  1  vice  the  default  of  S192 
Hz.  When  the  sampling  frequency  is  set  to  1,  the  time  scale  shows  the  actual  indices  of  the 
vector  being  displayed.  The  following  plot  should  be  seen  in  vectamia. 
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Three  sets  of  marks  are  available  in  vectarma.  During  this  first  example,  we  will  leave  the 
marks  in  their  initial  position  at  the  beginning  and  end  of  the  vectm*.  From  the  AR  order 
popupmenu,  Q,  select  a  model  order  of  four.  Also,  unselect  the  spectral  analysis  plot  by 
unchecking  the  Spectrum  checkbox.  Leave  all  other  settings  to  their  default  values  and 
apply  the  model  by  selecting  the  Apply  pushbutton.  This  will  generate  a  fourth-order  AR 
model  using  the  autocorrelation  method.  After  the  model  has  been  created,  two  gri4>hic 
windows  open  containing  the  following  two  plots. 


DaagwVhteaalM  Signal 


The  upper  plot  in  the  first  graphic  window  is  the  model  (the  solid  line)  overlayed  on  the 
desired  data  (the  dotted  line).  The  lower  plot  is  the  error  obtained  by  subtracting  the  model 
from  the  desire  data. 


tad  tocaaom 


The  second  graphic  window  contains  a  polar  plot  of  the  pole  locations. 
The  Matlab  command  window  will  contain  the  model  data; 


AR  Model  -  Autocorrelation  Method 
Poles 

0.7479 -I- 0.5865i 
0.7479  -  0.5865i 
-0.3110  •f0.4334i 
-0.3110  -  0.4334i 
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Magnitucle  Angle(D) 
0.9504  38.1002 
0.9504  <38.1002 
0.5335  125.6648 
0.5335-125.6648 


Model  coefficients: 
1.0000 
-0.8738 
0.2574 
0.1362 
0.2571 


Sum-of-squared  enors  « 1 .06 

Next,  we  will  generate  an  ARMA  model  using  the  Prony  method.  From  the  Model  Type 
popupmenu  select  ARMA  and  then  set  the  MA  model  order,  Q,  to  four  and  apply.  From 
the  model/polar  plots  shown  below,  it  can  be  seen  this  method  much  more  accurately  mod¬ 
els  the  data  set  TOl. 


OsaaiU/»%ntaiaii  Sttm 


Sometimes  it  is  not  necessary  to  use  all  available  data  to  sufficiently  model  an  AR/ARMA 
process.  Return  to  AR  modeling  by  selecting  AR  fi-om  the  Model  Type  poptqjmenu. 
Located  at  the  bottom  of  the  vectarma  dodad  are  pushbutums  to  set  the  Modd  Begin  and 
Model  End  marks.  These  marks  control  data  length  used  in  generating  the  model.  Until 
now,  they  have  been  set  to  the  beginning  and  end  of  the  entire  data  set.  Using  the  Model 
&id  pushbutton,  set  the  Model  End  mark  to  the  bottom  of  the  second  cycle  of  the  data  as 
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shown  below.  Whoi  the  model  is  ^>plied  diis  time,  only  the  data  between  the  start  of  tiie 
data  (where  the  Model  Begin  mark  is  still  located)  and  tte  Model  End  maiic  will  be  used  in 
ccMnputing  the  model  parameters.  Apply  the  model  at  this  time. 


Shown  below  are  the  results  of  using  only  this  first  portion  of  data.  Compare  these  results 
to  those  obtained  in  the  first  example. 


Hsli  and  iocuoas 

•0| 


To  continue  on  with  this  tutorial,  die  data  from  the  file  uuu.dat  has  to  be  loaded  into 
vectarma.  This  can  be  accomplished  by  either  closing  the  vectarma  dodad  and  restarting  it 
with  the  command  “vectanna(uuu)”  or  by  selecting  the  Load  pushbutton  and  entering  ‘^aiu” 
into  the  edit  box  which  temporarily  rqilaces  the  Save  and  Load  pushbuttims.  To  load  uuu 
by  this  method,  uuu  has  to  exist  already  within  the  Matlab  wm-kspace  and  die  edit  box  has 
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to  be  selected  with  the  mouse  before  the  name  can  typed  into  it.  If  you  load  uuu  using  die 
Load  pushbutton,  change  the  sampling  frequency  to  8192  Hz  using  the  popiqmienu. 

The  followup  plot  should  now  be  diqilayed  in  veclanna. 


Use  the  Begin  and  End  marks  to  zoom  in  on  at  least  eight  cycles  as  shown  below.  To  pre¬ 
pare  for  this  example,  set  the  Model  Begin  mark  to  the  highest  peak  of  the  first  whole 
cycle  and  set  the  Model  End  mark  to  the  same  point  seven  full  cycles  away  from  the 
Model  Begin  mark.  Set  the  Period  End  mark  to  the  peak  of  the  cycle  immediately  after  the 
Model  Begin  mark  as  shown. 


We  are  about  to  do  the  opposite  of  what  was  done  in  the  previous  example.  In  that  exam¬ 
ple.  a  model  was  generate  that  was  longer  than  the  data  used.  In  this  example,  we  will 
generate  a  modelled  signal  that  is  shcxler  than  the  data  used  to  generate  the  model  param¬ 
eters.  The  period  markers  are  used  to  set  the  length  of  die  modelled  signal  to  be  generated 
once  the  model  coefficients  are  computed.  When  the  Model  Begin  mark  is  set,  the  Period 
Begin  maiic  is  set  to  the  same  location.  This  is  usually  how  these  marks  are  used  aldiough 
the  Period  Begin  mark  can  be  set  AFTER  setting  the  Model  Begin  mark.  The  distance 
between  the  Period  marics  is  displayed  in  the  Chain  Period  edit  box.  In  this  example,  the 
Chain  Period  rqiresents  the  pitch  period  of  the  voiced  phoneme. 

In  order  to  reconstruct  the  phoneme  from  the  model,  several  periods  of  the  model  will 
have  to  be  ‘^chained”  together.  Check  the  Chain  checkbox  and  select  a  count  of  SO  periods 
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frcMn  tiie  Chain  Length  poptq)menu.  When  the  model  is  applied,  SO  periods  of  the  model 
are  chained  together  for  form  the  rectmstructed  phoneme.  This  reconstructicm  can  be  audi¬ 
bly  compared  to  the  (xiginal  by  selecting  Play  Desired  and  Play  Model  fr(»n  the  Play  pop- 
upmenu. 

Set  vectfilt  up  to  generate  an  ARMA  model  using  the  Prony  method  with  16  poles  and  16 
zeros.  The  order  16  is  ^)ecified  by  selecting  the  User  option  from  the  Q  and  P  popup- 
menus  and  by  entering  **16”  into  the  edit  box  that  temporarily  replaces  each  popt^menu. 
The  following  shows  the  miginal  signal,  die  modelled  signal  and  die  error  sig^  obtained 
after  applying  the  model. 


The  last  example  uses  an  AR  model  to  estimate  the  frequency  content  of  seatsit.  Load 
seatsit  into  vectarma  by  one  of  the  methods  previously  suggest^  for  loading  uw.  Set  the 
controls  in  vectarma  to  generate  a  14th  order  AR  mod:  ?  the  autocorrelation  method. 

Be  sure  to  turn  the  chain  off  and  turn  the  spectrum  plot  on.  Also,  you  might  want  to  turn 
off  the  overlay  and  error  plot  options  as  they  are  of  no  use  in  this  example.  Applying  this 
model  results  in  the  following  power  spectral  density  estimate. 


In  this  plot,  the  scattered  dots  represent  the  values  of  the  fourier  transfrxm  of  the  signal. 
The  solid  line  is  the  transfer  function  specified  by  the  model  coefficient  and  is  plotted 
using  the  freqz  fimction  from  the  Matlab  Signal  Processing  Toolbox.  The  peak  frequencies 
can  be  determined  by  comparing  the  power  spectral  density  plot  to  the  polar  plot  and  the 
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Speech  Tiiiie-Doinain  Analysis  Dodad 

This  tutorial  introduces  the  basic  operatioa  of  the  Speech  ‘nme-Domain  Analytis  Dodad 
(vectlime).  The  Speech  Hme-DiMnain  Analytis  Do^  provides  a  yyhical,  imeractive 
environment  to  ^>ply  time-domain  methods  to  determine  voice  and  unvoiced  phonemes  in 
speech.  Time-domain  methods  available  in  the  SPC  Toolbox  are  the  short-time  eneify. 
short-time  magnitude  and  the  short-time  zero  crossings  mediods.  Voioed  phonemes  are 
characterized  by  a  relatively  high  energy  content  while  unvmoed  frfionemes  are  character¬ 
ized  by  higher  frequency  content  (hence,  a  higher  number  of  zero  crossings). 

In  this  tutorial,  we  determine  the  voiced  and  unvoiced  phonemes  in  the  phrase  seat, 
sir  oHitained  in  the  file  seatsit.voc.  Load  the  file  saalsit.voc  as  described  in  the  vadarma 
tutmial  and  start  the  Speech  Hme-Domain  Analysis  Dodad  with  the  command  ‘Vect- 
tnn6(sealsit)”.  Note  once  again,  the  60  Hz  interference  contained  in  this  speech  sample.  To 
accurately  apply  time-domain  analysis  to  this  signal,  the  60  Hz  interference  needs  to  be 
removed. 


T)m»(Meond4 

Without  closing  vecttimo,  start  the  vectOT  filtering  dodad  using  the  command  “vectfit 
(saatsit)”.  Design  a  10th  order,  highpass,  Chebychev  Type  I  filter  with  a  lower  cutoff  fre¬ 
quency  of  approximately  80  Hz  and  apply  the  filter.  This  will  result  in  the  removal  of  the 
60  Hz  interference  as  shown  below. 


•or 


I  *  i  I  -  i  ■  *  - — - 1 

0  0.2  a4  0.6  0.0  I  1.2  14  1 « 

Tlmafaaronda) 
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V^thout  closing  either  the  vectfit  0€  vectliine  dodads,  load  the  filtered  signal  now  ctmtained 
in  vectfit.  into  vecttime.  This  is  accomplished  by  selecting  the  Common  pushbutton  in  vect* 
time.  The  Conunon  pushbutton  is  us«l  to  load  the  latest  change  made  to  a  signal  by  either 
the  vectfit  or  vectedH  dodads  into  any  of  the  dodads  in  the  SPC  Toolbox  without  having  to 
save  the  vectw  first  and  then  starting  a  dodad  with  the  saved  vector  as  the  input  aigument. 

With  the  “clean"  signal  now  in  vecttime,  check  the  Include  Signal.  Energy,  Magnitude  and 
Zero  Crossings  checkboxes.  Also  set  the  fi-ame  length  to  IS  milliseconds  and  the  frame 
overly)  to  20%  and  then  apply.  A  long  h^tlab  gr^>hics  window  opens  with  the  following 
plots  of  the  analysis  results.  Regions  with  high  energy  or  magnitude  correspond  to  voiced 
phonemes  in  the  real-time  signal  and  the  Zero-Crossing  correspond  to  imvoiced  pho¬ 
nemes. 


Kol-teupil 


Two  smoothing  methods  are  available  to  enhance  the  curves  resulting  from  the  analysis. 
Ideally,  sharp  transitions  in  the  curves  more  clearly  define  the  begirming  and  end  of  the 
phonemes  rq)resented.  Select  a  Median  filter  from  the  Smoothing  popupmenu  and  apply. 
Notice  that  Ae  seperation  between  voiced  and  unvoiced  phonemes  is  moe  clearly  seen 
afto-  ^>plying  the  median  filter. 
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Further  manual  analysis  can  be  used  to  combine  the  output  of  either  the  short-time  eneigy 
or  magnitude  analysis  with  the  short-time  zero-crossings  analysis.  Save  the  analysis  out¬ 
put  to  a  4nable  with  the  name  "ss”  using  tte  Save  pushbutton.  This  results  in  an  Nx4 
matrix  '  .  saved  to  the  Matlab  workspace  in  which  ^  first  column  is  the  time  scale,  the 
seconc  ui  is  the  short-time  energy  arudysis  curve,  the  third  column  is  the  short-time 
magnih  >  alysis  curve  and  the  fourtt  crdtimn  is  the  short-time  zero-crossings  arudysis 
curve.  Usv  die  following  commands  to  normalize  the  short-time  energy  and  zero-crossing 
curves  and  crmibine  them  on  the  same  grifih. 


asm  a  max(S8(:,2)); 
aaz  -  max(a8(:.4)); 
t  a  S8(:,1)! 

piol(t,M('.,2ys8m,y.t,a8(:.4)/!Baz,’:’); 
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tWeCNoniMiztd  (tfiortHimt  amigy  (•)  and  zaro-CRMSingB  (:)*); 
xlabainima(a)’): 

Note  the  period  of  sileme  between  the  words  “the”  and  “seat”  is  clearly  shown. 

RecOmmendeH  «p^-h  signal  pr«v>«^ing 

[1]  John  R.  Deller,  Jr,  Jdin  G.  R-oakis.  and  John  H.  L.  Hansen,  Discrete-Time  Processing 
of  Speech  Signals,  Macmillian,  1993. 

[2]  P.  J.  Owens,  Signal  Processing  of  Speech,  hfcGraw-Hill,  1993. 

[3]  L.  R.  Rabintf  and  R.  W.  Schafer,  DigiteU  Processing  of  Speech  Signals,  Prentice-Hall, 
1978. 
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Spectrum  Analyzer  Dodad 

This  tuUxial  introduces  the  basic  operatiai  of  the  Spectrum  Analyzer  Dodad  (sanalyzr). 
The  Spectrum  Analyzer  Dodad  provides  a  graphical,  interactive  environment  designed  to 
compute  a  signal’s  spectral  cmnponents.  Its  operation  is  similar  to  that  of  die  classical 
spectrum  analyzer.  IVo  indepen^t  cursors  are  provided  alcmg  with  a  digital  readout  of 
their  corresponding  frequencies  and  magnitudes  and  the  difference  in  frequency  and  mag> 
nitude  between  diem.  Ttie  cursors  may  be  moved  from  sample-to-sample  using  the  **<“  or 
pushbuttons  or  from  peak-to-peak  (hi^  or  low)  using  the  or  *’»”  pushbutUms. 
Single  "clicking’*  on  the  spectral  line  with  the  mouse  cursor  moves  the  closest  qiectrum 
analyzer  cursor  to  the  position  of  the  mouse  cursor.  A  spectrum  analyzer  cursor  can  also 
be  grabbed  with  the  mouse  and  dragged  m  a  new  posititm.  The  spectrum  analyzer  cursors 
can  also  be  moved  to  an  specific  frequency  be  entering  the  frequency  in  the  cursor’s  fre> 
quency  edit  box  and  pressing  return. 

Note:  Movement  of  the  spectrum  analyzer  cursors  with  the  mouse  under  the  Sun  operating 
system  using  Openwindows  requires  holding  the  mouse  button  down  until  the  cursor 
movement  is  completed.  Also,  ato  dragging  a  cursor  with  the  mouse,  the  cursor  must  be 
held  down  until  the  cursor  position  stabilizes.  Under  Microsoft  ^mdows,  a  quick  click  is 
all  that  is  required  to  move  the  closest  sp^trum  analyzer  cursor  to  the  mouse  cursor  loca¬ 
tion  and  no  pause  is  required  after  dragging  a  spectrum  analyzer  cursor  to  a  new  position. 

In  this  tutorial,  we  measure  the  largest  and  smallest  magnitude  and  corresponding  frequen¬ 
cies  in  the  file  seatsit.voc.  Load  seatsitvoc  accrx’ding  to  the  instructions  given  in  the 
vectarma  tutorial.  The  following  plot  should  be  displayed  after  the  spectrum  analyzer  has 
been  started  using  the  conunand  "sanalyzrfseaisit)”. 


The  spectrum  shown  above  uses  an  enhanced  version  of  the  spectnm  command  frmn  the 
Matlab  Signal  Processing  Toolbox,  named  speclnn2.ni,  allowing  user  defined  windows 
(vice  tiie  default  hannino  window).  The  default  window  in  apadrmZ  is  still  a  hanrang  win¬ 
dow  but  the  window  can  be  changed  after  sanalyzr  is  started  by  selecting  anotiier  window 
from  die  l^dow  pulldown  menu.  The  Matlab  spectrum  command  uses  the  Welch  method 
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of  spectral  estimation  and  is  used  here  with  no  window  overiap.  The  default  PFT  length  is 
1024  points.  Since  tmly  the  positive  frequencies  up  to  half  the  sampling  frequency  are  dis¬ 
played.  In  addititm.  the  E>C  ctmiponent  is  not  displayed. 

With  die  mouse,  drag  the  spectrum  analyzer  cursma  to  within  dose  proximity  of  the  high¬ 
est  and  lowest  peaks,  as  shown  below.  When  die  cursors  are  in  the  proximity  of  the  peaks, 
the  and  pushbuttons  may  be  used  to  align  the  cursor  exacdy  on  the  peaks.  The 
frequency  and  magnitude  measurements  are  di^layed  in  the  cursor  readouts. 


Note  a  linear  scale  is  also  available  from  the  pulldown  menu  to  display  the  spectrum  under 
study. 

> 
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Finally,  note  that  the  user  can  generate  a  variety  of  baseband  and  passband  communication 
signals  using  commands  available  in  the  SPC  Toolbox  and  then  study  their  spectrums 
using  sanalyzr. 
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BPSK  Modulation/Demodulation 

This  tutorisl  focuses  on  using  commands  frmn  the  SPC  Toolbox  to  generate  a  BPSK  sig¬ 
nal  containing  a  known  message,  simulate  transmission  of  diat  signal  through  a  noisy 
channel  and  then  demodulate  and  recover  the  transmitted  message  from  die  ‘^received** 
signal.  While  this  tutorial  uses  BPSK  as  the  modulation  method.  BPSK  or  OOK  modula¬ 
tion  could  easily  be  substituted  in  its  place. 

This  tutorial  will  use  the  following  parameters: 

BPSK  modulation, 

7-bit,  even  parity  ASCI  with  one  stq>  bit  coding, 

128  bits-per-sec<md  bit  rate. 

Sampling  frequency  of  8192  Hz, 

Carrier  frequency  of  1024  Hz, 

15  dB  in-band,  signal-to-noise  ratio. 

One  popular  way  to  transmit  binary  data  is  through  the  use  of  modems  such  as  those  found 
on  most  personal  computers.  A  modem  transmits  data  one  character  at  a  time.  Each  char¬ 
acter  is  transmitted  using  one  bit  to  mark  the  start  of  each  new  character  and  one  or  two 
bits  to  mark  its  end.  The  character  itself  is  transmitted  as  either  seven  m-  eight  bits  with  an 
qitional  parity  bit  as  a  means  of  perfrxming  some  error  checking.  The  function  str2masc 
converts  a  character  string  to  a  vectix  of  I’s  or  O’s  representing  the  binary  stream  that 
would  be  sent  if  the  string  was  transmitted  by  am*  .m. 

msg  -  ‘A  boy  and  his  dog’; 
m  -  str2ma^msg,7.'eM); 
m(1;10) 
ans- 

1  (start  bit) 

1  (data  bits) 

0 

0 

0 

0 

0 

1 

1  (parity  bit) 

1  (st(^  bit) 

The  unqpolar  command  can  be  used  to  generate  a  baseband  rendition  of  this  signal. 

fs  -  8192;  Rb- 128: 

msgwave  ■  unipolar(R),fs,m): 

pio^(0:iength(m8gwave)-1  )48, msgwave); 

sat(gca.’YUm’4-OJ!5 1.25]); 

title(‘  “A  boy  and  his  dog”  -  7E1  ASCII  Coding^; 

xIabeKTime  (sec)*); 
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A>CMC<awg 


THur^m) 

IWo  commands  are  available  to  directly  create  BPSK  signals.  The  first  is  die  bpak  com¬ 
mand  but  this  only  creates  random  messages.  The  second  command,  bpslenso.  is  used  to 
create  a  BPSK  sig^  containing  a  specific  binary  message. 

fc-1024: 

p  -  hp8lansg(Rb,fc.fs.m); 

To  simulate  a  noisy  transmission  channel,  the  signal  is  mixed  with  additive  Gaussian 
white  noise  such  that  the  in-band  signal-to-noise  ratio  is  IS  decibels.  The  seisnrtMv  com¬ 
mand  performs  this  function. 

SNR.  15; 

r .  6el8nr^p,randn(lenglh(p}.1),SNR,fc,Rb,f8); 

Taking  a  look  at  how  this  signal  appears  before  and  after  ‘*transmission’’  is  easy  to  do  with 
the  plottime  command. 


BPSK  signal  witfwut  noiss 


BPSK  signtf  wWi  IS  dB  SNR 


tiins(s«c) 


siJbplot(2.1 .1  ).plotlime(p,0.05), 
tWeCBPSK  signal  without  noise'); 
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xUbeKlinw  (sacn.'yWbeK'Ampitude*); 

•ubplo((2,1 ,2).piotftne(r,0.05); 

tille(rBPSK  signal  with  ‘  in12str(SNR)  *  dB  SNRD; 

xiabeKlima  (8ec)');ylabeK‘AmpiHude’): 


Also,  plotting  the  noise-free  and  noisy  power  spectral  densities  (periodograms)  is  easy 
with  the  Iparignn  ctmimand  fw  a  log  scaled  look  (decibels)  or  wperignn  for  a  linear  scaled 
look  (relative  watts). 


subplot(2.1 ,1  ),lperianii(p) 
titleCBPSK  spectrum  without  noise’); 
subplot(2,1 ,2),lperiorm(r) 

title(rBPSK  spectrum  with '  int2str(SNR)  ’  dB  SNRl); 


Thus  far,  the  following  has  been  simulated; 


n(t) 


c(0  -  cos 


First,  the  message  was  encoded  in  binary  using  str2masc  to  produce  m(t).  Second,  m(t)  was 
keyed  as  a  baseband  BPSK  signal,  p(th  and  modulated  onto  a  carrier,  c(t),  creating  the 
transmitted  signal  -  -  *  The  bpslensg  command  actually  did  this  operation  by  converting 
m(t)  into  a  polar  srg  .u,  using  the  antpodal  command  and  then  modulating  this  signal  onto  a 
sinusoidal  carrier  using  double-sideband,  suppressed-carrier  modulation  with  the  dsbsc 
command.  Third,  Gaussian  white  noise  was  added  to  the  transmitted  signal  to  simulate  a 
noisy  transmission  channel  using  the  setsnrtw  command. 
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The  block  diagram  below  outlines  the  system,  known  as  a  product  detector,  which  is  used 
to  demodulate  the  ‘‘received”  signal. 


c(0  -  cos 


The  reader  may  recall  that  demodulation  of  PSK  signaling  can  only  be  perfonned  using 
coherent  detection.  To  the  uninitiated,  this  means  that  the  sinusoidal  signal  shown  in  the 
diagram  above  must  be  in-phase  with  Um  carrier  signal.  This  is  easy  to  do  as  the  bpekmag 
command  generates  a  carrier  with  zero  phase  shift.  As  long  as  the  sinusoid  above  is  cre¬ 
ated  with  zero  phase  shift  (at  the  same  sampling  frequency),  the  sinusoid  and  the  carrier 
are  in-phase  and  hence,  coherent  detection  can  be  performed. 

t  -  0:1/(s:(length(r)-1  )^;  c  «  oos{2  *  pi  *  fc  .*  f); 
b  -  c  .*  r; 


R«c«iv«d  BPSK  signal  mixad  with  cosina 


*^0  500  1000  1500  2000  2500  3000  3500  4000  4500 

Raquancy  (Hz) 

8ubplot(2,1 ,1),plottinie(b.0.05): 
title(’Receiv6d  BPSK  signal  mixed  with  ooskw’) 

8ubplot(2,1 ,2),lperigmi(b} 
tWeCSpeclrum  of  same') 


The  above  plot  shows  that  the  spectrum  of  b(t)  has  an  image  at  baseband  and  at  twice  the 
carrier  frequency.  Also  note  that  the  magnitude  of  the  baseband  image  is  greater  than  the 
image  at  twice  the  carrier  frequency. 

The  next  step  requires  filtering  b(t)  using  a  lowpass  filter.  This  filtering  removes  the  image 
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at  twice  the  carrier  frequency  leaving  a  baseband  polar  signal.  A  number  of  factors  ccnne 
into  play  when  selecting  the  type  of  filter  and  its  cutoff  frequency  for  use  in  diis  stq>. 
These  include  the  roll-ofif  characteristics  of  the  filter  and  the  number  signal  side  lobes  to  be 
included  in  the  passband  of  the  filter.  Including  more  signal  side  lobes  will  include  more 
harmonic  frequencies  from  the  baseband  signal  (which  is  good)  but  it  will  also  include 
more  noise  (which  is  bad).  More  haraionics  frequencies  would  make  a  “squaro*”  wave¬ 
form  but  more  noise  could  offset  that  benefit.  Obviously,  trade-offs  have  to  be  made.  Here, 
a  tenth  order,  type  1  Chebychev  lowpass  filter  with  1  dB  of  ripple  in  the  passband  is  used. 
The  passband  include  the  main  lobe  and  one  side  lobe. 

P3b,aa]  -  chebyl(10,1.4*Rb/fs); 

d-fiter(bb,aa.b): 


Original  Matsaga  Signal 


subplot(2,1 ,1  ),piottime(msgwave,0.5); 
title(‘Original  Message  Signal’); 
subplot(2,1 .2),piottime(d..5); 
titie(‘Deniodulated  Message  Signal’); 

As  can  be  seen,  the  result  indeed  resembles  a  noisy  polar  waveform.  Comparing  this  plot 
to  the  first  plot  of  the  binary  waveform  reveals  a  phase  delay.  This  delay  is  caused  by  the 
non-causal  property  of  the  lowpass  IIR  filter  used  in  this  tutorial. 

The  final  step  is  to  recover  the  information.  To  do  this,  the  signal  d(t)  is  resampled  taking 
one  sample  per  bit.  To  properly  resample  d(t),  care  must  be  taken  to  sample  each  bit  only 
when  it  has  stabilized  after  a  bit  transition  but  before  a  transition  to  the  next  bit  starts. 
Because  of  the  filter  delay,  d(t),  can  simply  be  sample  at  times  when  bit  transitions  would 
have  occurred  (every  f^Rb  samples). 

k  ■  d(fsyRb:fs/Rb:length(d));  %  fe/Rb  -  samples-per-bit 

Once  resampled,  k(t)  still  represents  an  analog  “voltage.”  An  analog  system  would  use  a 
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threshold  detector  at  this  point  to  convm  the  analog  signal  to  a  binary  signal.  A  logical 
comparison  operator  in  I^tlab  can  be  used  to  simulate  a  threshold  detector.  In  diis  case, 
the  optimal  threshold  voltage  is  zero  and  everytime  the  signal  k(t)  is  greater  than  this 
value,  it  is  ctmsidered  a  binary  1. 

k-k>0;  %  convert  to  binary 

This  results  in  a  binary  representation  of  the  received  message.  This  biruuy  repiesentaticm 
can  be  compared  to  the  binary  representation  of  the  miginal  message  to  compute  the  num¬ 
ber  of  bit  errtvs. 

bitenor  ■  8um(k  m) 
biterror - 
0 

Last,  the  masc2str  command  can  be  used  to  convert  the  received  biiuuy  signal  into  a  read¬ 
able  fcHin. 


masc28tr(k,7,'e*.1) 

ans- 

A  boy  and  his  dog 


For  further  study: 

1.  Change  various  parameters  in  the  above  simulation  and  observe  the  results.  Candidate 
parameters  include; 

Signal-to-noise  ratio. 

Modulation  method. 

Lowpass  filter  type,  cutoff  frequency,  attenuation,  ripple,  etc. 

2.  Add  interfering  signals  (i.e.  voice,  BPSK,  etc.)  in  addition  to  the  Gaussian  white  noise 
and  attempt  to  recover  the  modulated  information.  This  enviromnent  can  be  generated  by 
first  creating  an  environment  containing  the  primary  signal  with  Gaussian  white  noise  (as 
in  this  tutorial),  and  then  generating  a  second  environment  by  creating  the  interfering  sig¬ 
nal  and  using  the  first  environment  as  the  noise  input  to  the  seisnrbw  command. 

3.  Pad  the  begiiming  and  end  of  the  signal  with  arbitrary  amounts  of  white  noise  (scaled 
appropriately  with  the  o  scaling  factor  returned  by  tiie  seisnrbw  command).  Write  a  com¬ 
mand  to  “synch-up”  with  the  signal  in  order  to  use  coherent  demodulation. 

4.  Experiment  using  error  correcting  codes  and  low  SNR’s. 

5.  Instead  of  subsampling  the  demodulated  signal,  try  using  an  integrator  w  some  other  bit 
detection  scheme. 
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ANTPODAL 


Purpose 

Generate  a  baseband,  antipodal  signal. 

Synopsis 

y  -  antpodal(ft),moo): 
y  -  antpoclaKfb,mao.fB); 

Description 

anlpodai(ft),msg)  •  Generates  a  baseband  antipodal  signal  at  bit  rate,/!).  Sampling 
frequency  defaults  to  8192  Hz.  msg  is  a  scalar,  a  message  is  generated  as  a  ran¬ 
dom  sequence  of  length  msg  widi  Pr(-l)sPr(l>O.S.  If  msg  is  a  vecto*,  it  must  be  a 
v^'^tor  olf  O’s  and  I’s.  The  length  of  the  ouq>ut  vector,  y,  can  be  computed  by  the 
f  .mula: 

length  m  nbrbaud  - floor 

antpodaKlb,msg.fs)  -  Sampling  frequency  is  set  to  fs. 

Note:  If  the  bit  ra.te,fb,  and  the  sampling  frequency,  fs,  are  both  equal  to  one,  the 
resulting  vector  will  contain  alternating  values  of  [-1,1]. 

Example 

Generate  30  bits  of  a  SO  bit-per-second  signal  sampled  at  8000  Hz.  Use  this  base¬ 
band  signal  as  the  modulating  signal  input  to  dsbac  to  generate  phase-shift  keyed 
signal. 


X  -  antpodaK50.30,8000); 
y-dsb6c(x.2000.8000); 


Limitations 

•  The  probability  1  occurs  will  not  be  exact  for  a  small  number  of  bits. 

•  If  the  sampling  frequency  is  not  an  exact  multiple  of  the  bit  rate,  unexpected 
results  may  occur.  Fw  example:  antpodal(9, 18,20)  will  produce  an  output  of  36 
samples  vice  the  expected  40  samples  (i.e.  18  /  9  =  2  seconds  at  20  samples-per- 
sec<^).  This  occurs  since  floor(fslfb}  =  2  samples-per-bit. 

See  Also 

unipolar,  aqwave 
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ARJURG _ 

Purpose 

Compute  pirameters  for  an  AR  model  using  the  Biug  method. 

Synopsis 

[a.bO.S.oammajioma]  •  ar_bino(x.P) 

Description 

ar_biJig(x,P)  -  Computes  the  coefficients  of  an  order  P  Auto-Regressive  Model 
(AR)  of  the  data  given  in  vector  x  via  the  Burg  method  as  put  forth  in  [1],  The  out¬ 
put  is  specified  by; 

a  «  filter  coefficients 
bo^gun 

S  s  prediction  etrw  variance 

The  gain  term,  %  is  not  computed  in  the  Buig  algorithm  itself.  It  is  computed 
using  ''brute  force”  by  generating  the  impulse  response  of  the  model  and  fiien 
Inlying: 

gain  * 

V  model) 

Example 

Estimate  the  parameters  for  a  second  rxder  AR  model  for  the  data  sequence  { 1,  -2» 
3,  -4, 5}.  (Exjunple  9.6  in  [1].) 

[a.bO,^  •  arjbuigdl  -23-4  5|^); 
a-  1.0000 

2.0000 
1.0000 
bO-  1.0000 

S-  2.84226-014 

Algorithm 

L  “  J  iff-i  J 

See  Also 

ar.conr.  ar.oovar,  ar.mdoov,  ar_8vd,  ar_prony.  arjduibn,  ar_shank 

Reference 

[1]  Charles  W.  Thenien,  Discrete  Random  Signals  and  Statistical  Signal  Process¬ 
ings  PP  545-548,  Prentice-Hall,  1992. 
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AR^CORR _ 

Purpose 

Compute  parameters  for  an  AR  model  using  the  autocorrelatimi  method. 

Synopsis 

Ia.bO.S.4  -  ar_ooff(x.P) 

Description 

arjoorr(x,P)  -  Computes  the  coefScients  of  an  order  P  Auto-Regressive  Model  (AR) 
of  the  data  given  in  vector  x  via  the  autocmrelation  method  as  put  forth  in  [1].  The 
ou^ut  is  specified  by: 

a  -  filter  coefficients 
bOsgain 

5  s  minimum  sum  of  squares 
s  =  estimate  fcH*  prediction  error  variance 

Example 

Estimate  the  parameters  for  a  second  order  AR  model  fm*  the  data  sequence  { 1,  -2, 
3,  -4, 5}.  (Example  9.6  in  [1].) 

[a,S.sl  -  ar_oonr(I1  -23-4  5l;2): 

a- 

1.0000 

0.8140 

0.1193 

S  ■ 

25.5404 

8« 

3.6486 

Algorithm 

»„  =  JS  x“xc  =  I 

See  Also 

ttrjoovar,  ar_mdoov,  ar_burg,  arjsvd,  ar_prony,  arjdurbn,  ar.shank 

Reference 

[1]  Charles  W.  Theirien,  Discrete  Random  Signals  and  Statistical  Signal  Process- 
ing,  pp  535-541,  Prentice-Hall,  1992. 
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AR,COVAR _ 

Purpose 

Compute  parameters  for  an  AR  model  using  the  covariance  mediod. 

Synopsis 

[a.bO.S.<  -  ■r_oovar(x,P) 

Description 

arjoovar(x.P)  •  Computes  the  coefficients  of  an  order  P  Auto-Regressive  Moctel 
(AR)  of  the  data  given  in  vector  jr  via  the  covariance  method  as  put  forth  in  [1]. 
The  ou^ut  is  specified  by: 

a  -  filter  coefficients 
bo- gun 

S  -  minimum  sum  of  squares 
5  -  estimate  for  prediction  error  variance 

The  gain  term,  is  not  computed  in  the  covariance  algorithm  itself.  It  is  com¬ 
puted  using  brute  f^ce  by  generating  the  impulse  response  of  the  model  and  then 
Implying: 

gain  *  sqrt(p^] 

VmodelJ 

Example 

Estimate  the  parameters  for  a  second  order  AR  model  for  the  data  sequence  { 1,  -2, 
3,  -4, 5}.  (Example  9.6  in  [1].) 

[a,bO,S.^  -  ar.oovarQI  -23-4  51.2); 
a-  1.0000 

2.0000 
1.0000 
bO-  1.0000 

S-  2.8422»013 

S-  5.6843«014 

Algorithm 

\pj 

See  Also 

arjoorr,  ar.mdoov.  arjxirg,  ar.avd,  arjxony,  arjAabn.  arjshank 

Reference 

[1]  Charles  W.  Therrien,  Discrete  Random  Signals  and  Statistical  Signal  Process¬ 
ing,  pp  535-541,  Prentice-Hall,  1992. 
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AR_DURBN 


Purpose 

Compute  parameters  for  an  ARMA  model  using  the  Durbin  method. 

Synopsis 

la^Ji  -  ar_duitn(x.P.Q.L) 

Description 

ar_durt)n(x,P,Q.L)  •  Computes  the  coefficients  of  an  order  P  Auto-Regressive,  order 
Q  Moving  Average  Model  of  the  data  given  in  vector  x.  The  L  parameter  can  be 
used  to  set  the  order  of  the  intermediate  AR  model  used  in  the  algorithm.  If  ncrt 
provided.  L  defaults  to  5  *  Q  of  tltt  length  of  x  minus  I,  whichever  is  smaller.  The 
intermediate  AR  model  is  first  found  using  the  covariance  mediod.  This  model  is 
them  used  to  generate  new  data  frmn  which  the  MA  parameters  are  found  using 
the  Durbin  method.  The  AR  parameters  are  found  using  the  covariance  method. 
The  ou^ut  arguments  are: 

a  s  AR  parameters 
b  =  MA  parameters 


Limitations 

The  Durbin  method  is  not  well  sui^  for  small  data  sets. 

See  Also 

arjoorr.  arjoovar,  ar_mdoov.  ar.buig.  ar_svd,  ar_prony,  ar_shank 

Reference 

[1]  Charles  W.  Theirien,  Discrete  Random  Signals  and  Statistical  Signal  Process¬ 
ing,  pp  558-560,  Plentice-Hall,  1992. 
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AR.LEVIN 


Purpose 

Compute  parameters  for  an  AR  model  using  the  Levinson  recursion. 

Synopsis 

lMl-ar_l«vin(r.P) 

Description 

arJlovin(r.P)  •  Computes  die  coefficients  of  an  order  P  Auto-Regressive  Model  from 
values  cf  the  autocorrelation  vector  r  =  fR(0)  R(I)  R(2) ...  R(P-1)J.  The  autocorre¬ 
lation  must  be  symmetric  {R(l)  =  R*(-l),  R(2)  =  R*(-2),  etc,.).  The  ou^ut  argu¬ 
ments  are: 


a  -  AR  parameters 
s  =  linear  prediction  error  variance 

Example 

Given  rs  [3  21000 ...],  find  the  second  order  AR  model.  (Example  8.1  in  [1]) 

[a.s]  -  ar_levin(P  2 1  Q,2); 
a-  1.0000 

-0.8000 
0.2000 
S-  1.6000 

Algorithm 

Initialization: 


Recursion: 


See  Also 


1-0  -^,[11 


-1  o^-a,[0j 


a  *7  ^ 


0 


-Y, 


0 


a 


^-1 J 


o  •  o*  ,  -T  A* 

p  p~i  'p  p 


arjocMT.  arjoovar,  ar.mdoov,  ar.burg,  ar_prony,  arjdurbn.  ar.shank 


Reference 

[1]  Charles  W.  Therrien,  Discrete  Random  Signals  and  Statistical  Signal  Process- 
ing,  pp  422-430,  Prentice-Hall,  1992. 
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AR.MDCOV 


Purpose 

C(xnpute  parameters  for  an  AR  model  using  the  modified  covarianoe  method. 

Synopsis 

[a.bO.S.i4  -  ar_mdoov(x.F0 

Description 

ar.mdoov(x,P)  •  Computes  the  coefficient  of  an  order  P  Auto>Regressive  Kfodel 
(AR)  of  the  data  given  in  vector  x  via  the  modified  covariance  method  as  put  forth 
in  [1].  The  ou^ut  is  specified  by: 


a  s  filter  coefficients 
bo  s  gain 

5  =  minimum  sum  of  squares 
s  =  estimate  for  prediction  errcM’  variance 


The  gain  term,  is  not  computed  in  the  modified  covariance  algmithm  itself.  It  is 
computed  using  brute  force  by  generating  the  impulse  response  of  the  model  and 
then  applying: 


gain  =  sqrt 


(A] 

V  model) 


Example 

Estimate  the  parameters  for  a  second  order  AR  model  for  the  data  sequence  { 1.  -2, 
3,  -4, 5}.  (Example  9.6  in  [1].) 

[abO,S.s]  -  ar.mdoovfil  -2  3-4  51,2); 
a-  1.0000 

2.0000 
1.0000 
bO-  1.0000 

S-  2.84226-014 

8  -  5.68436-015 

Algorithm 


)o 


p. 


See  Also 


ar.corr,  ar.oovar,  ar.buig.  ar.svd,  arjirony,  arjdurbn,  ar.shank 

Reference 


[1]  Charles  W.  Thenien.  Discrete  Raruhm  Signals  and  Statistical  Signal  Process- 
ing,  pp  542-544,  Prentice-Hall,  1992. 
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AR_PRGNY _ 

Purpose 

Compute  parameters  for  an  ARMA  model  using  the  Prony  method. 

Synopsis 

|a.tH  -  ar_prony(x.P.Q) 

Description 

ar_prony(x,P,Q)  -  Computes  the  coefficients  of  an  order  P  Auto-Regressive,  order  Q 
Moving  Average  Model  of  the  data  ^ven  in  vector  x.  The  AR  parameters  are 
found  using  the  covariance  method.  MA  parameters  are  found  using  the  Prony 
method.  The  output  arguments  are: 

a- AR  parameters 
b  =  MA  parameters 


See  Also 

ar_ooa,  ar_oovar,  ar_mdoov,  ar.burg.  ar_svd,  arjdurtxi,  ar_shank 

Reference 

[1]  Charles  W.  Thenien,  Discrete  Random  Signals  and  Statistical  Signal  Process- 
ing,  pp  550-555,  Prenticc-Hall,  1992. 
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AR^SHANK 


Purpose 

Compute  paramettt^s  for  an  ARMA  model  using  the  Shank  method. 

Synopsis 

[a.t4-«_shank(x.P,Q) 

Description 

arjBhank(x,P.Q)  -  Computes  the  coefficients  of  an  order  P  Auto-Regressive,  order  Q 
Moving  Average  Model  of  the  data  given  in  vector  jc.  The  AR  parameters  are 
found  using  the  covariance  method.  MA  parameters  are  found  using  the  Shank 
method.  The  output  aiguments  are: 

a  =  AR  parameters 
b  =  MA  parameters 


See  Also 

ar.oorr,  ar_oovar,  ar.mdoov,  ar_buig,  ar_svd.  arjxony,  arjdurbn 

Reference 

[1]  Charles  W.  Thcrrien,  Discrete  Random  Signals  and  Statistical  Signal  Process- 
ing,  pp  555-558,  Prentice-HaU,  1992. 
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AVSMOOTH,MDSMOOTH _ 

Purpose 

Snux^  a  curve  using  an  average  or  median  smoothing 

Synopsis 

y  -  avstnoolh(x.L); 
y  -  md8iiioolh(x.L): 

Description 

avsnKXMh(x,L)  -  Smooths  the  curve  specified  in  vector  x  using  an  average  smoothing 
filter  of  length!. 

md8moolh(x.L)  -  Smooths  the  curve  specified  in  vector  x  using  an  median  smooth¬ 
ing  filter  of  lengdi  L. 

Example 

Compute  the  short-time  energy  of  the  speech  signal  specified  in  x  and  enhance  the 
resulting  ciuve  using  an  average  smoother  of  length  7. 

y  -  8p_8teno(x.64,30, 'hamming'); 
yy  -  av8mooth(y.7)^ 
yyy  *  mdsmoo^.T); 

Unsmoothed  cuve 
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Algorithm 

A  non-causal  smoothing  algorithm  is  used.  If  L  is  odd,  samples  before  and 
after  the  current  sample  are  used  in  computing  the  average.  If  L  is  even,  L/2  sam¬ 
ples  before  and  ^  - 1  samples  after  the  current  sample  are  used  in  computing  the 
average.  It  is  recommended  an  odd  number  for  L  be  used  to  prevent  any  biasing  an 
even  number  for  L  can  cause. 


I  L-6  I 


I  L»7  I 


Limitations 

Zero  padding  is  used  on  either  end  of  the  input  vector  and  then  the  vector  is 
trimmed  to  produce  a  non-causal  filter.  Therefore,  samples  within  L/2  samples 
from  either  end  will  be  subject  to  a  transient  which  distorts  the  curve  toward  zero. 
This  distortion  becomes  more  exaggerated  the  further  the  absolute  value  of  the 
samples  in  this  region  are  from  zero. 

See  Also 

mdsmooth 
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BFSK,  BFSKMSG 


Purpose 

Generate  binary  frequency-shift  keyed,  bandpass  signals. 

Synopsis 

y  ■  bfek(Rb,shifl.fc); 
y  -  bfsk(Rb.8hifl,fc,duraiion); 
y  -  bf8k(Rb.8hift.fc.duration,fB); 
y  -  b(8km8g(Rb,8hift.fc.msg); 
y  •  bfskmsg(Rb,shift,fc,fS,m8g): 

Description 

b(sk(Rb,8hift,fc)  -  Generates  one  second  of  a  binary  frequency-shift  keyed  signal 
with  a  bit  rate  of  Rb  bits-per-second  and  a  frequency  shift  of  shift  centered  at  a  car¬ 
rier  frequency. /c.  Default  sampEng  rate  is  8192  Hz. 

bfsl<(Rb,8hift,fc, duration)  -  Generates  duration  seconds  of  the  binary  frequency-shift 
keyed  signal. 

bfsk(Rb,shift,fc,duration,f8)  -  Set  the  sampling  frequency  to  fs. 

bfsKmsg(Rb,shift.fc,m&g)  or  bfskmsg(Rb,shift.fc,fs,msg)  -  Generate  a  bandpass  binary 
frequency-shift  keyed  signal  with  a  bit  rate  of  Rb  bits-per-second,  a  frequency 
shift,  shift,  centered  at  a  carrier  frequency, /c,  sampled  at  sampling  frequency, /a:, 
and  containing  the  message,  msg.  U  fs  is  not  given, /5  defaults  to  8192  Hz. 

Example 

Generate  a  64  bit-per-second,  binary  frequency-shift  keyed  signal  at  a  carrier  fre¬ 
quency  of  2048  Hz  and  a  frequency  shift  of  128  Hz  lasting  one-and-one-half  sec¬ 
onds.  Play  the  signal  over  the  workstation’s  speaker. 

y  -  bfsk(64.128.2048.1.5); 
sound(y,8192); 

Algorithm 

The  binary  frequency-shift  keyed  signal  is  generated  by  first  creating  two  on-off 
keyed  signals  at  the  desired  biVsampling  rates  and  combining  them.  The  on-off 
keyed  signals  are  generated  such  that  the  mark  signal  is  generated  tAfc+shiftl2  and 
the  space  signal  is  generated  at  fc-shiftll.  After  the  signals  are  generated,  they  are 
added  to  form  a  BFSK  signal.  Note:  if  shift  is  an  integer  multiple  of  the  bit  rate 
{shift^n*Rb  where  n=!l,2,...),  the  resulting  BFSK  signal  is  coherent  (i.e.  has  con¬ 
tinuous  phase). 

Limitations 

See  limitations  imder  imipolar. 
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See  Also 

unipolar,  dsbac,  ook,  ooloTiag.  bpsk.  bpekiwg 

Reference 

[1]  Le<m  W.  Couch,  Digital  and  Analog  Communication  Systems,  pp.  332-344. 
351-356,  Macmillan.  1990. 


SPC  Toolbox  User’s  Guide 


45 


BPSK.BPSKMSG 


BPSK,  BPSKMSG 


Purpose 

Generate  binary  phase-shift  keyed,  bandpass  signals. 

Synopsis 

y  -  bpsk(Rb.fc); 
y  -  bp6k(re>,fc,duralion): 
y  •  bp6k(Rb,fc,duraten,fs); 
y  -  bp6k(Rb.fc;msg); 

y  -  bp6km8g(Rb,fc.fs.m8o): 

Description 

bp6k(Rb,fc)  -  Generates  one  second  of  a  binary  phase-shift  keyed  signal  with  a  bit 
rate  of  Rb  bits-per-second  centered  at  a  carrier  frequency,^.  Default  sampling  rate 
is  8192  Hz. 

bp6k(Rb,fc.duration)  -  Generates  duration  seconds  of  a  phase-shift  keyed  signal. 

bp6k(Rb,fc,duration.fs)  -  Sets  the  sampling  frequency  to/5. 

bp6kmsg(Rb,fc.msg)  or  bp6kmsg(Rb,fc,fs,msg)  -  Generates  a  binary  phase-shift  keyed 
signal  with  a  bit  rate  of  Rb  bits-per-second,  at  a  carrier  frequency, /c,  sampled  at 
sampling  frequency,  fs,  and  containing  the  message,  msg.  If  fs  is  not  given,  fs 
defaults  to  8192  Hz. 

Example 

Generate  a  64  bit-per-second,  binary  phase-shift  keyed  signal  at  a  carrier  fre¬ 
quency  of  2048  Hz  lasting  one-and-one-half  seconds.  Play  the  signal  over  the 
workstation’s  speaker. 

y-bpsk(64.2048,1.5): 

80und(y,8192): 

Algorithm 

The  binary  phase-shift  keyed  signal  is  generated  by  first  creating  a  binary  antqK)- 
dal  signal  at  the  desired  ^sampling  rates  and  then  applying  a  double-sideband, 
suppressed-cairier  modulation  scheme  to  the  resulting  baseband  signal  at  the 
desired  carrier  frequency. 

Limitations 

See  limitations  under  antpodai. 

See  Also 

bpakmag,  ook,  ootansg,  bfsk,  bfskmag 

Reference 

[1]  Leon  W.  Couch,  Digital  and  dialog  Communication  Systems,  pp.  332-336, 
Macmillan,  1990. 
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This  toolbox  contains  a  number  of  interactive  “dodads”  driven  by  pc^upmenus, 
pushbuttons  and  other  graphical  controls.  Several  of  these  controls  are  common  to 
a  number  of  dodads  and  are  discussed  below.  For  specific  instructions  on  the  use  of 
Matlab  pc^upmenus,  pushbuttons,  edit  boxes,  and  sliders,  refer  to  discussions  of 
“Hande  Graphics”  in  the  Matlab  User’s  Manual  and  the  “uiooritrol”  statement  in  the 
Matlab  Reference  Manual  for  your  system. 

Mark  Start/Mark  End  Pushbuttons 

Purpose 

•  Marks  the  beginning  and  end  of  a  vector  section  to  be  cut,  copied  or  otherwise 
acted  upon. 

•  Marks  cutoff  frequencies  for  filters  (lower  cutoff  =  Mark  Start,  upper  cutoff  = 
Mark  End). 

Use 

1.  Select  desired  mark  pushbutton  (Start  or  End). 

2.  After  selection,  the  cursor  changes  to  a  crosshair  and  the  plot  title  changes  to 
either  “Mark  start  with  cursor”  or  “Mark  end  with  cursor"  as  ^>propriate. 

3.  Move  the  cursor  to  the  desired  mark  location  and  click  the  left  mouse  button. 
When  placing  the  mark,  only  the  location  on  the  horizontal  (time  or  frequency) 
axis  is  of  importance.  Placing  a  mark  past  either  end  of  the  plot  axis  sets  the  mark 
to  the  first  or  last  value  displayed  on  the  horizontal  axis. 

Result 

After  placing  a  mark,  the  curscx*  changes  back  into  an  arrow  and  the  plot  title  is 
cleared.  The  start  mark  is  displayed  as  a  vertical  dashed  line  and  the  end  mark  is 
displayed  as  a  vertical  dash-dot  line. 

Reset  Marks  Pushbutton 

Purpose 

Reset  marks  to  the  extremes  of  the  current  display.  On  a  “Zoom  Full”  display. 
Reset  sets  the  marks  to  the  begiiming  and  end  of  the  vector.  On  a  “Zoom  Marked” 
display.  Reset  sets  the  marks  to  the  beginning  and  end  of  the  zoomed  section  of  the 
vector. 

Use 

1.  Click  the  left  mouse  button  while  the  cursor  is  over  the  Reset  Marks  control. 

Result 

The  respective  marks  are  changed  to  the  first  or  last  value  displayed  on  the  axis. 
When  marks  are  at  the  axis  extremes,  they  may  or  may  not  be  visible  depending 
up>on  whether  or  not  the  plot  border  obscures  them.  (This  is  a  problem  in  version 
4.0a  of  Matlab  for  Sun  workstations  that  has  been  fixed  in  version  4.1.) 

Play  Popupmenu 

Purpose 

Send  either  the  entire  vector  or  marked  section  thereof  to  the  workstation’s  audio 
output. 
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Use 

1.  If  desired,  mark  a  portion  of  the  vector  using  Mark  Start  and  Mark  End. 

2.  Place  the  cursor  on  the  Play  popupmenu  and  click  the  left  mouse  buttcm  once  to 
<^)en  the  menu. 

3.  Place  the  cursor  on  the  desired  function  (Play  Full  or  Play  Marked)  and  click  tiie 
left  mouse  button. 

Result 

Selecting  Play  Pull  outi>uts  the  entire  vector  to  the  audio  outi)ut.  Selecting  Play 
Maiked  ouq)uts  the  portion  of  the  vector  between  die  marks  to  the  audio  output, 
the  marks  have  not  been  set,  the  entire  vector  is  played. 

Notes 

•  On  a  Sun  wwkstation,  the  ou^ut  volume  and  device  (speaker  or  hea(4>h(»ie  jack) 
are  set  using  the  Sun  operating  system’s  AudioTool.  The  output  sampling  fre¬ 
quency  is  fixed  at  8,192  Hz. 

•  On  a  PC  under  Microsoft  \^dows,  the  ouq}ut  device  is  to  the  soundcard,  if  one 
is  installed  and  supported  by  Matlab.  Tools  which  accompany  the  soundcard  are 
used  to  set  the  output  volume.  The  output  sampling  frequency  is  set  by  the  sam¬ 
pling  frequency  popupmenu. 

Zoom  Popupmenu 

Purpose 

•  Zooms  the  display  in  on  a  marked  section  of  the  vector. 

•  Zooms  the  display  out  to  show  the  full  vector. 

Use 

1.  If  desired,  mark  a  portion  of  the  vector  using  Mark  Start  and  Mark  End. 

2.  Place  the  cursor  on  the  Zoom  popupmenu  and  click  the  left  mouse  button  once 
to  open  the  menu. 

3.  Place  the  cursor  on  the  desired  function  (Zoom  Full  or  Zoom  Marked)  and  click 
the  left  mouse  button. 

Result 

The  requested  acti  is  performed.  Zoom  Marked  can  be  repeatedly  used. 

Restore  Pushbutton 

Purpose 

Restores  the  (xiginal  vector  the  dodad  was  called  with  if  a  vector  argument  was 
supplied.  If  the  dodad  was  called  without  an  aigiiment,  the  common  vectOT  at  the 
time  the  dodad  was  called  is  restored  (see  Common  Pushbutton). 

Use 

1.  Place  the  cursor  on  the  Restixe  pushbutton  and  click  the  left  mouse  button. 

Result 

The.  cdginal  vectm*  is  restored.  All  prior  edits  are  lost  unless  previously  saved  to 
the  Matlab  wm-kspace  (see  Save  pushbutton). 

Common  Pushbutton 

Purpose 

A  "common”  vectOT  is  maintained  between  dodads  to  facilitate  loading  the  current 
vectm*  from  one  dodad  into  others  without  having  to  go  through  a  save,  quit,  and 
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restart  sequence.  This  ctunmon  vector  is  updated  after  any  operatkm  that  modifies 
a  vector.  For  example,  using  cut  and  paste  in  the  VKtor  Edit  Dodad  or  filtering  a 
vector  using  the  Vector  Filter  Dodad  are  examples  of  opeimtitms  that  change  a  vec> 
tor.  Suppose  a  vecux*  is  filtered  using  the  vectw  filter  dodad  and  the  vector  edit 
dodad  is  already  q>en.  After  the  filtering  is  complete,  selecting  the  Commcm  push¬ 
button  loads  the  filtered  signal  into  the  vector  edit  dodad. 

Use 

1 .  PerfcMm  an  q)eration  that  changes  a  vector  in  one  of  the  dodads. 

2.  Place  the  cursm*  on  the  Common  pushbutton  and  click  the  left  mouse  buttcm. 

Result 

A  c<^y  of  the  current  common  vector  is  loaded  into  the  dodad. 

Save  Pushbutton 
Purpose 

Saves  the  current  vector  to  the  Matlab  wOTkspace  under  a  new  variable  name.  This 
is  the  only  way  to  save  any  changes  made  to  a  vector. 

Use 

1.  Place  the  cursor  on  the  Save  Pushbutton  and  click  the  left  mouse  button. 

2.  The  Save  and  Load  pushbuttons  are  replaced  by  an  edit  box.  If  save  has  already 
been  used,  the  previously  entered  vector  name  are  contained  in  the  edit  box. 

3.  Place  the  cursor  on  the  edit  box  and  click  the  left  mouse  button.  Using  the  key¬ 
board,  backspace  over  any  previous  entry  and  then  type  in  the  desired  vector  name 
(following  standard  Matlab  variable  name  conventions). 

4.  Press  the  return  key. 

Result 

The  vector  is  saved  to  the  Matlab  workspace. 

Load  Pushbutton 
Purpose 

Loads  a  vector  from  the  Matlab  workspace  into  the  dodad. 

Use 

1.  Place  the  cursor  on  the  Load  Pushbutton  and  click  the  left  mouse  button. 

2.  The  Save  and  Load  pushbuttons  aie  replaced  by  an  edit  box.  If  load  has  already 
been  used,  the  previously  entered  vector  nameis  contained  in  the  edit  box. 

3.  Place  the  cursor  on  the  edit  box  and  click  the  left  mouse  button.  Using  the  key¬ 
board,  backspace  over  any  previous  entry  and  then  type  in  the  desired  vectcx  name 
(following  standard  Matlab  variable  name  conventions). 

4.  Press  the  return  key. 

Result 

The  vector  is  loaded  from  the  Matlab  workspace. 

Notes 

The  vector  has  to  exist  in  the  Matlab  wwlcspace.  If  the  vector  is  not  in  the  Matlab 
workspace  or  the  variable  name  is  misspelled,  the  error  message  "Invalid  variable 
name  or  vector  not  in  workspace!”  is  printed  in  the  Matlab  command  window  and 
the  Save  and  Load  pushbuttons  are  restored.  The  Load  pushbutton  needs  to  be 
reselected  to  retry. 
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Close  Pushbutton 
Purpose 

Qoses  the  dodsd  and  any  child  graphic  windows  that  may  be  open. 

Use 

1 .  Place  the  cursor  on  the  Close  Pushbutton  and  dick  the  left  mouse  butum. 

Result 

Any  child  gr^>hic  windows  are  closed  along  with  the  dodad  itself.  The  vector  is 
NOT  saved  upon  exit.  To  save  the  vector  before  exiting,  see  Save  pushbutUm. 

Snapshot  Pushbutton 

Purpose 

Opens  a  new  figure  window  and  recreates  the  current  dodad  plot,  effectivdy  taking 
a  ‘‘snapshot.”  This  new  figure  window  is  then  forgotten  by  the  dodad  and  remains 
(^n  after  the  dodad  is  dosed.  Once  the  new  figure  window  has  been  created,  die 
figure  window  can  be  treated  like  any  other  Matlab  figure  window  and  either 
dosed,  edited  or  printed.  Hiis  allows  the  user  to  take  a  snapshot  of  the  current 
dodad  plot,  add/change  its  tide  and  then  print  the  image  to  a  postscript  file  or 
printer,  as  desired. 

Use 

1 .  Place  the  cursor  on  the  Snapshot  pushbutton  and  click  the  left  mouse  button. 

Result 

The  current  plot  is  recreated  in  a  new  figure  window  as  discussed  above.  Before 
printing  or  otherwise  handling  this  new  figure,  it  should  be  made  the  current  figure 
with  the  use  of  the  command  ligureih)  here  h  is  the  figure  window  number.  It  is  up 
to  the  user  to  close  this  figure  window.  The  number  of  sn^ishots  displayed  at  one 
time  is  only  limited  by  the  workstations  memory. 

Sampling  Frequency  Popupmenu 

Purpose 

Set  the  sampling  frequency. 

Use 

1.  Place  the  cursor  on  the  Sampling  Frequency  Popupmenu  and  click  the  left 
mouse  button  once  to  qien  the  menu. 

2.  Place  the  cutsot  on  the  sampling  firequency  and  click  die  left  mouse  button.  The 
popupmenu  contains  sampling  frequencies  of  1,  4096,  8000,  8192,  11024,  22048 
and  44096  Hz.  If  the  sampling  frequency  is  not  one  of  diose  mentioned  earlier,  the 
desired  value  can  be  entered  by  selecting  the  “User”  menu  item.  In  this  case,  an 
edit  box  i^pears  in  place  of  the  popupmenu.  Enter  the  desired  sampling  frequency 
into  this  edit  box  and  press  return. 

Result 

The  display  is  redrawn  with  a  new  time  axis  acox’ding  to  the  sampling  frequency. 

Notes 

In  the  popupmenu,  the  “User”  menu  item  is  rqilaced  by  the  sampling  firequency  if 
it  is  entered  using  the  edit  box.  To  subsequently  change  the  sampling  firequency 
after  using  the  edit  box,  the  last  sampling  frequency  entered  has  to  be  selected. 
This  is  always  the  last  menu  item  on  the  popupmenu.  If  a  sampling  frequency  of  1 


SPC  Toolbox  User’s  Guide 


50 


cx>iiA#0Ncxmnt0LS 


Hz  is  seleded,  tiie  horizontal  shows  the  vector  indioes. 

Hme/Spectrum  lUdiolHittoiis 
Purpose 

When  jx’ovided,  changes  the  vector  display  between  the  time  and  frequency 
domains. 

Use 

1.  Place  tte  cursor  on  the  desired  radiobutton  optitm  and  click  the  left  mtnise  but- 

ttm. 

Result 

The  display  is  redrawn  apprt^riately. 

Notes 

•  If  the  display  is  zoomed  before  switching  dcnnains,  the  display  is  still  zotHned 
after  switching  back. 

•  The  spectrum  displayed  is  generated  by  the  Signal  Processing  Toolbox  spectrum 
command  using  1024  points  (512  points  displayed). 
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COSWAVE,SlNWAVE _ 

Purpose 

Generate  a  baseband,  cosine  or  sine  wave. 

Synopsis 

y  >  ooawave(lb); 

y  >  ooewave^duralion): 

y  •  oocwave(ft)^duration.fB); 

y  -  ooawave(1b.dunition.1a,phaaMhilt); 

y  >  ainwave^): 

y  •  sinwave(tb,duralion); 

y  -  sinwave(lb.durition.(s); 

y  >  sinwaye(tb.dtralion,fS4)haseshtt); 

Description 

oo8wave(1b)  or  sinwave(t))  -  Generates  one  second  of  a  sinusoidal  wave  at  cycle  ire* 
quency,^,  and  a  sampling  frequency  of  8192  Hz. 

ooswave(fb.duralion)  or  sinwawe(fb,duration)  •  Generates  duration  seconds  of  a  sinuso* 
idal  wave  at  cycle  frequency,  fb^  and  a  sampling  frequency  of  8 192  Hz. 

ooswave(fb,(luration.fs)  or  sinwave(lb,duration.fs)  >  Generates  duration  secmids  of  a 
sinusoidal  wave  at  cycle  frequency,^,  and  sampling  frequency, /r. 

oo6W8ve(fb.dur«lion,fs,pha8e8hift)  sinwave(ft),cluration.fsj)hase6hif0  •  Generates 
duration  seconds  of  a  sinusoidal  wave  at  cycle  frequency,  ft),  and  sampling  fre- 
quency,/5.  with  phase  shift,  phaseshift. 

Examples 

Use  a  10  Hz  sinusoidal  modulating  signal  to  generate  an  AM  signal  with  a  modula¬ 
tion  index  of  0.5  and  a  carrier  frequency  of  SO  Hz,  sampled  at  128  Hz. 

y  >  ooswave(10, 128,60); 
yy-d8blo(y.0.5.50.128); 

Generate  3425  samples  of  a  800  cycle-per-second  sinusoidal  wave  sampled  at 
7000  Hz. 

y  -  ooswave(800,342SmXX).7000) 

See  Also 

sqwave,  triwave,  sawwave 
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Purpose 

Generate  a  double-sideband.  lai:ge-caiTier  amplitude  modulated  signal. 

Synopsis 

y  >  d8falc(modsig.m.tc.ls) 

Description 

d8blc(modsig.m.fc.f8)  -  Generates  a  double-sideband,  laige-carrier  amplitude  modu¬ 
lated  signal,  y.  The  signal  contaiimd  in  the  vector  modsig  is  modulated  tmto  a  sinu¬ 
soidal  carrier  of  amplitude  1  at  carrier  frequency,  fc.  The  maximum  deviatitm  of 
the  modulating  signal  must  vary  evenly  about  zero  [i.e.  (max(modsig)  «>  ab6(min(- 
modsiQ))].  modsig  is  multiplied  by  the  modulation  index,  /n.  before  being  modu¬ 
lated  onto  the  carrier.  The  length  of  the  signal  returned  is  equal  to  the  lengdi  of 
modsig.  The  carrier  is  generated  as  a  cosine  with  zero  phase  shift.  The  sampling 
frequency,  fs,  must  be  the  same  sampling  frequency  as  that  used  to  generate  the 
modulating  signal. 

Example 

Modulate  one  second  of  a  60  cycle  triangular  wave  using  a  carrier  of  1024  Hz  with 
a  modulation  index  of  0.3. 

X  « triwave(60.30); 

y-dsblc(x.0.3.1Q24): 

piottim6(y,0.04): 


Algorithm 

The  standard  equation  for  double-sideband,  suppressed-carrier  amplitude  modula¬ 
tion  with  zero  phase  shift 

y(r)  =  (1  +  mfit}  cos  {2nf^t) 

See  Also 

dsbac,  modindex,  envelope 
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DSBSC 


Purpose 

Geno^ate  a  dcmble-sideband,  suppressed-carrier  amplitude  modulated  signal. 

Synopsis 

y  •  d8bsc(mod8ig,fc.fa) 

Description 

d8bac(modsio.fc,fs)  -  Generates  a  double-sideband,  suppressed-canio'  amplitude 
modulated  signal,  y.  The  signal  ctmtained  in  the  vector  modsig  is  modulated  oaxo  a 
sinusoidal  carrier  of  amplitude  1  with  a  carrier  frequency  fc.  The  length  of  the 
modulated  signal  returned  is  equal  to  the  length  of  the  modng  vector.  Tlw  carrier  is 
generated  as  a  cosine  wave  wiA  zero  phase  shift.  The  sampling  frequency,  A  tnust 
be  the  same  sampling  frequency  used  to  generate  the  modulating  signal. 

Example 

Generate  16  bits  of  an  8  bit-per-second  signal  sampled  at  128  Hz.  Use  this  base¬ 
band  signal  as  the  modulating  sigiud  input  to  dsbac  to  generate  an  on-off  keyed  sig¬ 
nal. 

X  -  unipolar(8.I1  0  0 1  0 1  0 1  1  0  0 1 1  0 1].128); 
y«dsb6c(x,32.128): 


Tm«  domain  pM 


Algorithm 

The  standard  equation  for  double-sideband,  suppressed-carrier  amplitude  modula¬ 
tion  with  zero  phase  shift  is  given  by: 

.v(/)  s  m(Ocos(27i/j.O 


See  Also 

dsbac 
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ENVELOPE 


Purpose 

Demodulate  an  amplitude  modulated  signal  using  envelope  (noncoherent)  detec- 
dim. 

Synopsis 

[y.m]  -  onvelope(x); 

Description 

env6iope(X)  -  Returns  die  AM  envelc^  of  the  input  signal,  x,  and  the  modulation 
index,  m. 

Example 

Demodulate  using  envelope  detection  the  DSBLC  modulated  triangular  wave  used 
in  the  DSBLC  command  example. 

X  -  triwave(60,30): 
y  -  dsblc(x.0.3,1024) 
z  -  enveiope(y); 
plotlime(z,0.04); 


(Note  the  transients  due  to  the  FIR  mpiementation  of  the  Hiibeit  Transform.) 

Algorithm 

The  complex  envelc^  of  die  signal  x  is  obtained  by  computing  the  Hilbert  trans¬ 
form  of  X.  The  resulting  real  enveli^  is  the  magnitude  of  the  complex  envelope. 

See  Also 

dsbsc,  dsbic,  modindex 
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GRAFILTR 


The  Graphical  Filter  Design  Dodad  provides  an  interactive,  graphical  envircximent 
within  h^tlab  to  design  digital  filters  through  individual  (c<»nplex-c<mjugate  pair) 
placement  of  poles  and  zeros  rni  the  pole-zero  plot  of  the  filter  transfer  fiincticm. 
>^th  the  Graphical  Filter  Design  Dod^,  you  can: 

•  Place  complex,  real  or  purely  imaginary  poles  and  zeros  onto  the  pole- 
zero  plot. 

•  Place  complex,  real  or  purely  imaginary  poles  and  zeros  onto  the  pole- 
zero  plot  by  specifying  the  real  and  imaginary  parts. 

•  Place  complex,  real  or  purely  imaginary  poles  and  zeros  onto  the  pole- 
zero  plot  by  specifying  the  magnitude  and  angle  (degrees  (x  radians). 

•  Move  pole  and  zero  locations  using  the  mouse. 

•  Delete  poles  ot  zeros  from  the  design  using  the  mouse. 

•  Force  die  filter  transfer  function  to  minimum  or  maximum  phase  condi¬ 
tion. 

•  View  the  filter  frequency  and  phase  response  as  the  design  progresses. 

•  Generate  the  impulse  response  of  the  filter. 

•  Print  the  filter  parameters. 

•  Save  the  filter  transfer  function  coefficients. 


With  Mouse  With  Mouse  Phese  Condition  Megnitude  Scele 


Starting  the  Graphical  Filter  Design  Dodad 

The  Graphical  Filter  Design  Dodad  can  be  started  in  one  of  two  different  ways. 
The  first  one  is  to  start  the  grafitr  with  two  vectors  as  input  aiguments.  In  this  usage 
C‘grt>fiKr(b,a)”)>  the  aiguments  a  and  b  specify  the  coefficients  of  a  pre-existing 
transfer  function  following  Madab  filter  transfer  function  conventions.  The  argu¬ 
ments  a  and  b  are  the  same  as  those  return  by  any  of  the  Madab  Signal  Processing 
filter  design  programs  such  as  butter  or  chebyl .  See  the  discussion  f<x  the  Madab  fi- 
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ter  crmimand  for  t  cOTiplete  description  of  the  transfer  functitm  specificatioa.  Once 
a  pre-existing  transfer  functitm  has  been  loaded  into  Qrafllr,  poles  and  zeros  can  be 
deleted,  added  and/or  moved  about  in  order  to  “fine  tune”  or  otherwise  manipulate 
the  filter  response. 

The  second  way  to  start  the  filter  dodad  is  to  execute  the  “grafUr”  command  with¬ 
out  any  arguments.  This  effectively  starts  grafiltr  with  a  transfer  function  of  1  (no 
poles  (X  zeros).  Filter  design  then  start  from  a  “clean  slate’  by  adding  poles  and 
zeros. 

No  output  argtiments  are  supported.  Completed  filter  designs  must  be  saved  using 
the  Save  pushbutton  as  described  below. 

Note:  Due  to  a  unique  usage  of  the  Matlab  polar  command,  this  program  can  only 
be  used  with  Matlab  v.4.1  on  Sun  workstations. 

Graphical  Filter  Design  Dodad  Specific  Controls 

In  addition  to  controls  discussed  under  common  controls,  the  following  controls 
are  specific  to  the  Grtqihical  Filter  Design  Dodad  (grafitr). 

“With  Mouse”  X  and  O  Pushbuttons 

Purpose 

Place  a  pole  or  zero  on  the  pole-zo-o  plot. 

Use 

1.  Set  the  prt^rty  of  the  pole  or  zero  using  the  “With  Mouse”  property  radiobut¬ 
tons  describe  below. 

2.  Select  the  pole,  “X”,  or  zero,  “O”,  pushbutton.  After  selection,  the  cursor  turns 
into  a  CTOsshair  and  the  x-axis  label  of  the  pole-zero  plot  changes  to  “Mark  pole 
location  with  cursor”  or  “Mark  zero  location  with  cursor”  as  ^propriate. 

3.  Move  the  crosshair  to  the  desired  pole  or  zero  location  on  the  polar  plot  and 
press  the  left  mouse  button. 

Result 

The  pole  or  zero  is  placed  on  the  pole-zero  plot  and  the  frequency  response  is 
recomputed  and  displayed.  The  mouse  cursor  changes  back  into  an  arrow  and  the 
x-axis  label  is  cleared. 

Notes 

•  Poles  and  zeros  located  off  the  displayed  pole-zero  plot  can  be  entered  using  one 
of  the  coordinate  entry  groups  describe  below. 

‘‘With  Mouse”  Property  Radiobuttons 
Purpose 

Select  whether  poles  and  zeros  placed  on  the  pole-zero  plot  are  complex,  real  or 
purely  imaginary. 

Use 

1.  Place  the  mouse  cutsch'  on  the  desired  property  and  press  the  left  mouse  button. 
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Result 

The  radiobutton  is  selected  and  the  chosen  property  is  ^plied  the  next  tune  a  pole 
or  zero  is  placed  with  the  mouse. 

Notes 

Complex  and  purely  imaginary  poles  and  zeros  are  placed  as  complex-conjugate 
pairs.  Real  poles  and  zeros  are  placed  as  a  single  pole  or  zero  on  the  real  axis.  In 
placing  real  poles  and  zeros,  only  the  location  of  the  cursor  with  respect  to  the  real 
(horizontal)  axis  is  impm^t.  In  placing  imaginary  poles  and  zeros,  only  die  loca¬ 
tion  of  the  cwsor  with  respect  to  the  imaginary  (vertical)  axis  is  important. 

Rectangular  Coordinate  Entry  Group 

Purpose 

Enter  a  pole  or  zero  location  using  rectangular  (real  and  imaginary)  coordinates. 

Use 

1.  Enter  the  real  and  imaginary  parts  of  one  pole  or  zero  in  the  complex-conjugate 
pair.  Enter  only  the  real  or  only  the  imaginary  part  of  a  purely  real  or  imaginary 
pole  or  zero. 

2.  Select  the  group  “X”  pushbutton  to  place  a  pole  (pair).  Select  the  group  “O” 
pushbutton  to  place  a  zero  (pair). 

Result 

The  desired  pole  or  zero  is  added  to  the  pole-zero  plot  and  the  frequency  response 
is  recomputed  and  displayed. 

Notes 

•  Numeric  expressions  (i.e  “1/2”,  “3*2/5”,  etc.,)  can  be  used  in  the  edit  boxes. 

•  Blank  edit  boxes  are  assumed  to  be  zero. 

Polar  Coordinate  Entry  Group 

Purpose 

Enter  a  pole  or  zero  location  using  polar  (magnitude  and  angle)  coordinates. 

Use 

1.  Enter  the  magnitude  and  angle  from  the  pole-zero  plot  origin  of  one  pole  or  zero 
in  the  complex-conjugate  pair. 

2.  Select  the  group  “X”  pushbutton  to  place  a  pole  (pair).  Select  the  group  “O” 
pushbutton  to  place  a  zero  (pair). 

Result 

The  desired  pole  or  zero  is  added  to  the  pole-zero  plot  and  the  frequency  response 
is  recomputed  and  displayed. 

Notes 

•  Numeric  expressions  (i.e  “pi/4”,  “45/2”,  etc.,)  can  be  used  in  the  edit  boxes. 

•  Angle  entry  is  tied  to  the  Angle  Scale  checkbox.  If  the  Angle  Scale  checkbox  is 
checked,  angle  entry  is  in  degrees.  Otherwise,  angle  entry  is  in  radians. 

•  Blank  edit  boxes  are  assumed  to  be  zero. 

Phase  Condition  Radiobuttons 

Purpose 

Select  between  a  miiumum,  maximum  or  “mixed”  phase  transfer  function. 
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Use 

1.  Place  the  mouse  cursor  on  the  desired  condititm  and  press  the  left  mouse  butttm. 

Result 

A  new  transfer  function  is  computed  and  displayed  according  to  the  desired  selec¬ 
tion. 

Notes 

•  Minimiun  phase  systems  have  all  their  polynomial  zeros  located  inside  the  unit 
circle.  Maximum  phase  systems  have  all  their  polynomial  zeros  located  outside  the 
unit  circle.  “Mix«l”  phase  systems  have  zeros  located  both  inside  and  outside  the 
unit  circle. 

■  Selecting  Minimum  or  Maximum  phase  does  not  affect  the  internal  storage  of  the 
zero  locations.  Selecting  ‘^Mixed**  after  selecting  Minimum  or  Maximum  returns 
the  zeros  to  their  original  location. 

•  The  state  of  the  Phase  Condition  radiobuttons  affects  the  (^)eration  of  the  Print 
and  Save  pushbuttons. 

Magnitude  Scale  Radiobuttons 

Purpose 

Select  the  scaling  of  the  frequency  response  magnitude  plot  Gogarithmic  or  linear). 

Use 

1.  Place  the  mouse  cursor  on  the  desired  scale  and  press  the  left  mouse  button. 

Result 

The  frequency  response  if  redisplayed  with  the  desired  magnitude  scale. 

Angle  Scale  Checkbox 
Purpose 

•  Select  the  scaling  of  the  frequency  response  phase  plot  (degrees  or  radians). 

•  Select  the  entry  option  for  the  angle  e^t  box  in  the  Polar  Coordinate  Entry  group 
Use 

1.  Place  the  mouse  cursor  on  the  checkbox  and  select  or  unselect  as  desired. 

Result 

When  selected,  the  phase  response  is  displayed  in  degrees  and  the  entry  option  for 
the  angle  edit  box  is  in  degrees.  When  unselected,  both  of  these  are  in  radians. 

Impulse  Pushbutton 

Purpose 

Generate  a  plot  of  the  first  100  points  of  the  current  filter  impulse  response. 

Use 

1.  Place  the  mouse  cursor  over  the  pushbutton  and  press  the  left  mouse  button. 

Result 

A  child  graphics  window  is  opened  containing  a  plot  of  the  impulse  response. 

Notes 

•  The  filter’s  stability  can  be  determined  from  the  impulse  response. 

Print  Pushbutton 

Purpose 

Print  the  filter  parameters  into  the  Matlab  command  prompt  window. 
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Use 

1.  Place  the  mouse  cursor  over  the  pushbutton  and  press  the  left  mouse  button. 

Result 

Current  pole  and  zero  locations  are  printed  in  both  rectangular  and  polar  cowdi- 
nates  along  with  the  transfer  function  coefficients  in  the  Matlab  command  window. 

Notes 

•  The  state  of  the  Phase  Condition  radiobuttons  determines  whether  the  parameters 
of  a  minimum,  maximum  or  “mixed”  phase  system  are  printed. 

•  The  filter  parameters  can  be  saved  to  a  text  file  by  using  the  diary  command 
before  printing  or  by  cutting  from  the  conunand  window  using  the  cut  and  paste 
facility,  if  available. 

•  The  pole  or  zero  locations  can  be  saved  to  a  variable  by  first  saving  the  transfer 
function  coefficients  using  the  Save  pushbutton  and  then  using  the  Matlab  roots 
command  to  compute  the  pole  and  zero  locations  from  the  transfer  function  poly¬ 
nomials. 

Snapshot  Pushbutton 

Purpose 

Recreate  the  current  pole-zero  and  frequency  response  plots  in  a  new  graphics 
window  as  a  preliminary  step  to  obtaining  hardcopy. 

Use 

1.  Place  the  mouse  cursor  over  the  pushbutton  and  press  the  left  mouse  button. 

Result 

The  current  plots  are  reaeated  in  a  new  graphics  window.  The  user  can  then  treat 
this  window  as  any  other  Matlab  graphics  window.  Hardcopy  can  be  obtained 
using  the  Matlab  print  command. 

Delete  Pushbutton 

Purpose 

Delete  a  pole  or  zero  complex-conjugate  pair  using  the  mouse. 

Use 

1.  Place  the  mouse  cursor  over  the  pushbutton  and  press  the  left  mouse  button. 

2.  The  mouse  cursor  will  change  into  a  crosshair  and  the  x-axis  label  will  change  to 
“Click  cursor  on  item  to  delete.” 

3.  Move  the  crosshairs  over  one  of  the  poles  or  zeros  in  the  complex-conjugate  pair 
and  select  the  left  mouse  button. 

Result 

The  complex-conjugate  pair  is  deleted  and  the  transfer  function  is  recomputed  and 
displayed.  The  mouse  cursor  changes  back  into  an  arrow  and  the  x-axis  label  is 
cleared. 

Notes 

•  Real  poles  and  zeros  are  deleted  one-at-a-dme. 

Restore  Pushbutton 

Purpose 

Restore  the  original  filter  design. 


SPC  Toolbox  User’s  Guide 


60 


catAFiL'ni 


Use 

1.  Place  the  mouse  cunor  over  the  pushbutton  and  press  the  left  mouse  button. 

Result 

If  grafiltr  was  called  with  input  arguments,  the  transfer  function  specified  by  the 
calling  arguments  is  restored.  All  edits  are  lost  if  not  previously  saved  using  the 
Save  pushbutton.  If  no  input  arguments  were  used,  the  Restore  pushbutton  has  the 
same  effect  as  the  Clear  pushbutton. 

Clear  Pushbutton 

Purpose 

Clear  all  poles  and  zeros  from  the  filter  design,  effectively  setting  the  transfer  func¬ 
tion  to  equal  the  constant  1 . 

Use 

1.  Place  the  mouse  cursor  over  the  pushbutton  and  press  the  left  mouse  button. 

Result 

All  poles  and  zeros  are  cleared.  All  edits  are  lost  if  not  previously  saved  using  the 
Save  pushbutton. 

Save  Pushbutton 

Purpose 

Saves  the  filter  transfer  function  coefficients  to  the  given  variable  named  suffixed 
with  an  underscore  and  an  “a”  for  the  zero  polynomial  parameters  and  a  “b”  for  the 
pole  polynomial  parameters. 

Use 

1.  Place  the  cursor  on  the  Save  Pushbutton  and  click  the  left  mouse  button. 

2.  The  Save  pushbutton  is  replaced  by  an  edit  box.  If  Save  has  already  been  used, 
the  previously  entered  vector  name  is  contained  in  the  edit  box. 

3.  Place  the  cursor  on  the  edit  box  and  click  the  left  mouse  button.  Using  the  key¬ 
board,  backspace  over  any  previous  entry  and  then  type  in  the  desired  name  (fol¬ 
lowing  standard  Matlab  variable  name  conventions). 

4.  Press  the  return  key. 

Result 

The  transfer  function  is  saved  as  specified  above  to  the  Matlab  workspace.  The  edit 
box  is  replaced  by  the  Save  pushbutton. 

Notes 

Once  saved  to  the  Matlab  workspace,  these  variable  can  be  treated  like  any  other 
“a”  or  “b”  variable  returned  by  the  Matlab  Signal  Processing  Toolbox  filter  design 
functions. 

Close  Pushbutton 
Purpose 

Qose  the  dodad. 

Use 

1.  Place  the  mouse  cursor  over  the  pushbutton  and  press  the  left  mouse  button. 

Result 

The  dodad  window  is  closed  and  any  global  variables  o-eated  by  the  dodad  are 
cleared. 
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LD8BIT,LD16BITJ.D32BIT _ 

Purpose 

Load  a  file  of  samples  stored  as  8-bit  signed  integers. 

Synopsis 

y  -  Id8bil(‘name'): 
y  -  ld8t)il(‘nanM’,numbar) 
y  •  ld8bit(‘name'.number.ollset) 
y  >  USbitCname’.O.offset); 

Description 

USbitCname’)  -  Load  the  entire  file  as  8-bit  signed  samples. 

hJSbitCname'.number)  -  Loads  the  first  number  of  samples  starting  firtnn  the  begin¬ 
ning  of  the  file. 

IdSbitCname'.number, offset)  -  Loads  number  of  samples  beginning  at  offset  from  the 
beginning  of  the  file.  If  number  s  0.  loading  starts  from  offset  and  continues  to  the 
EOF. 

The  LD16BIT  and  LD32BIT  commands  take  the  same  arguments  as  the  LD8BIT 
command. 

Example 

Datafile.dat  contain  four  seconds  of  samples  sampled  at  8192  Hz.  Load  the  data 
collected  during  the  third  second. 

y  .  Id8bit(‘datafie.dar.2*81 92-1-1.8192); 

See  Also 

Uiebt.  I(i32bit 
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LOADSSPl,SAVESSPI _ 

Purpose 

Load  data  stored  in  the  fwmat  specified  by  the  Cyclic  Spectral  Analysis  Software 
Package  from  Statistical  Signal  Accessing.  Inc. 

Synopsis 

y  -  loadsspKViame'): 

y  -  loaclsspi('name’.*binary’); 

sayesspKx.'name’); 

Description 

loadsspi(*nama')  •  Loads  SSPI  software  generated  data  from  the  file  name.  Data  is 
stored  in  ASCII  fmmat. 

Ioadsspi('nafne','binary0 '  Data  is  stored  in  binary  format. 

8avesspi(x,’name’)  -  Saves  the  vector  x  to  the  file  name  in  SSPI  data  file  format. 
Data  is  written  in  floating  point  ASCII. 

savesspi(x,‘name', ’binary*)  -  Data  is  stored  in  binary  format. 

ASCn  SSPI  format  is: 

first  line  -  type  nbrsamples 

rest  of  lines  -  one  or  two  columns  of  ascii  numeric  data 


where, 

type:  1  s  real,  2  =  complex 

nbrsamples  =  number  of  samples  (one  sample  per  line). 


2  6 

3.226078e-01  2.978590e+00 
5.433282e-01  8.382938e-01 
3.492872e+00  9.382934e+00 
3.2938426-01  2.394829e+00 
4.5928396+01  2.2939406+00 
6.2034826+00  2.3829206+01 

Binary  format  follows  the  same  line  except  data  is  stOTe  serially. 

Reference 

[1]  Stephan  V.  Schell  and  Chad  M.  Spooner,  Cyclic  Spectral  Analysis  Software 
Package  User’s  Manual,  Statistical  Signal  Processing,  Inc.,  1991. 


SPC  Toolbox  User’s  Guide 


63 


LOADVOC^AVEVOC 


LOADVOC,SAVEVOC _ 

Purpose 

Load  and  save  data  to  and  from  Soundblaster  Voice  files. 

Synopsis 

(y.fa]  -  loadvoc(name): 
sav«voc(nama,x,f8); 

Description 

[y,f4  -  ioadvoc(name)  -  Loads  die  data  samples  from  a  Soundblaster  Creative  Voice 
Pile  (*.voc)  into  y  and  stores  the  sampling  frequency  into  fs.  The  function  checks 
the  file  to  ensure  it  is  a  versimi  1.10  Soundblaster  file.  This  function  only  supports 
raw.  8-bit.  unpacked  voice  files.  The  filename  extensimi  “.voc**  is  $ppeaded  to 
name  if  no  extension  is  supplied.  Any  mean  is  removed  after  being  read  frtnn  the 
voice  file. 

savevoc(name,x.fs)  •  Saves  the  data  in  x  to  a  Soundblaster  Creative  Voice  File 
(*.voc)  named  “name.voc”.  Default  sampling  frequency  is  8192  Hz.  This  function 
saves  the  file  as  a  version  1.10  file  and  stores  the  data  in  the  raw.  8-bit  unpacked 
format. 

Example 

Load  a  voice  file,  cut  the  sampling  ratt  in  half  and  store  to  a  new  file; 

[s,f^  •  kMdvocCaeaisiO: 
y  ■  decimate(s,2); 
savevocCSiowfsear.y,fs/2); 

See  Also 

loadwave,  savewave,  readau,  writeau 

Reference 

[1]  The  Developer  Kit  for  Sound  Blaster  Series,  pp  4-6  to  4- 1 2,  Creative  Labs,  Inc., 
Nov  1991. 
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Purpose 

Diq>lay  the  periodognun  of  •  signal  on  a  logihthmic  scale. 

Synopsis 

gienonn(x} 

iperigrm(x.fB) 

lperionn(x.niindb) 

lperionm(x.fB,mindb) 

lpenonn(x.fs,‘phaBe*) 

lporignn(x,f6,nindb.‘phaae’) 

Description 

lperignn(x)  -  Computes  the  periodpgram  of  x  and  displays  die  result  on  a  logarith¬ 
mic  (decibel)  scale.  Only  the  p<»itive  frequencies  are  displayed.  The  default  sam¬ 
pling  frequency  is/s  -  8192  Hz. 

lperignn(x.fB)  -  Sets  the  sampling  frequency  to  fs. 

Iperigrm(x,mindb)  •  Cuts  the  plot  off  below  mifidb  for  a  better  display  of  the  data 
above.  If  mindb  is  positive,  it  is  interpreted  as  /r  in  the  above  usage.  Therefore, 
mindb  must  be  negative  for  this  usage. 

lperiorm(x,fe.mindb)  •  Used  to  set  fs  and  mindb  to  other  than  their  default  values  or  to 
set  mindb  when  mindb  is  positive. 

lperiorTn(x,fe,'phase')  and  lperiomi(x,fs.iTiindb, ’phase’)  -  Splits  the  graphics  window 
and  displays  both  the  magnitude  and  phase  information. 

Example 

Generate  a  frequency-shift  keyed  signal  and  display  it’s  periodogram.  Display  the 
magnitude  down  to  -1(X)  dB. 

y  -  bfsk(S12, 512.2048,0.1); 
lperigrm(y,-100); 
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Algorithm 

The  periodogram  is  given  by: 

where  is  the  fourier  transform  of  i(r) : 

t*T/2 

f  x(u)e'^^*‘*du 

t-T/2 

The  plot  magnitude  is; 

201pg(|X^r./)|) 

Limitations 

The  size  of  the  PPT  used  to  compute  the  fouri^  transform  is  die  closest  power-of> 
two  greater-than  or  equal-to  the  length  of  x.  Long  input  vectors  require  kmger 
compute  times.  If  jc  is  short  and  a  finer  firequency  resolution  is  desirul.  the  iiqiut 
vector  X  can  be  zero  padded  befcve  calling  the  function. 

Iperigrm([x  zero6(300.1)]); 

See  Also 

wperigrrn.  plotlime 

Reference 

[1]  Williaffl  A.  Gardner,  Statistical  Spectral  Analysis,  A  Nonprobabilistic  Theory, 
pp.  5-7,  Prentice-Hall,  1988. 
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LRS 


Purpose 

Generate  a  maximal-length,  linear  recursive  sequence. 

Synopsis 

y  -  lis(r.n); 
y  -  Ii8(r,n,fl,tape): 

Description 

lrs(r.n)  -  Returns  n  bits  of  a  binary,  maximal-length,  linear  recursive  sequence  using 
an  r  bit  shift  register.  Uses  a  randmnly  generated  fill  with  randmnly  genmted  taps. 

Ir8(r,n,fl.tap8)  -  Sets  the  initial  fill  and  ti^  locations  to  those  specified  in  the  vectors 
fill  and  taps.  These  vectors  may  be  q)ecified  by  a  vector  cmitaining  a  binary  r^e- 
sentation  of  the  register  and  the  taps  or  by  a  list  of  locations  for  each  cell  set  to  one 
and  each  t^)  location.  Pw  example,  an  initial  fill  of  [1 0  1 0  1]  can  also  be  specified 
by  [1  3  5]  and  taps  of  [1  0  0  1 0]  can  be  specified  as  [1  4].  If  the  binary  representa¬ 
tion  is  specified,  its  length  must  equal  r.  If  a  list  is  specified,  it  largest  element  must 
be  less  than  r.  The  ouq)ut  is  taken  from  a  tap  off  cell  one.  Therefore,  a  tap  off  stage 
one  is  always  included  whether  not  it  is  specified. 

Example 

Generate  IS  bits  of  an  RS(1,4)  LRS  with  an  initial  fill  of  all  ones. 
y-lrs(5.15.I1 1 1 1 11.W); 

y- 

[1001  10100100001] 


[T 

2 

4 

3 

MOD  2 
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Purpose 

Compute  the  filter  coefficients  for  a  least-squares  optimal  filter  using  singular 
value  decmnposititm. 

Synopsis 

[h.S.Px]  -  ar_svd(d.x,P) 

Description 

arjsvd(d.x.P)  *  Given  the  observed  data  sequence  x,  compute  the  filto*  coefficients 
for  a  least-squares  q)timal  filter  of  order  P,  that  produce  the  desired  data  sequence 
d  firom  the  observed  data  sequence.  Ouq>ut  arguments  are: 

h  =  filter  coefficients 
S  =  sum  of  squared  errors 
Px  =  the  projection  matrix 

Example 

Given  the  desired  data  sequence  [1  -2  3  -4  5],  compute  the  filter  coefficients  of  the 
second  (xder  least-squares  filter  that  will  produce  the  desired  data  sequence  [1-11 
-1  1].  (Example  9.4  firom  [1]) 

d  -  (1  -1 1  -1 11; 
x-11  -23-451: 

[h.Si  -  ls_8vd(d.x.2); 
h-  1.0000 
1.0000 

S-  3.5527e-015 
y-f«er(1,h,[1  00001); 

y 

1-11-11 

See  Also 

Is.svd 

Reference 

[1]  Charles  W.  Theirien,  Discrete  Random  Signals  and  Statistical  Signal  Process¬ 
ing,  pp  518-528,  Prentice-Hall,  1992. 
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LS.WHOPF 


Purpose 

Compute  the  filta*  coefficients  for  a  least-squares  optimal  filter  using  the  Menor- 
Hopf  equation. 

Synopsis 

[h.S,Px]  -  ar_«vtK)pf(d.x,P) 

Description 

ar_whopf(d,x.P)  -  Given  the  observed  data  sequence  x,  compute  die  filter  co^- 
dents  for  a  least-squares  (^timal  filter  of  order  P,  that  will  {M-oduce  die  desired  data 
sequence  d  from  the  observed  data  sequence.  Ouq>ut  aiguments  are: 

h  =  filter  coefficients 
S  -  sum  of  squared  eiros 
Px  =  the  projection  matrix 

Example 

Given  the  desired  data  sequence  [1  -2  3  -4  5],  compute  the  filter  coefficients  of  the 
second  (xder  least-squares  filter  that  will  produce  the  desired  data  sequence  [1-11 
-1  1].  (Example  9.3  from  [1]) 

d  -  [1  -11  -1 IJ; 

X-I1  -23-451; 

[b.S5  ■  ls_whod(d,x,2): 
h-  1.0000 
1.0000 
S-  0 

y-fiter(1.h,[1  0000]); 

y- 

1-11-11 

Algorithm 

Wiener-Hopf  equation  solved  for  h: 

h  =  (X^'^X) 

Sum  of  squared  errors: 

S  *  d*''d-d*^Xh 

See  Also 

ISJBVd 

Reference 

[1]  Charles  W.  Theirien,  Discrete  Random  Signals  and  Statistical  Signal  Process- 
PP  518-528,  Prentice-Hall,  1992. 
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MINPHASE 


Purpose 

Return  a  polynomial  dut  is  in  minimum-phase  form. 

Synopsis 

y  >  minph88e(P); 

Description 

minpha8e(P)  -  Returns  a  minimum-phase  form  polynrMnial  for  F  if  P  is  not  mini¬ 
mum-phase  already.  If  P  already  is  minimum  phase.  P  is  simply  returned.  A  mini¬ 
mum-phase  polynomial  is  one  whose  roots  are  inside  the  unit  circle. 

Example 

Return  the  minimum  phase  form  of  2z^-«-z-t-4.  The  roots  of  this  polynomial  are: 

It 

r,  2  -  1.414« 

The  roots  of  the  resulting  minimum  phase  polynomial  are  : 

ft 

r,  j-  0.7071e  V 

It  can  be  verified  using  the  fireqz  command  that  the  magnitude  of  the  impulse 
response  for  both  these  polynomials  are  the  same  and  that  the  second  polynomial 
has  minimum  phase  lag. 

b  -  minphase(I2 1  4]); 
b-  4.0000 
1.0000 
2.0000 

[h.w]-freqz(b,1.50); 

mao  "  phase  «  angleCh); 

8einilooy(w,mag);  semilooy(w,phase); 

Polynomial  b  is  4z^ + z + 2 . 

Minmuni  phiM  Mnnum  phast 
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Algorithm 

The  fijnrithm  multiplies  the  polynomud  by  the  following  factor  for  each  root  diat 
has  a  magnitude  greater  than  1: 

ZZq  >  1 


See  Also 

oonv.  deoonv,  freqz,  poly 


References 

[1]  Charles  W.  Thenicn,  Discrete  Random  Signals  and  Statistical  Signal  Process¬ 
ing,  pp  250-253,  Prentice-HaU.  1992. 

[2]  Alan  Y  Oppenheim  &  Ronald  W.  Schafer,  Digital  Signal  Processing,  pp  345- 
353,  Prentice-HaU,  1975. 
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MODINDEX 


Purpose 

Computes  the  modulation  (deviation)  index  of  an  AM  signal. 

Synopsis 

y  -  modindex(x); 

Description 

modindex(x)  -  Returns  the  modulation  (deviation)  index  of  x. 

Example 

Return  the  deviation  index  of  a  DSBLC  modulated  triangular  wave  used  in  the 
DSBLC  command  example. 


X  -  triwave(60,30); 
y  -  dsblc<x.0.3,1024) 
m  -  modindex(y); 

m« 


3.0000 


Algorithm 

Finds  the  complex  envelope  and  from  its  extrema  computes; 

max  -  min 

m  = - 

max  +mm 


See  Also 

dsbsc,  dsbic,  envelope 
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NORMALEQ _ 

Purpose 

Solves  a  system  of  Normal  equations  used  in  linear  predictive  filtering. 

Synopsis 

[e.a]  >  normaleq(R) 

Description 

nomialeq(R)  -  Returns  the  prediction  error  variance,  e.  and  the  linear  predictive  fil¬ 
ter  coefficients,  [i  a,  aj  - 

Example 

Use  nonrnaleq  to  solve  the  system  of  nmmal  equations  used  in  Example  7.1  of  [1]. 

R  >  toeplitzdl  0.5  0.2^); 

[e.a]  -  nonnaleq(R); 

e- 

0.7S00 

a« 

-0.5000 

0 

Algorithm 

This  function  is  a  direct  implementation  of  equation  (7.23)  given  in  [1].  Refer  to 
[1]  for  a  detailed  discussion.  In  references  on  linear  predictive  filtering,  these  equa¬ 
tions  are  often  referred  to  as  the  augmented  normal  equations. 

Normal  equations  for  p  =  2: 


1 

a"! 

/!*,(-!]  «,[0]  /?,[!) 

= 

£  1 
o| 

Loj 

where  o*  is  the  prediction  error  variance  and  [n]  is  the  auto-correlation  of  x  at  n. 

Reference 

[1]  Charles  W.  Therrien,  Discrete  Random  Signals  and  Statistical  Signal  Process¬ 
ing,  eq  7.23,  p  345,  Prentice-Hall,  1992. 
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OOK,OOKMSG 


Purpose 

Generate  an  on-off  keyed,  bandpass  signal. 

Synopsis 

y  -  ool<(Rb.fc); 
y  -  ook(Rb.fc,duralion); 
y  -  ook(Rb,fc,duration.fs); 
y  -  oolvn8g(Rb.fc.m8g); 
y  -  ooton8g(Rb,fc.f8.m8g); 

Description 

ook(Rb,fc)  -  Generates  one  second  of  an  on-off  keyed  signal  with  a  bit  rate  of  Rb 
bits-per-second  centered  at  a  carrier  frequency,  fc.  Default  sampling  rate  is  8192 
Hz. 

ook(Rb,fc.duration)  -  Generates  duration  seconds  of  the  on-off  keyed  signal. 
ook(Rb,fc,duration,fs)  -  Set  the  sampling  frequency  to  fs. 

ookmsg(Rb,fc,msg)  or  ookmsg(Rb,fc,fs,msg)  •  Generate  a  bandpass  binary  on-off 
keyed  signal  with  a  bit  rate  of  Rb  bits-per-second,  at  a  carrier  frequency, /c,  sam¬ 
pled  at  sampling  frequency, /s,  and  containing  the  message,  msg.  Mfs  is  not  given, 
fs  defaults  to  8192  Hz. 

Example 

Generate  a  64  bit-per-second,  on-off  keyed  signal  at  a  carrier  frequency  of  2048  Hz 
lasting  one-and-one-half  seconds.  Play  the  signal  over  the  workstation’s  speaker. 

y  -  ooK(64.2048,1.5); 
sound(y,8192); 

Algorithm 

The  on-off  keyed  signal  is  generated  by  first  creating  a  unipolar  signal  at  the 
desired  bit/sampling  rates,  and  then  applying  a  double-sideband,  suppressed-car- 
rier  modulation  to  the  resulting  b^eband  signal  at  the  desired  earner  fluency. 

Limitations 

See  limitations  under  unipolar. 

See  Also 

unipolar,  dsbac,  bpak,  bpakmsg,  bfsK  bfskmsg 

Reference 

[1]  Leon  W.  Couch,  Digital  and  Analog  Communication  Systems,  pp.  332-335, 
Macmillan,  1990. 
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PEAKS 


Purpose 

Find  the  peaks  within  a  vecttv. 

Synopsis 

M-peaks(x); 

[y.il-peaks(x); 

Description 

PEAKS(x)  -  Returns  a  vector  containing  the  magnitude  of  the  peaks  in  vector  x.  A 
"peak”  is  defined  as  a  vector  member  that  is  either  greato-  or  smaller  than  both  of 
its  immediate  neighbors  depending  upon  whether  it  is  a  maximtim  or  minimum 
"peak”.  Optional  output  argument  i  returns  a  vector  cttresponding  to  the  indices  of 
the  peaks  returned  in  y.  If  the  ou4)ut  aigument  /  is  not  requested,  ouq>ut  argument  v 
will  be  the  same  length  as  x  with  non-peak  members  set  to  zero.  If  ou^ut  aigument 
/  is  requested,  output  argument  y  will  contain  only  the  peak  values.  There  is  a  one- 
to-one  cmrespondence  between  members  of  y  and  i. 

Examples 

Find  the  peaks  of  [1  5  4  2  6  7  -1  0  3]. 

x-[1  54267-1  031: 
y  -  peaks(x); 
y-[050207-103] 

Return  just  the  peaks  and  their  indices. 

|y.i]  -  peaks(x); 
y- [52  7-1  31 
i-|246791 

Limitations 

PEAKS  makes  extensive  use  of  for  loops  and  therefore  will  be  slow  when  the 
input 
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PLOTTIME 


Purpose 

Plot  the  time-domain  display  of  a  signal. 

Synopsis 

piottime(x) 

plottime(x,f8) 

plottime(x.duration) 

plottiine(x.fs.duraten) 

plottime(x,duration,on8et) 

plottime(x,fs.duration,offset) 

Description 

piotlime(x)  -  Plots  the  time-domain  representation  of  x  with  a  time  axis  scaled  fcx  a 
sampling  frequency  of  8192  Hz.  Default  duration  is  one  second. 

plottime(x.f8)  -  Sets  the  sampling  frequency  to  fs  Hz.  The  variable  fs  must  be  greater 
than  or  equal  to  100  Hz  for  this  usage. 

plottime(x, duration)  -  Plot  up  to  duration  seconds  of  x.  The  duration  must  be  less 
than  100  seconds  for  this  usage. 

plottime(x,fs,duration}  •  Plots  up  to  duration  seconds  of  x. 

pk)ttim6(x,(hiration.offset)  and  piottime(x,fs, duration, offset)  -  Starts  the  plot  offset  sec¬ 
onds  from  zero  seconds. 

Example 

Generate  one  second  of  a  binary  phase-shift  signal.  Display  1/40  seconds  of  the 
signal  beginning  at  0.3  seconds. 

y  -  bpsk(256.256); 
plottime(y,1/40,.3); 


See  Also 

iperigrm,  wperigmi 
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SANALYZR 


The  Spectrum  Analyzer  Dodad  provides  a  gr^>hical,  interactive  environment  from 
which  to  take  spectral  measurements  of  a  signal.  Its  operation  is  similar  to  that  of 
the  lab  equipment  of  the  same  name.  IVo  independent  curstM^  are  provided  along 
with  a  digiud  readout  of  their  ct^esponding  frequencies  and  magnitudes  and  the 
difference  in  frequency  and  magnitude  between  them.  The  curs(»^  may  be  moved 
from  sample-to-sample  using  the  or  pushbuttons  or  from  peak-to*peak 
(high  or  low)  using  the  “«“  w  “»”  pushbuttons.  Single  “clicking”  on  the  spec¬ 
tral  line  with  the  mouse  cursor  moves  the  closest  spectrum  analyzer  cursor  to  the 
position  of  the  mouse  cursor.  A  spectrum  analyzer  cursor  can  also  be  grabbed  with 
the  mouse  and  dragged  to  a  new  position.  The  spectrum  analyzer  curstx^  can  also 
be  moved  to  an  specified  frequency  be  entering  the  frequency  in  the  cursor’s  fre¬ 
quency  edit  box  and  pressing  return. 

Starting  the  Spectrum  Analyzer  Dodad 

The  Spectrum  Analyzer  Dodad  can  be  started  in  two  different  ways.  The  first  is  to 
start  the  dodad  with  a  vector  name  as  an  argument.  A  vector  is  defined  as  a  IxN  or 
a  Nxl  Matlab  variable  and  it  must  exist  in  the  M  lab  workspace  before  starting 
the  Spectrum  Analyzer  Dodad.  To  analyze  a  vector  named  “myvector”,  the  com¬ 
mand  “sanalyzr(  myvector)”  needs  to  be  executed. 


No  ouq)ut  arguments  are  supported. 


Notes 

•  The  spectrum  displayed  is  the  output  of  an  enhanced  version  of  the  spectrum 
command  from  the  Matlab  Signal  Processing  Toolbox  named  spectrm2.m  which 
allows  for  user-defined  windows. 
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•  The  default  window  is  still  a  hanntng  window  but  the  window  can  be  changed 
after  sanalyzr  is  started  by  selecting  another  window  from  the  Window  pulldown 
menu. 

•  The  spectrum  command  uses  the  Welch  method  of  spectral  estimation  and  is  used 
here  with  no  window  overlap.  The  default  FFT  length  is  1024  points.  Only  the  pos¬ 
itive  frequencies  up  to  half  the  sampling  frequency  are  displayed.  In  addition,  die 
DC  component  is  not  displayed. 

•  Movement  of  the  spectrum  analyzer  cursors  with  the  mouse  under  the  Sun  cqier- 
ating  system  requires  holding  the  mouse  button  down  until  the  cursor  movement  is 
complete.  Also,  after  dragging  a  cursor  with  the  mouse,  the  cursor  must  be  held 
down  until  the  cursor  position  stabilizes.  This  applies  to  both  versions  4.0a  and  4. 1 
of  Matlab  running  under  Sun  Open  Windows.  Under  Microsoft  Windows,  a  quick 
click  is  all  that  is  required  to  move  the  closest  spectrum  analyzer  cursor  to  the 
mouse  cursor  location  and  no  pause  is  required  after  dragging  a  spectrum  analyzer 
cursor  to  a  new  position. 
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SAWWAVE 


Purpose 

Generate  a  baseband,  sawtoodi  wa^.  This  function  can  also  be  used  to  generate  a 
ramp. 

Synopsis 

y  ■  sawwaveffb): 

y  -  8aNVwave((l>.’andpodan: 

y  >  savAvave(fb,duration): 

y  >  saMAvave(1b,duration, ’antipodal*); 

y  -  saMVwava(fb.duraiion,fs); 

y  >  8aMvwave(1b.duration.fs„‘antipodal'); 

Description 

sawwaveffb)  -  Generates  one  second  of  a  baseband,  sawtooth  wave  at  cycle  fre¬ 
quency,^.  and  sampling  frequency  of  8192  Hz. 

sawwaveffb, duration)  -  Generates  duration  seconds  of  a  baseband,  sawtooth  wave  at 
cycle  frequency,^,  and  sampling  frequency  of  8192  Hz. 

sawwave(fb,duralion,fs)  -  Generates  nbrcycles  of  a  baseband,  sawtooth  wave  at  cycle 
frequency,^,  and  sampling  frequency, /s. 

The  argument  ‘antipodal’  changes  the  amplitude  range  of  the  output  wavef(»rm 
from  [0.1]  to  [-1.1]. 

Examples 

Generate  30  cycles  of  a  sawtooth  wave  with  a  cycle  frequency  of  25  Hz  sampled  at 
8000  Hz,  varying  between  0  and  5  volts. 

X  -  5  *  sawwave(25, 30/25,8000); 
plottimefx, 8000, 0.2); 


T«n«  uom«n  pMC 
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Generate  a  ramp  frcmi  -1  to  3  with  100  samples: 
y-4*sawwave(1.1,100)<1; 

Limitations 

The  waveform  is  created  by  first  creating  a  single  cycle  of  the  wavefwm  and  then 
repeating  this  cycle  the  required  number  of  times.  Hence,  when  the  cycle  rate  is  not 
evenly  divisible  into  the  sampling  frequency,  truncation  errors  occur.  The  length 
of  the  ouq>ut  vecto*,  y,  can  be  computed  by  &e  formula: 

length  ■  flooriduration  floor 

The  exact  cycle  frequency  can  be  computed  from: 


See  Also 

sqwave,  triwave,  coswave.  sinwave 
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SETSNR 


Purpose 

Mix  two  vectors  (one  rq>resenting  a  signal,  one  representing  noise),  such  that  die 
ou^ut  signal’s  total  bandwidth  signal-to-noise  ratio  is  set  to  die  specified  value. 

Synopsis 

[y.sigfnid  *  setsnr(8ignal,SNR): 

[y.sioma]  >  8etsnr(signal,noise,SNR); 

Description 

SETSNR(signal.noi8e,SNR)  -  Returns  the  additive  signal  such  that  the  total  band¬ 
width  signal-to-noise  ratio  is  equal  to  SNR  decibels.  The  term  “total  bandwidth” 
implies  the  noise  power  in  all  frequencies  from  0  to  fsl2  is  used  in  the  computation 
of  the  SNR.  Sigma  is  the  parameter  a  described  below. 


SETSNR(signal,SNR)  -  Gaussian  white  noise  is  added  by  the  function. 

Example 

Create  a  frequency-shift  keyed  signal  and  then  add  Gaussian  white  noise  such  that 
the  total  bandwidth  SNR  is  15  decibels. 


X  -  bfsk(64.2048); 
y  -  setsnr(x,randn(length(x).1),15); 


Algorithm 


The  ouqiut  signal  is  defined  by 

y  [/i]  =  5  [/?]  +  a  •  Q)  [n] 
for  which  a  is  computed  such  that 


SNR  :=  10  log 


Var{s[n])  ~ 
a^V’ar((o[/i] ) 


where  SNR  is  the  desired  signal-to-noise  ratio.  The  amplitude  of  the  noise  is 
adjusted  before  the  noise  is  added  to  produce  the  desired  SNR.  In  the  case  of  Gaus¬ 
sian  white  noise,  relationship  between  a  and  the  noise  power  spectrum: 


^0 

2 


See  Also 

setsnibw 

Reference 

[1]  Leon  W.  Couch.  Digital  and  Analog  Communication  Systems,  pp.  103-105. 
Macmillan.  1990. 
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SETSNRBW 


Purpose 

Mix  two  vectors  (one  representing  a  signal,  one  representing  noise),  such  that  Ae 
ou^ut  signal’s  noise-equivalent  (in-band)  signal-to-noise  ratio  is  set  to  the  speci¬ 
fied  value. 

Synopsis 

[y, sigma]  -  setsnft)w(6ignal,SNR,fc.bw); 

[y,8i(^j  -  8etsnrbw(8ignal, noise, SNR,fc.tMtf,fB); 

Description 

SETSNRBW(sional.noise,SNR,fc,bw,fs)  -  Returns  the  additive  signal  such  that  the 
noise-equivalent  bandwidth  (in-band)  signal-to-noise  ratio  is  equal  to  SNR  deci¬ 
bels.  The  term  “noise-equivalent  bandwidth”  implies  the  noise  power  only  in  fre¬ 
quencies  within  the  signal  bandwidth  is  used  in  computing  the  SNR.  The  signal 
bandwidth  is  centered  at  4  and  is  equal  to  2*bw.  This  definition  requires  the  value 
actually  supplied  for  bw  be  one-half  the  desired  bandwidth  specification  (see 
example).  The  sampling  frequency,  fs,  is  required  to  find  the  power  spectral  den¬ 
sity  of  the  noise.  Sigma  is  the  parameter  a  described  below. 

SETSNRBW(signal,SNR,fc,bw)  -  Gaussian  white  noise  is  added  by  the  function  and 
the  sampling  frequency  defaults  to  8192  Hz. 

Example 

Create  a  coherent,  frequency-shift  keyed  signal  and  then  add  Gaussian  white  noise 
such  that  the  noise-equivalent  bandwidth  SNR  is  15  decibels.  Use  the  distance 
from  the  center  frequency  to  the  first  null  as  one-half  the  bandwidth. 

X  -  bpsk(256,2048); 

y  -  setsnrbw(x,15, 2048,256.81 92); 

lperigrm(y); 


Algorithm 


The  ouq)ut  signal  is  defined  by 


y[n]  a  5[n]  +a-  Q)[n] 
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for  which  a  is  computed  such  that 


SNR  - 10  log 


Var(5[yi]) 
a  [n] ) 


where  SNR  is  the  desired  signal-to-noise  ratio.  The  amplitude  of  the  noise  is 
adjusted  before  the  noise  is  added  to  produce  the  desired  SNR.  In  the  case  of  Gaus¬ 
sian  white  noise,  the  relationship  between  a  and  the  power  spectral  density  of  the 
noise  is: 


See  Also 

setsnr 


Reference 

[1]  Leon  W.  Couch,  Digital  and  Analog  Communication  Systems,  pp.  103-105, 
Macmillan,  1990. 
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SP_STENG,SP,STMAG,SP,STZCR _ 

Purpose 

Compute  the  short-time  energy,  magnitude  and  zero-crossing  curves  of  qjeech  sig¬ 
nals. 

Synopsis 

[y.Q  •  sp_steng(x,franM,overlap,fs,\vindow'); 

[y.Q  -  8p_slmaQ(x,frame,overt^,fe,>Mindo«/): 

[y.Q  *  spJStzcr(x,frame.overlap,fs,^Mndow’): 

Description 

SP_STENG(x,frame,overlap,fs,'vMndo«0  -  Compute  the  short-time  energy  of  x  using  a 
window  size  of  frame  length  and  a  percentage  overlap  between  successive  win¬ 
dows  using  a  ‘Svindow”  prefiltering  data  window.  Available  windows  are  ’rectan¬ 
gular’,  ’hamming’,  ’harming’,  ’blackman’  and  ’bartlett’.  A  rectangular  window  is 
the  default  if  the  ‘^window”  argument  is  not  given.  The  ou^ut  arguments  are: 

y  -  short-time  fimcuon  curve. 
t  -  time  indices  for  each  cmresponding  sample  in  y. 

SP_STMAG(x, frame, overlap, fs, ’window*)  -  Compute  the  sh(»t-time  magnitude. 
SP_STMAG(x,frame,overtap,fs,’window')  -  Compute  the  short-time  zero-crossing  rate. 

Examples 

Compute  and  display  the  short-time  energy  using  a  30  millisecond  frame,  a  20% 
frame  overlap,  a  Hamming  window  and  a  median  smoothing  filter.  Note  the  st^s 
necessary  to  plot  the  grtqjhs  on  the  same  scale  and  aligned  in  time. 

[y,f|  -  sp_8teno(x,0.03,20,fs,'hamming’}: 

yy  ■  mdsmoothfy,?); 

subplot(2,1 ,1  );plol((0:length(x)-1  )/fs,x); 

subplot(2,1,1);plot(t,yy) 


0  0.2  0.4  0.0  00  1  1.2  14  1.6 
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Algorithm 

Short-time  eneigy: 


^li  =  \[x{m)w(n-m)]f 

m 

Short-dme  magnitude: 

m 

^n-  |ac(/n)|w(n-m) 

m  ■-“ 

Short-time  zero  crossings: 

Z„  =  |sgn[j:(m))l -sgn[jf(m- l)]lw(n-m) 

See  Also 

avsmooth,  mdsmootti 

Reference 

L.  R.  Rabiner  &  R.  W.  Schafer,  Digital  Processing  of  Speech  Signals,  pp  120-130, 
Prentice  Hall,  1978. 
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SQWAVE _ 

Purpose 

Generate  a  baseband,  unipolar  [0.1]  (x  antipodal  [-1.1]  square  wave. 

Synopsis 

y  >  sqwave(fb); 

y  -  sqwave(fb,'ant4X)dar): 

y  ■  8C|wave(fb,(&jnition); 

y  •  8qwave(fb.duration.'antipodtf ): 

y  •  8qwave(fb,f8,nbrcydas); 

y  >  sqwave(fb,fe,nbrcycles,'antipodai’); 

Description 

8awwave(fb)  -  Generates  one  second  of  a  baseband,  sawtooth  wave  at  cycle  fre¬ 
quency,^,  and  sampling  frequency  of  8192  Hz. 

sawwave(fb,duralion)  -  Generates  duration  seconds  of  a  baseband,  sawtooth  wave  at 
cycle  frequency,^,  and  sampling  frequency  of  8192  Hz. 

sqwave(fb.duralion.f8)  -  Generates  duration  of  a  baseband  unipolar  [0,1]  square 
wave  at  cycle  frequency,  yb,  and  sampling  frequency, /5. 

The  argument  ‘antipodal’  changes  the  amplitude’  range  of  the  ouq>ut  waveftMrm 
from  [0,1]  to  [-1,1].  Note:  If  the  (^le  frequency,^,  and  the  sampling  frequency, 
fs,  are  both  equal  to  one,  the  resulting  vector  will  contain  alternating  values  of  [0,1] 
or  [-1.1]. 

Examples 

Generate  30  cycles  of  a  square  wave  with  a  cycle  frequency  of  25  Hz  sampled  at 
8000  Hz,  varying  between  0  and  5  volts. 

X  «  5  *  sqwave(25.30/25.8000); 


I  «««•  0*01 


Generate  30  bits  of  a  50  bit-per-second  signal  sampled  at  8000  Hz.  Use  this  base¬ 
band  signal  as  the  modulating  signal  input  to  dsbac  to  generate  a  phase-shift  keyed 
signal. 

X  ■  8qwave(50/2.2*30/50.8000, 'antipodal'); 
y-dsb6C(x.2000.8000); 
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Generate  a  vector  with  100  elements  alternating  between  1  and  0. 
y-aqwave(1 ,1,100) 

y- 

1 

0 

1 

0 

1 

0 

1 

(and  so  forth) 

Limitations 

The  waveform  is  created  by  first  creating  a  single  cycle  of  the  waveform  and  then 
repeating  this  cycle  the  required  number  of  times.  Hence,  when  the  cycle  rate  is  not 
evenly  divisible  into  the  sampling  frequency,  tnmcation  errors  occur.  The  length 
of  the  output  vector,  y,  can  be  computed  by  die  formula: 

length  >  flooriduratioH  floor 

The  exact  cycle  frequency  can  be  computed  from: 


flot 

See  Also 

antpodal,  unipolar,  tiiwave,  sawwave 
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STR2MASC,MASC2STR 


Purpose 

Convert  a  string  into  its  binary  ASCII  representatitm  as  generated  by  a  modem. 

Synopsis 

y  •  8tr2ma8C(‘strino‘) 

y  -  str2maso(‘sirina’.databits,’perity’.stopbits) 
y  -  ma8c28lr(x) 

y  -  masc2str(x,datab(i8,'parity’,stopbits): 

Description 

str2masc(‘strinjj‘)  -  Converts  string  to  a  vector  of  1  ’s  and  O’s  rq>resenting  raw,  eight- 
bit  ASCn  coding. 

8lr2masc(‘string',dalabit8.’parity’.6topbits)  •  Converts  string  to  a  vecUM*  of  I’s  and  O’s 
representing  modem  encoded  ASCII.  The  number  of  databits  can  be  7  or  8.  Parity 
can  be  *n’  for  none,  ’o’  for  odd,  <x  ’e’  for  even.  The  number  of  stopbits  can  be  1  or 
2.  One  start  bit  is  always  used.  Valid  combinations  are; 

7nl,  7el,  7ol,  7n2, 7e2, 7o2, 8nl,  8n2 

The  output  argument,  y,  is  nN  by  1  vector  where  N  is  based  upon  the  length  of 
string  and  the  combination  of  the  parameters  bits,  parity,  and  stopbits. 

masc2str(x)  •  Converts  a  vector  of  1  ’s  and  O’s  representing  raw  eight-bit  ASCII  cod¬ 
ing  into  a  string. 

masc2slr(x, databits, 'parit/, stopbits)  -  Converts  a  vector  of  I’s  and  O’s  rq>resenting 
the  binary  modem  ASCII  into  a  string.  Valid  argument  combinations  are  the  same 
as  those  listed  above  for  str2masc.  Parity  and  start/stq}  bit  errm^  are  detected  and 
shown  in  the  ou^rut  string  by  the  following  special  characters; 

startbit  error,  ~,  first  stt^bit  error 
’,  second  stopbit  error,  +,  parity  error 

Parity  errcns  take  precedence  over  start/stop  bit  errms  in  the  ou^ut  string. 

Examples 

Generate  the  binary  rqrresentation  of  “BZ”  sent  at  7  bits-per-character  with  even 
parity  and  2  stop  bits.  Convert  back  to  a  string. 

y  >  sli2ma8c(’BZ’,7,’e’.2); 

y- 

{110000 1011111011010111] 

X  •  masc28lr(y.7,’e’,2); 

Xm 

BZ 
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TRIWAVE 


Purpose 

Generate  •  baseband,  triangular  wave. 

Synopsis 

y  ■  lriwave(lb); 

y  -  triwave^'anlipodar): 

y  ■  triwave(fb,clurtfion): 
y  ■  tnw«ve^.duralion.’antipodar); 
y  >  triwav6(lb, duration, fa): 
y  ■  triwaye(fb,duration.fs.’antipodar): 

Description 

tiiwaye(fb)  -  Generates  one  second  of  a  baseband,  triangular  wave  at  cycle  fre> 
quency.yb,  and  sampling  frequency  of  8192  Hz. 

triwa/e(lb,duralion)  •  Generates  duration  seconds  of  a  baseband,  triangular  wave  at 
cycle  frequency./^,  and  sampling  frequency  of  8192  Hz. 

tiiwave(fb, duration, fa)  -  Generates  nbrcycles  of  a  baseband,  triangular  wave  at  cycle 
frequency, and  sampling  frequency, /s. 

The  argument  ‘antipodal’  changes  the  amplitude  range  of  the  ouq>ut  waveform 
from  [0,1]  to  [-1,1]. 

Example 

Generate  30  cycles  of  a  triangular  wave  with  a  cycle  frequency  of  25  Hz  sampled 
at  8000  Hz,  varying  between  0  and  5  volts. 

X  .  5  *  tiwave(2S,3Q/25,8000) 


Limitations 

The  wavefonn  is  created  by  first  creating  a  single  cycle  cd  the  waveform  and  dien 
rq)eating  this  cycle  the  required  number  of  times.  Hence,  when  the  cycle  rate  is  not 
evenly  divisible  into  the  sampling  frequency,  truncation  errors  occur.  The  lengfh 
of  the  ouq>ut  vector,  y,  can  be  computed  by  Ae  formula: 
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Itiigik  m  flootidurtttloH  -floor 


The  exact  cycle  frequency  can  be  conputed  from: 


A  "mathematically  perfect”  triangular  wave  is  generated  such  that  successive 
cycles  fit  perfectly.  Mthout  getting  into  a  long  discussion,  tte  bottom  line  is  dut 
when  the  number  of  samples-per-cycle  is  odd,  the  maximum  value  of  y  will  not  be 
exactly  1,  as  occurs  when  the  number  of  samples-per-cycle  is  even.  In  addition,  die 
last  cycle  does  not  return  exactly  to  zero  in  either  case  (die  beginning  of  die  next 
cycle  would  have  started  at  zero).  The  "floor”  term  in  the  above  equation  is  the 
number  of  samples-per-cycle. 

See  Also 

aqwave.  sawwave 
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UNIPOLAR 


Purpose 

Generate  •  btsebend,  unipdar  [0,1]  signtl. 

Synopsis 

y  ■  unipoler(fb,mig 
y  -  unipoiar(lb.mig,f») 

Description 

un4)ol«<flMn80)  -  Generates  e  baseband  unipolar  signal  at  bit  rate,^.  Sampling  fre¬ 
quency  defaults  to  8192  Hz.  If  msg  is  a  scalar,  a  message  will  be  generated  as  a 
randmn  binary  sequence  of  length  msg  with  Pr(0)=Pr(l)sO.S.  If  msg  is  a  vector,  it 
must  be  a  vector  of  0*s  and  I’s.  The  length  of  the  ouq>ut  vector,  y,  can  be  computed 
by  the  formula: 

length  a  nbrbaud  -floor 
unipoiar(fb.m8a,fe)  -  Sets  the  sampling  frequency  to/s. 

Note:  If  the  bit  rate,)Z>,  and  the  sampling  frequency,  /s,  are  both  equal  to  one,  the 
resulting  vector  will  cmitam  alternating  values  of  [0,1]. 

Example 

Generate  30  bits  of  a  SO  bit-per-second  signal  sampled  at  8000  Hz.  Use  this  base¬ 
band  signal  as  the  modulating  signal  input  to  dsbec  to  generate  an  on-off  keyed  sig¬ 
nal. 

x-unipolar(50.30.8000); 

y-d8b6C(x,2000.8000); 

Limitations 

•  The  {M'obability  1  occurs  is  not  exact  for  a  small  number  of  bits. 

•  ff  die  sampling  frequency  is  not  an  exact  multiple  of  the  bit  rate,  uneiqiected 
results  may  occur.  For  example:  unipolar(9,18,20)  induces  an  ouq>ut  of  36  samples 
vice  die  expected  40  sample  (i.e.  18  /  9  »  2  sectmds  at  20  samples-per-seco^) 
This  occurs  since  floor(fiifl))  =  2  samples-per-bit. 

See  Also 

anlpodal.  aqurave 
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VECTARMA 


The  Vector  ARMA  Dodad  provides  an  interactive,  graphical  environment  within 
Matlab  to  model  Auto-Regressive  (AR)  and  Auto-Regressive/Moving  Average 
(ARMA)  processes.  With  the  Vector  ARMA  Dodad,  you  can: 

•  Model  AR  fmxesses  using  the  autocorrelation,  covariance,  modified 
covariance  or  Burg  methods. 

•  Model  ARMA  process  using  the  Prony,  Durbin  or  Shank  methods. 

•  Generate  minimum-phase  ARMA  models. 

•  Drive  the  model  with  an  impulse  function  or  Gaussian  white  noise. 

•  Plot  the  impulse  response. 

•  Plot  the  error  signal. 

•  Generate  a  chain  of  models.  This  feature  is  useful  in  modeling  parte  of 
speech. 

•  Play  the  model  and  errcx'  signal  on  the  workstation’  speaker. 


Starting  the  Vector  Filter  Dodad 

The  Vector  ARMA  Dodad  can  be  started  in  one  of  two  different  ways.  The  first 
one  is  to  start  the  dodad  with  a  vector  name  as  an  argument.  A  vecux*  is  defined  as 
a  IxN  or  a  Nxl  Matlab  variable  and  it  must  exist  in  die  Matlab  work^ace  before 
starting  the  Vecuv  ARMA  Dodad.  To  model  a  vector  named  ’’myvector”.  the  com¬ 
mand  *\eGUma(myvoGtor)”  needs  to  be  executed.  Supplyiitg  the  command  with  a 
secoiul  argument  of  "1”  starts  the  dodad  with  a  sampling  firequency  of  one. 

No  output  arguments  are  supported.  The  model  must  be  saved  using  the  save  push¬ 
button  as  described  under  common  controls. 

The  second  way  to  start  die  dodad  is  to  execute  the  ‘Vectarma”  command  without 
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any  arguments.  This  option  can  be  chosen  only  after  the  vet^edh  or  vedfil  dodad 
has  been  used  to  mod^  a  vector.  Every  time  a  vector  is  modified  using  these 
dodads,  the  modified  vector  is  saved  to  a  cmnmon  vector  in  the  global  area  of  the 
Matlab  workspace.  When  no  vecto'  argument  is  supplied,  this  common  vector  is 
loaded  into  the  vectarma  dodad. 

Vector  ARMA  Dodad  Specific  Controls 

In  addition  to  controls  discussed  under  common  controls,  the  following  controls  are 
specific  to  die  Vector  ARMA  Dodad  (vectamui).  For  the  most  part,  these  control 
simply  set  the  filter  parameters.  No  filtering  actually  takes  place  until  the  Apply 
pushbutton  is  selected.  Additionally,  not  all  of  the  following  controls  are  available 
at  any  one  time.  Only  controls  pertinent  to  the  currendy  selected  model  type  are 
displayed. 

Model  lype  Popupmenu 

Purpose 

Selects  between  AR  and  ARMA  modeling. 

Use 

1.  Place  the  cursor  on  the  Model  lypc  popupmenu  and  click  the  left  mouse  button 
once  to  q)en. 

2.  Place  die  cursor  on  the  desired  model  type  (AR  or  ARMA)  and  click  the  left 
mouse  button. 

Result 

The  desired  model  type  is  selected  and  displayed  on  the  pt^upmenu.  Additional 
controls  are  added  tx*  removed  dependent  upon  the  model  type  chosen. 

Model  Method  Radiobuttons 

Purpose 

Selects  the  modeling  method  to  be  implemented. 

Use 

1.  Place  the  curstx  on  the  radiobutton  of  the  desired  modeling  method  and  dick  the 
left  mouse  button  once  to  select. 

Result 

The  radiobutton  of  the  desired  modeling  method  is  selected  and  all  others  are 
turned  off. 

Min  Phase  Checkbox 
Purpose 

Tbms  minimum  phase  ARMA  modeling  on  or  off. 

Use 

1.  Place  the  curscx  on  the  checkbox  and  click  the  left  mouse  button  until  the  check¬ 
box  is  turned  on  or  off  as  desired. 

Result 

The  ARMA  model  transfer  fimction  polynomials  are  forced  to  a  minimum  phase 
condition  when  generated. 
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Notes 

•  This  centred  is  enly  available  when  the  current  modeling  m^iod  it  set  to  ARMA. 

•  A  minimum-phase  model  may  or  may  not  best  fit  the  moekled  data. 

P  Order  Pbpupmenu 

Purpose 

Sets  the  order  of  die  AR  model  (poles). 

Use 

1.  Place  the  cursor  on  the  wder  popupmenu  and  click  the  left  mouse  butten  once  to 
open  the  menu. 

2.  Place  the  cursor  m  the  desired  ordtf  (1.2,3.4,6.8,10.12,14)  and  dick  die  left 
mouse  button.  If  the  desired  order  is  not  one  of  the  default  opdmis,  it  can  be 
entered  manually  by  selecting  the  ‘"User”  option.  In  this  case,  an  edit  box  ^^lears 
in  place  of  the  popupmenu.  Enter  die  desired  order  into  the  edit  box  and  {»ess 
return. 

Result 

The  AR  order  is  set. 

Q  Order  Popupmenu 

Purpose 

Sets  the  order  of  the  MA  model  (zeros). 

Use 

1.  Place  the  cursor  on  the  order  popupmenu  and  click  the  left  mouse  button  once  to 
open  the  menu. 

2.  Place  the  cursor  on  the  desired  order  (1,2,3,4,6,8,10,12,14)  and  click  the  left 
mouse  button.  If  the  desired  wder  is  not  one  of  the  default  options,  it  can  be 
entered  manually  by  selecting  the  ‘TJso'*’  option.  In  this  case,  an  edit  box  ^ipears 
in  place  of  the  popupmenu.  Enter  the  desired  order  into  the  edit  box  and  press 
return. 

Result 

The  MA  (xder  is  set. 

Notes 

This  control  is  only  available  when  die  current  modeling  method  is  set  to  ARMA. 

Plots  Checkboxes 
Purpose 

Selects  the  plots  to  be  displayed. 

Use 

1.  Place  the  cutsot  over  the  desired  plot  and  dick  the  left  mouse  button  to  turn  the 
plot  on  or  off. 

Result 

The  desired  plots  are  displayed  when  the  model  is  generated. 

Driving  Source  Radiobuttons 

Purpose 

Chooses  between  an  impulse  and  Gaussian  white  noise  to  drive  die  model. 
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Use 

1.  Place  tl»  cursor  over  the  desired  driving  source  and  dick  the  Idt  mouse  button 
to  select. 

Result 

The  desired  source  is  used  when  the  model  is  applied. 

Model  Mark  Start/Mark  End  Pushbuttons 
Purpose 

Marks  the  beginning  and  end  of  aportion  of  the  vector  to  be  used  in  generating  the 
model. 

Use 

1.  Select  desired  mark  pushbutum  (Begin  or  End). 

2.  After  selection,  the  cursor  changes  to  a  crosshair  and  the  plot  title  changes  to 
either  “Mark  start  of  model  with  cursor”  or  “Mark  end  of  model  with  cursor”  as  ai^opri- 
ate. 

3.  Move  the  cursor  to  the  desired  mark  location  and  click  the  left  mouse  butkm. 
When  placing  the  mark,  only  the  location  on  the  time  axis  is  of  importance.  Plac* 
ing  a  mark  past  either  end  of  the  plot  sets  the  mark  to  the  first  or  last  value  dis¬ 
played  on  the  horizontal  axis. 

Result 

After  placing  a  mark,  the  cursor  changes  back  into  an  arrow  and  the  plot  title  is 
cleared.  The  model-begin  mark  is  displayed  as  a  vertical  dashed  line  and  the 
model-end  mark  is  displayed  as  a  vertical  dotted  line. 

Notes 

The  default  option  uses  the  entire  vector  to  generate  the  model. 

Chain  Period  Mark  Start/Mark  End  Pushbuttons 
Purpose 

Marks  the  begiiining  and  end  of  the  Chain  period. 

Use 

1.  Select  desired  mark  pushbutton  (Start  or  End). 

2.  After  selection,  the  curstx'  changes  to  a  crosshair  and  the  plot  title  is  changed  to 
either  “Mark  start  of  period  with  cursor”  or  “Mark  end  of  period  with  cursor”  as  ^jprtqjri- 
ate. 

3.  Move  the  cursor  to  the  desired  mark  location  and  click  the  left  mouse  button. 
When  placing  the  mark,  only  the  location  on  the  time  axis  is  of  impcatance.  Plac¬ 
ing  a  mark  past  either  end  of  the  plot  sets  the  mark  to  the  first  or  last  value  dis¬ 
played  on  the  hmizontal  axis. 

Result 

After  pladng  a  mark,  the  cursm*  changes  back  into  an  arrow  and  the  plot  title  is 
cleared.  The  period-begin  mark  is  displayed  as  a  vertical  dashed  liM  and  the 
period-end  mark  is  displayed  as  a  vertiad  dotted  line. 

Notes 

•  The  default  option  uses  the  entire  vector  length  as  the  period. 

•  The  period-end  cannot  extend  past  the  length  of  the  vector  from  the  period  begin 
mark. 
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•  The  period  b^in  maik  is  set  coincidental  with  the  model  b^in  mark  by  the 
Model  Begin  pushbutton. 

Chain  Checkbox 

Purpose 

Tiims  the  gmieraticm  of  a  chain  of  models  on  or  off. 

Use 

1.  Place  the  cursor  on  the  Qiain  checkbox  and  click  tlm  left  mouse  button  to  turn 
the  chain  on  or  off. 

Result 

The  model  is  chained  the  number  of  times  specified  by  the  Chain  Length  popup- 
menu  when  the  Chain  checkbox  is  checked. 

Chain  Length  Popupmmiu 

Purpose 

Sets  the  number  of  times  the  model  is  rq)eated  in  the  chain. 

Use 

1.  Place  the  cursOT  on  the  Chain  Length  pq}upmenu  and  dick  die  left  mouse  but¬ 
ton  once  to  open  the  menu. 

2.  Place  the  cursor  on  the  desired  chain  length  (10,20,30,40,50  or  60)  and  click  the 
left  mouse  button.  If  the  desired  chain  length  is  not  one  of  the  default  optitms,  the 
chain  length  can  be  entered  manually  by  selecting  the  “User”  cation.  In  this  case, 
an  edit  box  i^ipears  in  place  of  the  popupmenu.  Enter  the  desired  chain  length  into 
the  edit  box  and  press  return. 

Result 

The  desired  chain  length  is  set. 

Notes 

This  feature  is  useful  for  modeling  parts  of  speech. 

Chain  Period  Edit  Box 

Purpose 

•  Displays  the  time  between  the  Period  Markers. 

•  Allows  the  operator  to  manually  enter  the  chain  period. 

Use 

1.  Place  the  cursor  over  the  edit  box  and  click  the  left  mouse  button. 

2.  &iter  the  desired  chain  period  length  and  press  return. 

Result 

The  chain  period  is  displayed  in  the  edit  box. 

Notes 

The  period  end  cannot  extend  past  the  length  of  the  vector  from  the  period-bogin 
mark.  When  setting  the  period  past  the  length  of  the  vectcx,  the  period  is  adjusted 
so  that  the  period-end  mark  is  set  to  the  end  of  the  vector. 

Play  Desired,  Play  Model,  Play  Error  Popupmenu  Menu  Items 

Purpose 

Sends  the  desired  data  vectw,  the  model  vector  or  the  erro*  vector  to  the  worksta- 
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ti(»*s  audio  ou^ut. 

Use 

1.  Place  the  cursor  on  the  Play  p<^upmenu  and  click  the  left  mouse  buttcm  to  open 
the  menu. 

2.  Select  the  Play  All  menu  item. 

Result 

The  data  is  sent  to  the  audio  ou^ut. 

Note 

This  option  is  ideal  for  testing  speech  models. 

Save  Pushbutton 
Purpose 

Saves  die  modeled  data,  model  and  error  signals  to  the  Matlab  workspace  under  a 
new  variable  name,  if  a  model  was  generated.  In  addition,  the  model  parameters 
are  saved  to  the  given  variable  named  suffixed  with  an  underscore  and  an  “a”  for 
the  AR  parameters  and  a  “b”  for  the  MA  parameters. 

Use 

1.  Place  the  cursor  on  the  Save  Pushbutton  and  click  the  left  mouse  button. 

2.  The  Save  and  Load  pushbuttons  are  replaced  by  an  edit  box.  If  save  has  already 
been  used,  the  previously  entered  vector  name  is  contained  in  the  edit  box. 

3.  Place  the  cursor  on  the  edit  box  and  click  the  left  mouse  button.  Using  the  key¬ 
board,  backspace  over  any  previous  entry  and  then  type  in  the  desired  vector  name 
(following  standard  Matlab  variable  name  conventions). 

4.  Press  the  return  key. 

Result 

The  data  is  saved  into  an  Nx4  matrix  where  the  first  column  is  the  time  indices,  the 
second  column  is  the  modelled  data,  the  third  column  is  the  model  and  the  fourth 
column  is  the  error  signal. 
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VECTEDIT _ 

The  Vector  Edit  Dodad  provides  an  interactive,  graphical  environment  within  Mat- 
lab  to  edit  vectors.  This  is  ideally  stiited  fw  editing  vectors  that  represent  digitized 
signals  such  as  speech.  V^th  the  Vector  Edit  Dodad,  you  can: 

•  Cut  unwanted  sections  of  a  signal. 

•  Crop  unwanted  sections  of  a  signal. 

•  Rearrange  sections  of  a  signal  by  cut  and  paste. 

•  Adjust  the  amplitude  (volume)  of  all  or  a  portion  of  a  signal. 

•  Introduce  periods  of  silence  into  a  signal. 

•  Simply  zoom  in  on  portions  of  a  signal  for  a  better  look. 


Starting  the  Vector  Edit  Dodad 

The  Vector  Edit  Dodad  can  be  started  in  one  of  two  ways.  The  first  one  is  to  start 
the  dodad  with  a  vector  name  as  an  aigument.  A  vector  is  defined  as  a  IxN  (m*  a 
Nxl  Madab  variable  and  it  must  exist  in  the  Matlab  workspace  before  starting  the 
vector  edit  dodad.  To  edit  a  vector  named  “myvector”,  the  command  ‘VectedH(- 
myvector)”  needs  to  be  executed. 

No  output  arguments  are  supported.  Edited  vectors  must  be  saved  using  the  save 
pushbutton  as  described  under  common  controls. 

The  second  way  to  start  the  Vector  Edit  Dodad  is  to  execute  the  ‘Vectecfit”  com¬ 
mand  without  any  arguments.  This  can  be  done  only  after  the  vectedit  or  vedfit 
dodad  has  been  used  to  modify  a  vector.  Everytime  a  vector  is  modified  using  these 
dodads,  the  modified  vector  is  saved  to  a  common  vecttx'  in  the  global  area  of  the 
Matlab  workspace.  When  no  vector  argument  is  supplied,  this  common  vector  is 
loaded  into  the  Vector  Edit  Dodad. 
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Vector  Edit  Dodad  Specific  Controls 

In  addition  to  ctmtrols  discussed  in  the  Common  Controls  chiq>ter,  the  following 
controls  are  specific  to  the  vectcr  edit  dodad  (vectedh). 

Cut  Pushbutton 
Purpose 

Removes  a  marked  section  of  the  vectm  placing  the  cut  section  into  the  cut  and 
paste  buffer. 

Use 

1.  Mark  the  section  to  be  cut  using  the  Mark  Start  and  Maiic  &id  pushbuttons. 

2.  Place  the  cursor  on  the  Cut  pushbutton  and  click  the  left  mouse  button. 

Result 

The  marked  section  is  cut  and  the  display  is  redrawn. 

Notes 

•  The  removed  portion  can  be  restored  by  immediately  selecting  the  Cut  pushbut¬ 
ton  a  second  time.  This  option  is  void  if  either  of  the  marks  are  moved  before 
selecting  the  Cut  pushbutton  for  the  second  time  or  any  other  editing  options  is 
performed. 

•  If  the  beginning  and  end  marks  are  located  at  the  beginning  and  end  of  the  vector 
(i.e.  the  entire  vector  is  marked),  selecting  the  Cut  pushbutton  has  no  effect. 

•  The  common  vector  is  modified. 

Copy  Pushbutton 

Purpose 

Copies  a  marked  section  of  the  vector,  placing  the  copied  portion  into  the  cut  and 
paste  buffer. 

Use 

1.  Mark  the  section  to  be  copied  using  the  Mark  Start  and  Mark  End  pushbuttons. 

2.  Place  the  cursor  on  the  Cc^y  pushbutton  and  click  the  left  mouse  button. 

Result 

The  copied  section  is  placed  into  the  cut  and  paste  buffer. 

Paste  Pushbutton 

Purpose 

Inserts  the  contents  of  the  cut  and  paste  buffer  into  the  vector  at  a  point  selected 
with  the  mouse. 

Use 

1.  Place  the  cursor  on  the  Paste  pushbutton  and  click  the  left  mouse  button. 

2.  The  cursOT  changes  into  a  crosshair  and  the  display  title  changes  to  "Mark  paste 
insertion  point  with  cursor.” 

3.  Move  the  cursor  to  the  desired  insertion  point  and  elide  the  left  mouse  button. 

Result 

The  contents  of  the  cut  and  paste  buffer  are  inserted  into  the  vector  and  the  vector 
display  is  redrawn. 

Notes 

•  If  the  cut  and  paste  buffer  is  empty,  no  action  is  performed  (a  cut  or  copy  has  to 
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be  performed  first  to  fill  the  cut  aiid  paste  buffer). 

•  The  contents  of  the  cut  and  paste  buffer  are  not  cleared.  Multiple  copies  of  die  cut 
and  paste  buffer  can  be  insert  by  reselecting  the  Paste  Pushbutton. 

•  The  cmnmon  vector  is  modified. 

Crop  Pushbutton 

Purpose 

Deletes  sections  of  the  vector  tnitside  the  currently  mark  region. 

Use 

1.  Mark  the  region  to  be  retained  using  the  Mark  Start  and  Mark  End  pushbuttons. 

2.  Place  the  cursor  on  the  Crop  pushbutton  and  click  the  left  mouse  button. 

Result 

The  sections  of  the  vector  outside  the  marked  region  are  deleted. 

Notes 

•  The  deleted  sections  are  not  saved. 

•  The  cut  and  paste  buffer  is  not  modified. 

•  The  common  vecttM*  is  modified. 

Mean  User  Edit  Box 

Purpose 

•  Display  the  ciurent  vector’s  mean  value. 

•  Remove  the  mean. 

•  Change  the  mean  (if  the  vecttur  rq;>resents  a  voltage,  this  is  equivalent  to  chang¬ 
ing  the  DC  offset). 

Use 

Remove  the  mean: 

1.  Select  the  edit  box  by  clicking  the  left  mouse  button  while  the  cursw  is  located 
over  the  edit  box. 

2.  Enter  zero  into  the  edit  box  and  press  return. 

Change  the  mean: 

1.  Select  the  edit  box  by  clicking  the  left  mouse  button  while  the  cursor  is  located 
over  the  edit  box. 

2.  Enter  the  desired  offset  value  into  the  edit  box  and  press  return. 

Result 

The  mean  is  adjusted  and  the  display  is  redrawn. 

Notes 

•  After  the  mean  is  changed,  the  mean  is  recomputed  and  redisplayed  in  the  edit 
box.  Note  that  after  removing  the  signal  mean,  the  resulting  value  may  not  be 
exactly  zero  due  to  round-off  errcxs. 

•  The  common  vecUx*  is  modified. 

Zero  Marked  Pushbutton 

Purpose 

Zeroes  out  the  vecUx  in  the  marked  section.  This  option  is  useful  fn-  creating  abso¬ 
lute  periods  of  silence. 

Use 

1.  Mark  the  section  to  be  zeroed  out  using  the  Mark  Start  and  Mark  End  pushbut- 
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tons. 

2.  nace  the  cursor  on  the  Zero  Marked  pushbutton  and  dick  the  left  mouse  button. 

Result 

The  marked  vector  section  is  set  to  zero  and  the  display  is  redrawn. 

Notes 

The  common  vector  is  modified. 

Volume  Full/Marked  Popupmenu  and  Volume  Slider 

Purpose 

Adjust  the  amplitude  (volume)  of  the  vectm-  or  a  marked  section  0%  to  200%  of 
it’s  current  amplitude. 

Use 

1.  If  desired,  mark  a  section  to  be  adjusted  using  the  Mark  Start  and  Mark  End 
pushbutums. 

2.  Adjust  the  volume  slider  to  the  d^ed  amplitude  change  percent. 

3.  Place  the  cursor  on  the  Volume  popupmenu  and  click  the  left  mouse  button  once 
to  <^n  the  pt^upmenu. 

4.  Place  the  cursor  on  the  desired  function  (Volume  Full  or  Volume  Marked)  and 
click  the  left  mouse  button. 

Result 

The  amplitude  is  adjusted  as  requested  and  the  display  is  redrawn. 

Notes 

•  If  the  vector  is  to  be  saved  in  an  audio  file  fcamat,  care  must  be  taken  to  ensure 
the  new  range  of  values  do  not  exceed  the  limits  of  the  target  format.  For  example, 
Microsoft  Window  “wave”  files  and  Soundblaster  “voice”  files  usually  contain 
only  8-bit  values  (-128  to  +127).  Values  outside  this  range  are  clipped  when  saved 
to  Ae  audio  file.  The  anqrlitude  can  be  properly  adjusted  to  range  between  +127 
and  -127  by  the  command  “» tnyvector  >  127  *  myvector  /  max(myvectoi);”. 

•  The  common  vector  is  modified. 
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VECTFILT _ 

The  Vector  Hlter  Dodad  provides  an  interactive,  gn4>hical  environment  within 
Matlab  to  filter  vectors.  Tliis  is  ideally  suited  for  filtering  vectors  that  refnesent 
digitized  signals  such  as  speech  or  digital  bandpass  signals.  With  the  Vector  Filter 
Dodad,  you  can: 

•  PerfcMm  lowpass.  highpass.  bandpass  and  stopband  filtering. 

•  Filter  using  a  Finite  Impulse  Response  filter. 

•  Filter  using  Butterwmth,  Chebychev  or  Elliptical  Infinite  Impulse 
Response  (IIR)  filters. 

•  Vary  filter  parameters  such  as  FIR  window,  filter  order  and  stopband 
attenuation. 

•  View  the  filter  transfer  function  before  applying. 

•  Inspect  the  signal  by  zooming  in  on  the  time  domain  display. 

•  Inspect  the  signal  by  zooming  in  on  the  power  spectral  density  display. 

•  Determine  the  audible  effects  of  filtering  by  playing  filtered  signals  on  the 
speaker. 


Starting  the  Vector  Filter  Dodad 

The  VectOT  Filter  Dodad  can  be  started  in  one  of  two  ways.  The  first  one  is  to  start 
the  dodad  with  a  vector  name  as  an  aigument.  A  vector  is  defined  as  a  IxN  or  a 
Nxl  Matlab  variable  and  it  must  exist  in  the  Matlab  workspace  befcxe  starting  the 
Vectw  Filter  Dodad.  To  filter  a  vector  named  "myvedor”,  the  command  *Vectfit( 
myvector)”  needs  to  be  executed. 

No  ou^ut  arguments  are  supptated.  Filtered  vectors  must  be  saved  using  the  Save 
pushbutton  as  described  imder  common  oorrtrols. 


SPC  Toolbox  Urn’s  Guide 


102 


VECIPILT 


The  second  way  to  start  the  filter  dodad  is  to  execute  the  “vedfit”  command  with¬ 
out  any  afguments.  This  can  be  done  tmly  after  the  vededit  or  vedtime  dodad  has 
been  used  to  modify  a  vector.  Every  time  a  vector  is  modified  using  these  dodads, 
the  modified  vector  is  saved  to  a  common  vector  in  the  global  area  of  the  Matlab 
workspace.  When  no  vector  argument  is  supplied,  this  common  vector  is  loaded 
into  the  >^tor  Filter  Dodad. 

Vector  Filter  Dodad  Specific  Controls 

In  addition  to  controls  discussed  under  common  controls,  the  following  controls 
are  specific  to  the  Vector  Filter  Dodad  (vedfit).  Fm*  die  most  part,  these  ccmtrols 
simply  set  the  filter  parameters.  No  filtering  actually  takes  place  until  the  ^>ply 
pushbutton  is  selected.  Additionally,  not  all  of  the  following  controls  are  available 
at  any  one  time.  Only  controls  pertinent  to  the  currently  selected  filter  and  filter 
type  are  displayed.  Fix'  example,  there  is  no  need  to  set  an  upper  cutoff  frequency 
when  a  low  or  highpass  filter  type  is  selected. 

Filter  Popupmenu 

Purpose 

Selects  the  filtering  method  to  be  used. 

Use 

1.  Place  the  curscx  on  the  Filter  pqsupmenu  and  click  the  left  mouse  button  once  to 
open. 

2.  Place  the  cursor  on  the  desired  filter  (FIR,  Butterworth,  Chebychev  fype  1,  Che- 
bychev  2,  or  Elliptical)  and  click  the  left  mouse  button. 

Result 

The  desired  filtering  method  is  selected  and  displayed  on  the  popupmenu.  Addi¬ 
tional  controls  are  added  or  removed  dependent  upon  the  filtering  method  chosen. 

Notes 

The  filter  are  implemented  using  functions  from  the  Matlab  Signal  Processing 
Toolkit  which  is  required  to  use  this  dodad. 

Filter  l^pe  Popupmenu 

Purpose 

Selects  the  filter  type  to  be  implemented. 

Use 

1.  Place  the  cursor  on  the  Filter  Type  popupmenu  and  click  the  left  mouse  button 
once  to  (^>en. 

2.  Place  the  curstx  on  the  desired  filter  type  (Lowpass,  Higlq>ass,  Bandpass,  <x 
Stopband)  and  click  the  left  mouse  button. 

Result 

The  desired  filtering  method  is  selected  and  displayed  oa  the  popiqmienu.  Addi¬ 
tional  controls  are  added  or  removed  dq)endent  upon  tte  filter  type  chosen. 

Window  Popupmenu 

Purpose 

Selects  die  type  of  data  prefiltering  window  to  be  applied  when  applying  a  FIR  fil- 
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tenag. 

1.  Place  die  cunor  on  die  window  p<^upmenu  and  click  dM  left  mouse  buttm  mice 
to  open. 

2.  Place  the  ctffsor  on  the  desired  window  (Rectangular.  Triangular.  Hamming. 
Hanning.  Blackman,  or  Bartlett)  and  click  the  left  mouse  button. 

The  desired  window  is  selected  and  displayed  on  the  popupmenu. 

•  The  windows  are  implemented  using  functions  from  the  Matlab  Signal  Process¬ 
ing  Toolkit  which  is  required  to  use  this  dodad. 

•  This  ccmtrol  is  only  available  when  the  current  filtering  method  is  set  to  FIR. 


Order  Popupmenu 

Purpose 

Sets  the  filter  order. 


Result 


Notes 


1.  Place  the  cursor  on  the  Order  pt^upmenu  and  click  the  left  mouse  button  once 
to  open  the  menu. 

2.  Place  the  cursor  on  die  desired  filter  order  and  click  the  left  mouse  buttm.  The 
popupmenu  ccmtains  orders  of  11,  21,  31,  41,  51  and  61  when  the  FIR  filter  is 
selected  and  orders  of  2,  4,  6,  8,  10  and  12  when  an  IIR  filter  is  selected.  If  the 
desired  order  is  not  am<mg  these  options,  the  required  order  can  be  entered  iF/'jii'u- 
ally  by  selecting  the  "User”  opdmi.  In  such  a  case,  an  edit  box  appears  in  place  of 
the  popupmenu.  Enter  the  desired  ordo’  into  the  edit  box  and  fvess  return. 

The  filter  order  is  set. 

•  The  filter  order  for  bandpass  and  stopband  filters  is  twice  the  order  shown  oa  the 
popupmenu.  This  is  consistent  widi  the  input  aigumoit  spedficatimis  for  the  Mat- 
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lab  Signal  Processing  Toolbox  routines. 

•  Higl^>ass  and  stopband  filters  requires  an  even  order.  Odd  orders  are  automati¬ 
cally  rounded  up  one  to  the  next  hi^er  even  order. 

•  The  filters  are  implemented  using  functions  from  the  Matlab  Signal  Processing 
Toolkit  which  is  required  to  use  this  dodad.  Some  restrictions  on  the  filter  order 
apply  when  using  these  functicms.  See  the  signal  processing  toolkit  documentatitm 
for  specifics. 

•  In  general,  FIR  filters  require  a  high  order  to  obtain  a  sharp  rolloff.  IIR  filters  can 
achieve  the  same  or  better  results  using  a  much  lower  order. 

Lower  Cutoff  Frequency  Edit  Box 

Purpose 

•  Displays  the  lower  cutoff  frequency  when  it  is  selected  using  the  mouse. 

•  Allows  the  operator  to  manually  enter  the  lower  cutoff  frequency. 

Use 

When  selecting  by  mouse: 

1.  Display  the  spectrum  of  the  signal  using  the  Spectrum  radiobutton. 

2.  Mark  Ae  lower  cutoff  frequency  using  the  beginning  mark. 

Manual  entry: 

1 .  Place  the  cursor  over  the  edit  box  and  click  the  left  mouse  button. 

2.  Enter  the  desired  lower  cutoff  frequency  and  press  return. 

Result 

The  lower  cutoff  frequency  is  displayed  in  the  edit  box. 

Notes 

•  If  the  lower  cutoff  frequency  is  (1)  less  then  or  equal  to  zero  or  (2)  greater  than  or 
equal  to  fJ2  or  the  upper  cuUdf  frequency,  the  error  message  **lnvaiid  lower  cutoff  fre¬ 
quency  entered.  Try  again...”  is  displayed  in  the  Matlab  command  window  and  die 
edit  box  is  cleared.  This  cycle  continues  until  a  valid  lower  cutoff  frequency  is 
entered. 

•  After  a  cutoff  frequency  has  been  set,  the  transfer  function  of  the  currently 
defined  filter  is  drawn  against  the  background  of  the  signal  spectrum.  An  unstable 
filter  can  generally  be  recognized  from  the  transfer  function. 

•  Setting  the  lower  cutoff  frequency  too  close  to  zero,^/2  or  the  upper  cutoff  fre¬ 
quency  can  lead  to  an  iinstable  filter. 

•  The  effect  of  an  unstable  filter  is  immediately  recognizable  after  iqpplying  the  fil¬ 
ter. 

Upper  Cutoff  Frequency  Edit  Box 
Purpose 

•  Displays  the  upper  cutoff  frequency  when  it  is  selected  using  the  mouse. 

•  Allows  the  operator  to  manually  enter  the  upper  cutoff  frequency. 

Use 

When  selecting  by  mouse: 

1.  Display  the  spectrum  of  die  signal  using  the  Spectrum  radiobutton. 

2.  Maik  toe  upper  cutoff  frequency  using  toe  end  mark. 

Manual  entry: 
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1.  Place  the  cursor  ow  the  edit  box  and  click  the  left  mouse  buttoa. 

2.  Eiiter  the  desired  iqjper  cutoff  frequency  and  press  return. 

Result 

The  ui^)er  cutoff  frequmicy  is  displayed  in  the  edit  box. 

Notes 

•  If  the  upper  cutoff  frequency  is  (1)  less  then  or  equal  to  zero  or  (2)  greater  than  or 
equal  to  ff/2  or  the  lower  cutoff  frequency,  the  error  message  **lnvaid  upper  oulof 
frequency  entered.  Try  again...”  is  displayed  in  the  Matlab  command  window  and  die 
edit  box  is  cleared.  This  cycle  cmitinues  until  a  valid  upper  cutoff  frequency  is 
entered. 

•  Setting  the  uj^ier  cutoff  frequency  too  close  to  zero,  f/2  or  die  lower  cutoff  fre¬ 
quency  can  lead  to  an  unstable  filter. 

•  The  effect  of  an  unstable  filter  is  Immediately  recognizable  after  iq^lying  the  fil¬ 
ter. 

Passband  Ripple  Popupmenu 

Purpose 

Sets  the  amount  of  allowable  ripple  in  the  passband  of  Chebychev  or  elliptical  fil¬ 
ters. 

Use 

1.  Place  the  cursm*  on  the  Passband  Ripple  popupmenu  and  click  the  left  mouse 
button  once  to  qien  the  menu. 

2.  Place  the  cursor  on  the  desired  passband  ripple  and  click  the  left  mouse  button. 
The  popupmenu  contains  ripple  values  of  0.01, 0.05, 0.1, 0.5, 1.0  and  5.0  decibels. 
If  the  desired  ripple  is  not  contained  amoung  the  available  values,  the  ripple  can  be 
entered  manually  by  selecting  the  "User”  option.  In  this  case,  an  edit  Iwx  appears 
in  place  of  the  popupmenu.  Enter  the  desired  ripple  into  the  edit  box  and  press 
return. 

Result 

The  passband  ripple  is  set. 

Notes 

This  control  is  only  available  when  the  current  filtering  method  is  set  to  Cheby¬ 
chev  Type  1  or  E114>tical. 

Stopband  Attenuation 

Purpose 

Sets  the  amount  frequencies  in  the  stopband  are  attenuated. 

Use 

1.  Plme  the  cursor  cm  the  Stopband  Attenuaticm  popupmenu  and  dick  the  left 
mouse  button  cmce  to  open  the  menu. 

2.  Place  the  curscx  cm  the  desired  stopband  attenuaticm  and  click  tiie  left  mouse 
buttcm.  The  popupmenu  contains  attenuation  values  cS  10, 15, 20. 30. 40  and  50  dB 
decibels,  ff  the  desired  attenuaticm  is  not  ccmtained  amount  the  available  values, 
the  attenuaticm  can  be  entered  manually  by  selecting  the  "User”  opticm.  to  this 
case,  an  edit  box  iq>pears  in  place  of  the  pc^upmenu.  Enter  the  desir^  attenuaticm 
into  the  edit  bc«  and  press  return. 
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Result 

The  stopband  attenuatum  is  set. 

Notes 

This  control  is  only  available  when  the  current  filtoing  method  is  set  to  Cheby- 
chev  l^pe  2  or  Elliptical. 

Linear/LogariUimic  RadiotNiUon 
Purpose 

Shows  the  power  ^>ectral  density  display  scaled  linearly,  or  logarithmi¬ 

cally  in  decibels.  20  iot//(/) . 

Use 

1.  Race  the  cursor  cm  the  desired  radiobutton  option  and  click  the  left  mouse  but¬ 
ton. 

Result 

The  display  is  redrawn  apprq)riately. 

Notes 

This  control  is  available  only  when  dre  Spectrum  radiobutton  is  active. 

Apply  Pushbutton 
Purpose 

Applies  the  filter. 

Use 

1.  Click  the  left  mouse  button  while  the  cursor  is  over  the  ^ply  pushbutton. 

Result 

The  filter  specified  by  the  other  controls  is  i^>plied  to  the  vector  and  the  result  is 
displayed. 

Notes 

The  result  of  an  unstable  filter  is  immediately  recognizable.  Selecting  Restme 
returns  the  original  signal. 
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The  Speech  'Ilme-Domain  Analysis  Dodad  provides  and  grt^hical,  interactive 
environment  to  ^>ply  time-domain  methods  of  determining  voice  and  unvoiced 
phonemes  in  speech.  Available  time-domain  analysis  algtxithms  are  the  short-time 
energy  and  magnitude  methods  for  determining  voiced  phonemes  and  the  shmt- 
time  zero  crossings  method  fm*  determining  unvoiced  phonemes.  Voiced  phonemes 
are  characterized  by  a  relatively  high  energy  content  while  unvoiced  phonemes  are 
characterized  by  higher  frequency  content  (hence,  a  higher  number  of  zero  cross¬ 
ings). 

Starting  the  Speech  Time-Domain  Analysis  Dodad 

The  Speech  Tune-Domain  Analysis  Dodad  can  be  started  in  one  of  two  different 
ways.  The  first  is  to  start  the  dodad  with  a  vector  name  as  an  argument.  A  vecttM*  is 
defined  as  a  IxN  or  a  Nxl  Matlab  variable  and  it  must  exist  in  the  Matlab  wOTk- 
space  before  starting  the  Speech  Time-Domain  Analysis  Dodad.  To  analyze  a  vec¬ 
tor  named  “myvector”,  the  command  “vecttimef  myvector)”  needs  to  be  executed. 


0.$  1 
Time  (seconds) 


No  output  arguments  are  supported.  Speech  time-domain  analysis  curves  must  be 
saved  using  the  Save  pushbutton  described  below. 

The  second  way  to  start  the  Speech  Time-Domain  Analysis  Dodad  is  to  execute 
the  “vecttime”  command  without  any  arguments.  This  can  be  done  only  after  the 
vectedit  or  vectfit  dodad  has  been  used  to  modify  a  vector.  Every  time  a  vector  is 
modified  using  these  dodads,  the  modified  vector  is  saved  to  a  common  vector  in 
the  global  area  of  the  Matlab  workspace.  When  no  vector  argument  is  suj^lied, 
this  common  vector  is  loaded  into  the  Speech  Time-Domain  Analysis  Dodad. 
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Speech  Time-Domain  Anaiysis  Dodad  Specific  Controls 

In  addition  to  ccmtrols  discussed  under  common  controls,  the  following  controls 
are  specific  to  the  Speech  Hme-Domain  Analysis  Dodad  (vectlime).  These  ctmtrols 
setup  the  desired  analysis  method,  the  analysis  frame  length  and  overly),  ami  out¬ 
put  ciuve  smoothing  optitms. 

Dme  Analysis  Checkboxes 

Purpose 

Selects  the  analysis  method  to  be  performed. 

Use 

1 .  Place  the  cursor  over  the  desired  aiudysis  and  click  the  left  mouse  button  to  turn 
the  analysis  method  on  or  off. 

Result 

The  desired  is  perfmmed  when  die  ^ply  pushbutton  is  selected. 

Include  Signal  Checkbox 
Purpose 

Includes  a  plot  of  the  time-domain  signal  in  the  ou^ut  graphic  window. 

Use 

1.  Place  the  cursor  over  the  Include  Signal  checkbox  and  click  the  left  mouse  but¬ 
ton  to  turn  the  inclusion  option  on  or  off. 

Result 

The  time-domain  signal  is  included  or  excluded  when  the  ^ply  pushbutton  is 
selected. 

Smoothing  Popupmenu 

Purpose 

Selects  the  smoothing  method  to  be  applied  to  the  output  analysis  curves,  if  any. 

Use 

1.  Place  the  cursor  on  the  Smoothing  popupmenu  and  click  the  left  mouse  button 
once  to  open. 

2.  Place  Ae  cursor  on  the  desired  smoothing  method  (None,  Average  or  Median) 
and  click  the  left  mouse  button. 

Result 

The  desired  smoothing  method  is  selected  and  displayed  on  the  pqjupmenu. 
Smoothing  Filter  Length  Popupmenu 
Purpose 

Selects  the  length  of  the  smoothing  filter  employed  on  the  ou4>ut  analysis  curves. 
Use 

1.  Place  the  cursor  on  the  Length  popupmenu  and  click  the  left  mouse  button  once 
to  open. 

2.  Race  the  cursw  on  the  desired  length  (3,  5,  7,  9,  11,  <x  13)  and  click  the  left 
mouse  button.  If  the  length  is  not  one  of  those  available,  the  desired  value  can  be 
entered  by  selecting  the  ‘'User”  menu  item.  In  this  case,  an  edit  box  appears  in 
place  of  the  popupmenu.  Enter  the  desired  length  into  this  edit  box  and  press 
return. 
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Result 

The  desired  length  is  selected  and  diq>layed  oa  tl^  popufMnenu. 

Time  M^ndow  Frame  Length  Popupmenu 

Purpose 

Selects  the  length  of  the  analysis  window  frame. 

Use 

1 .  Place  cursor  cm  die  Hme  ^^dow  popupmenu  and  click  the  left  mouse  but¬ 
ton  once  to  open. 

2.  Place  the  cursor  on  the  desired  length  (5  ms.  10  ms,  IS  ms,  20  ms,  30  ms,  cr  SO 
ms)  and  click  the  left  mouse  button.  If  the  length  is  not  one  of  those  available,  the 
desired  value  can  be  entered  by  selecting  the  “User”  menu  item.  In  this  case,  an 
edit  box  qipears  in  place  of  the  pr^upmenu.  Enter  the  desired  length  in  millisec¬ 
onds  into  this  edit  box  and  press  return. 

Result 

The  desired  length  is  selected  and  displayed  on  the  popupmenu. 

Time  Window  Overlap  Popupmenu 

Purpose 

Selects  the  length  of  overlap  of  consecutive  analysis  windows. 

Use 

1.  Place  the  cursor  on  the  Window  Overlap  popupmenu  and  click  the  left  mouse 
button  once  to  open. 

2.  Place  the  cursor  on  the  desired  overlap  (10%,  20%,  30%,  50%,  75%  or  90%) 
and  click  the  left  mouse  button.  If  the  desired  overlap  is  not  one  of  the  available 
default  choices,  the  desired  value  can  be  entered  by  selecting  the  “User**  menu 
item.  In  such  a  case,  an  edit  box  appears  in  place  of  the  popupmenu.  Enter  the 
desired  overlap  into  this  edit  box  and  press  return. 

Result 

The  desired  overlap  is  selected  and  displayed  on  the  popupmenu. 

Apply  Pushbutton 

Purpose 

^plies  the  analysis. 

Use 

1.  Click  the  left  mouse  button  while  die  curso*  is  over  the  Apply  pushbutton. 

Result 

The  analysis  method  is  qiplied  as  specified.  A  gr^ihic  window  qiens  to  display 
the  results. 

Save  Pushbutton 
Purpose 

Saves  the  short-time  energy,  magnitude  and  zero-crossing  curves  to  the  Madab 
workspace,  if  the  analysis  was  generated. 

Use 

1.  Place  the  cursor  on  the  Save  Pushbutton  and  click  die  left  mouse  button. 

2.  The  Save  and  Load  pushbuttons  are  rqilaced  by  an  edit  box.  If  save  has  already 
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been  used,  the  previously  entered  vector  name  are  contained  in  the  edit  box. 

3.  Place  the  cursor  on  the  edit  box  and  click  the  left  mouse  button.  Using  the  key¬ 
board,  backspace  over  any  previous  entry  and  then  type  in  the  desired  vector  name 
(following  standard  Matlab  variable  name  conventions). 

4.  Press  the  return  key. 

Result 

The  data  is  saved  into  an  Nx4  matrix  where  the  first  column  is  the  time  indices,  the 
second  column  is  the  sh(»t-time  energy,  the  third  column  is  the  short-time  magni¬ 
tude  and  the  fourth  column  is  the  short-time  zero-crossings. 
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WPERIGRM _ 

Purpose 

Display  the  periodogram  of  a  signal  using  a  linear  scale  (relative  magnitude). 

Synopsis 

vvperignn(x) 

wperigrm(x,fs) 

wperigrm(x,fs.'phase’) 

Description 

wperigrm(x)  -  Computes  the  periodogram  of  x  and  displays  the  result  using  a  linear 
(relative  magnitude)  scale.  Only  die  positive  frequencies  up  to  half  the  sampling 
frequency  are  displayed.  The  default  sampling  frequency  is  8192  Hz. 

wpengrm(x,fs)  •  Sets  the  sampling  frequency  to  fs. 

wperiorm(x.fs, ‘phase’)  -  Splits  the  graphics  window  and  displays  both  the  magnitude 
and  phase  information. 

Example 

Generate  a  frequency-shift  keyed  signal  and  display  it’s  periodogram. 

y-bfsk(512, 512.2048.0.1): 
wperign^y.’phase’);; 


Pow«r  Sp«cti«l  D«n«)ty 


Fraquvney  (Hz) 

Algorithm 

The  periodogram  is  given  by: 
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where  XjO.f)  is  the  fourier  transfran  of  x(t): 

i*T/2 

/-r/2 

Limitations 

The  size  of  the  FFT  used  to  compute  the  fourier  transfmm  is  the  closest  power-of- 
two  greater-than  or  equal-to  the  length  of  x.  Long  input  vectors  require  kmger 
compute  times.  If  jc  is  short  and  a  finer  frequency  resolution  is  desired,  the  input 
vectCM'  X  can  be  zero  padded  before  calling  the  function, 
wperigrmdx  zeros(300.1 )]); 

See  Also 

Iperigrm.  pkrtlime 

Reference 

[1]  >^^am  A.  Gardner,  Statistical  Spectral  Analysis.  A  Nonprobabilistic  Theory, 
pp.  5-7,  Prentice-Hall,  1988. 
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